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(I) 


IntTQductiQ-n 

Our  primary  objective  is  to  report  on  the  investigations  of  the  behavior  of  small- 
scale  gravity  waves;  their  effects  on  the  fluctuations  of  mesospheric  and  lower 
thermospheric  airglow  such  as  the  OH  infrared  bands  and  the  oxygen  green  lines. 
This  research  was  carried  out  under  the  ARC-AARC  airglow  modelling  program  and 
forms  part  of  the  MAPSTAR  program  which  includes  the  following  organizations:  The 
Brigham  Young  University,  the  University  of  Cincinnati,  the  University  of 
Southhampton  (England).  Utah  State  University,  the  University  of  Western  Ontario 
(Canada)  and  the  Phillips  Laboratories. 

The  theoretical  analysis  of  our  part  of  the  investigation  have  resulted  in  four 
papers,  two  of  which  have  already  appeared  in  the  Journal  of  Geophysical  Research, 
the  third  has  just  been  accepted  by  the  same  Journal,  while  a  fqurth  has  been  written 
up  and  is  now  undergoing  its  final  revisions  before  being  submitted,  again  to  the  same 
Journal.  The  experimental  part  of  our  program  has  also  resulted  in  several  papers 
which  have  been  submitted  and  accepted  by  Geophysical  Research  Letters. 

In  Section  (11),  we  present  two  papers  which  deal  with  small-scale  gravity 
waves  (G.W.)  and  their  effects  on  spatial-temporal  variations  in  airglow  structure. 

The  first  paper,  "Brunt-Ooppler  ducting  of  small-period  gravity  waves". 


1 


investigates  the  simultaneous  effects  of  two  possible  ducting  mechanisms  on  small- 
period  G.W.;  mechanisms  which  enhance  spatial-temporal  fluctuations.  The  first 
mechanism  is  the  ducting  produced  by  variations  of  the  Brunt  period  with  altitude. 
At  low  altitudes  (below  mesopause)  the  Brunt  period  is  of  the  order  of  5  minutes  and 
it  increases  to  10-15  minutes  in  the  thermosphere.  A  gravity  wave  with  period  less 
than  say  10  minutes  will  become  evanescent  above  a  certain  height  level.  The 
downward  reflection  from  such  height  levels  is  the  only  mechanism  in  the  atmosphere 
that  can  produce  fully  guided  modes.  Other  reflection  mechanisms  such  as  the 
inhomogeneous  background  density  and  temperature  variations  can  only  produce 
partially  guided  modes.  The  second  mechanism  is  the  downward  reflection  from 
horizontal  winds  which  Doppler  shifts  the  frequency  of  the  wave.  We  do  not 
consider  the  case  for  critical  layers  when  the  situation  is  far  more  complicated  in  this 
report.  The  combined  Brunt-Ooppler  ducting  can  produce  both  the  fully  and  the 
partially  guided  modes. 

The  dispersion  of  these  modes  with  respect  to  frequency  and  wind  propagation 
have  been  examined.  In  general,  the  frequencies  of  each  set  of  guided  modes 
increase  with  the  phase  velocity  in  a  complicated  way.  With  wind  dispersions,  if  the 
direction  of  the  wave  propagation  is  along  the  wind  direction,  the  wave  period  js 
Doppler-shifted  to  large  values  and  vice  versa.  This  is  expected  from  the  formula  for 
the  Doppler-shifted  frequency. 
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Whilst  the  effects  of  viscosity  have  been  negligible  over  nearly  all  the  guided 
modes,  the  effects  of  instability  are  not  negligible  for  most  of  the  modes.  In  fact,  only 
in  the  lowest  mode  (corresponding  to  lowest  eigenvalue  {  1/v^,),  or  the  highest 
horizontal  phase  velocity  v,),  i.e.  the  Lamb  mode,  can  the  instability  be  neglected. 

After  the  publication  of  this  paper,  we  have  found  an  error  in  our  wind  profiles, 
which  would  make  a  difference  for  the  higher  modes  but  would  make  no  difference 
for  the  lower  modes.  Since,  owing  to  instability  problems  with  the  higher  modes,  the 
results  would  in  any  case  be  suspect,  there  is  therefore  no  'jnod  reason  for  elaborate 
re-calculations  with  the  corrected  profile;  the  error  occurs  only  at  relatively  high 
altitudes. 

Simultaneous  observations  of  the  lo  .'*  attitude  ( <  100  km)  airglow  such  as  the 
557.7  nm  01  and  the  OH  bands  seem  to  reveal  considerably  more  short  period 
structures  (of  the  order  of  minutes)  than  at  higher  altitudes  (say  the  630.0  nm  red 
line).  In  this  part  of  the  report  we  have  shown  that  such  short-period  structures  can 
just  as  easily  be  produced  by  the  prevalence  of  short-period  ducted  wave  modes  at 
tower  altitudes  as  by  instabilities  from  larger  scdle  waves. 

The  second  paper  (Section  II)  that  deals  with  small-scale  spatial-temporal 
fluctuations  is  entitled  ”On  the  importance  of  the  purely  gravitationally  induced 
density,  pressure  and  temperature  variations  in  gravity  waves;  their  application  to 


airglow  observations”.  The  paper  shows  that  for  small-scale  gravity  waves  where 
v/c  <  <  1 ,  the  variations  in  pressure,  density  and  temperature  are  produced  primarily 
by  the  changes  in  the  background  gravitational  and  buoyancy  fbrces  as  the  air  parcel 
moves  up  and  down  in  altitude.  For  such  waves  the  "acoustic  component",  i.e.  the 
pressure,  density  and  temperature  variations  generated  from  the  pressure  gradient  in 
"purely  acoustic  wave  propagation"  play  a  very  small  role.  The  first  and  primary 
effect  will  be  hence  referred  to  as  G.I.C.  (gravitationally  induced  compression  and 
expansion). 

We  begin  by  showing  that  the  three  fluctuating  variables  (pressure,  density  and 
temperature)  can  be  explicitly  expressed  as  a  sum  of  two  terms,  one  describing  the 
G.I.C.,  and  the  other  the  "wave  motion”.  The  term  describing  G.I.C.  vanishes  if  we 
let  g  0  where  g  is  the  acceleration  due  to  gravity,  leaving  only  the  "wave  motion”. 

The  dependence  of  three  variables  on  v./c  are  then  plotted  for  different  G.W. 
periods.  The  results  show  that  for  small-scale  G.W.  (e.g.  v/c  <  18%)  the  GJ.C.  term 
by  itself  accounts  for  95%  of  the  temperature  and  density  fluctuations.  This  permits 
us  to  make  quick  estimates  of  say  height  variations  in  the  airglow  corresponding  to 
any  observed  temperature  variation.  For  instance,  a  10*  temperature  variation 
corresponds  to  a  height  variation  of  1  km;  the  simple  conversion  readily  establishes 
a  relation  between  small-scale  temperature  and  heig.ht  fluctuations. 


In  summary,  for  Section  II,  we  have  examined  mechanisms  such  as  ducting 
which  can  enhance  small-scale  temporal  fluctuations  in  airgfow  and  we  have  also 
separated  the  G.I.C.  from  acoustic  wave  compression  to  show  that  the  small-scale 
G.W.  is  mainly  dominated  by  G.I.C.  which  can  lead  to  the  simple  relationships  in 
small-scale  fluctuations  mentioned  above. 

In  Section  III  we  also  present  two  papers  which  deal  with  the  nonlinear 
integrated  and  local  response  to  a  strictly  linear  G.W.  Tnese  papers  are  important  to 
the  study  of  airglow  fluctuations  because  G.W.  behavior  has  often  been  inferred  from 
the  airglow  response  through  assuming  that  a  linear  G.W.  produces  a  linear  response. 
If  such  inferences  are  in  fact  unjustified,  we  neied  to  study  G.W,  by  means  other  than 
their  indirect  effect  on  sonr>e  minor  species  with  say  a  layered  background  structure. 
The  presence  of  additional  structures  in  the  airglow  produced  from  the  nonlinear 

response  should  add  to  our  general  study  of  spatial-temporal  airglow  fluctuations;  the 

'  '  ! 

fluctuations  in  this  case  are  produced  from  nonlinear  effects. 

We  were  led  to  the  first  work  in  Section  ill  by  the  papcir  of  Hines  &  Tarasick 
(Planets  Space  Sci.  851  (1987)]  who  showed  that  unlike  radar  and  lidar 
observations  which  measure  the  iocal  fluctuations  at  each  height  level,  the  ground- 
based  photorfieter  and  radiometer  measure  the  integrated  eirglow  that  can,  at  least 
in  first  order;  allow  the  large  nonlinear  local  Eulerian  response,  produced  from  the 
steep  density  gradient  of  the  bonom  side  of  the  minor  species,  to  be  trahsformed 
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away  with  a  Lagrangian  transformation.  Since  the  large  nonlinear  local  response  from 
the  bottom-side  of  the  minor  species  structure  has  often  been  blamed  for  the 
importance  of  the  higher  order  response  terms,  such  a  transformation,  which  can  be 
realized  in  the  integrated  airglow,  would  seem  to  indicate  two  things:  (1)  in  the 
integrated  airglow  the  higher-order  terms  are  filtered  out  so  that  a  linear  G.W.  can  only 
elicit  a  linear  response;  (2)  in  studying  G.W.  the  measurement  of  integrated  airglow 
would  therefore  enjoy  an  advantage  over  the  measurement  of  local  responses,  since 
the  former  filter  out  the  undesirable  nonlinear  effects,  whilst  the  latter  would  show 
such  effects  even  if  they  do  not  exist  in  the  linear  G.W.  itself. 

In  this  report  we  have  computed  the  higher  order  response  terms  in  the 
Lagrangian  frame  and  found  that  they  remain  too  large  to  be  ignored.  This  means  that 
the  higher  order  harmonics  should  exist  in  ti>9  response  even  if  the  G.W.  is  linear  and 
monochromatic  and  irrespective  of  whether  the  local  or  the  integrates  airglow  is 
observed. 

We  were,  however,  unable  to  determine  if  the  large  magnitude  of  the  higher 
order  terms  is  due  to  the  steep  vertical  density  gradient  just  mentioned.  Nevertheless, 
the  assumption  that  the  higher  order  terms  are  unimportant  and  that  a  linear  G.W.  can 
only  produce  a  linear  response  has  to  be  questioned. 

In  the  second  paper  in  Section  III,  we  have  developed  a  perturbation  expansion 
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for  computing  the  higher-order  response  terms  in  the  purely  Eulerian  system.  The 
results  again  show  that  the  integrated  higher-order  response  terms  are  important  and 
that  they  cannot  be  transformed  jway  either  by  a  Langrangian  transformation  or  by 
integ.'ation  by  parts.  In  this  case,  the  importance  of  the  higher-order  terms  appears 
to  be  connected  with  the  steep  density  gradient,  ever  thougii  such  a  connection 
carinot  be  established  in  the  Lagrangian  system. 

We  snould  mention  tnat  in  general  tha  gravity  wave  from  a  ground  source 
always  consists  of  a  linear  combination  o^  two  waves:  one  is  a  purely  travelling  wave 
along  both  the  vertical  and  horizcntal  directions,  while  the  other  is  stationary  along 
the  vertical  but  travelling  along  the  horizontal  direction.  The  linear  combination  takes 
into  account  both  the  purely  ttavelling  part  of  the  gravity  wave  and  the  partial  (which 
may  sometimes  be  very  small  and  sometimes  be  very  large)  reflection.  The  purely 
travelling  part  can  only  be  computed  in  the  usual  first  order  response  without 
encountering  the  problem  of  secularity  which  requires  renormalization  techniques  of 
the  K6M  (Krylov-Bogolinbov-Mitropolsky),  (See  also  "Averaging  Mathods  in  Non-Linear 
Dynamical  Systems,  Appl.  Math.  Sci.  (Springer- Veriag)  by  Sanders  and  Verhulst).  For 
the  present  problem,  we  shall  consider  only  the  stationary  vertical  propagation  which 
does  not  encounter  the  secularity  problem  in  any  order  of  response. 

Our  overall  theoretical  conclusion  for  this  part  of  the  report  is:  (1)  so  long  as 
a  minor  species  is  used  for  determining  gravity-wave  behavior,  the  nonlinear  response 


is  always  important  for  an  average  wave  amplitude  whether  one  uses  photometer,  a 
riometer,  a  lidar  or  radar  as  probes.  In  this  case  only  from  direct  measurements  oi 
some  major  species  in  hydrostatic  equilibrium  can  one  expect  the  behavior  of  the 
response  to  strictly  reflect  the  G.W.  behavior;  i2)  in  measuring  the  local  minor 
species  response  (with  say  either  the  radar  or  lidar),  the  density  Gradient  at  the 
bottom-side  of  the  minor-species  layer  plays  a  crucial  role  in  both  the  linear  and  the 
nonlinear  response  terms,  since  the  gradient  cannot  be  transformed  away  by  an 
integration  by  parts  as  it  can  be  for  the  first  order  integrated  airg  jw;  |3)  for  the 
integrated  airglow  we  are  not  able  to  determine  if  the  large  nonlinear  response  in  the 
Lagrangian  formulation  is  caused  by  the  steep  gradient,  since  a  nonlinear 
transformation  is  required  to  go  from  the  Eulerian  to  the  Lagrangian  system.  Further 
work  is  needed  to  clarify  this  point. 

In  Section  IV  we  present  the  imaging  of  the  OH  airglow  wave  structure 
measured  in  late  March  and  early  April  1990  at  Maui,  Hawaii.  A  second  set  of 
measurements  on  April  18-28,  1990  were  taken  at  Haleakala  following  the  success 
of  the  first.  The  results  have  been  submitted  to  Geophys.  Res.  Lett,  and  preprints  of 
these  publications  are  presented  in  Section  IV.  The  principal  investigators  for  the 
MAPSTAR  portion  of  the  ALOHA  1990  campaign  are  Drs.  M.  J.  Taylor  (Univ.  of 
Southhampton)  and  R.  P.  Lowe  (Univ  of  Western  Ontario). 
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SMALL-SCALE  GRAVITY  WAVES 


PART  A:  BRUNT  DOPPLER  DUCTING  OF  SMALL  PERIOD  GRAVITY  WAVES 


D.  Y.  Wang 
and 

T.  F.  Tuan 


Physics  Department,  University  of  Cincinnati 
Cincinnati,  Ohio  45221 


'Permanent  Address:  Institute  of  Geophysics,  Academia  Sinica,  Beijing,  PRC 


ABSTRACT 


The  variation  of  the  Brunt  period  with  height  lenc*s  itself  to  a 
natural  ducting  and  filtering  mechanism  for  low  altitude  short-period  , 
gravity  waves.  We  investigate  this  mechanism  in  combination  with  Dcppler- 
ducting  produced  by  the  variation  in  horizontal  winds.  Both  the  frequency 
dispersion  at  fixed  propagation  direction  and  the  direction  dispersion  at 
fixed  frequency  have  been  examined  in  a  OOSPAR  badoground  atno^iieie  with 
zcnal  and  meridional  winds. 

(Xir  results  shew  that  indeed  the  lew  altitude  short-period  gravity 
waves  cure  not  only  ducted,  but  unlike  the  usual  ducting  mechanisms  due  to 
uneven  structure  and  dissipation  vAiich  cxily  produce  the  partially  guided 
modes,  this  mechanism  produces  primarily  guided  modes  in  the  absence  of 
winds  and  a  mixture  of  fully  and  peurticilly  guided  modes  with  the  winds. 

The  wind  effects  eure  very  large  on  the  higher  and  less  significant  on 

the  few  lowest  modes,  including  the  lamb  mode. 

Investigations  of  viscous  dissipation,  non-linearity  and  instability 
have  shown  ttet  viscosity  is  uninportant  for  most  altitudes  of  interest  auid 
that  non-linearity  and  instability  can  play  a  role  for  all  but  the  lowest 
guided  modes. 

We  propose  that  siaultaneous  continuous  observation  of  airglow  at 
mesoqjheric  and  ionospheric  altitudes  be  nede  to  verify  not  only  the  low- 
altitude  Brunt-f  c^ler  ducting  for  short-period  gravity  waves,  but  the 
vertical  energy  distributicxi  of  the  medium  and  large-scale  TID's. 


1.  nnmxjcncN 


•niGre  has  been  cx)nsiderable  interest  in  the  small  to  medium  scale 
atmospheric  fluctuations  with  periods  in  the  rainge  of  a  few  minutes  to  lew 
tens  of  minutes  and  horizontal  phase  velocities  in  the  range  of  lew  tens  to 
lew  hundreds  of  meters  per  seexxid. .  They  have  been  detected  through  nany 
different  teehniques  including,  for  instance,  the  mesospheric  and  Icwer 
thermospheric  airglow  observations  (e.g.  Okuda,  1962;  Silverman,  1962; 
Krassovsky  and  Shagaev,  1974a, b,  and  1977;  Myreibo  et  ad.,.  1983,  1984; 
Peterson,  1979;  Freund  and  Jacka,  1979;  Hapgood  and  Taylor,  1982;  Peterson 
and  Adams,  1983;  Clairanidi  et  al.,  1985;  Jacob. and  Jacka,  1987;  etc.),  or 
the  different  radar  systems  (e.g.  Manson  and  Meek,  1980;  Kansen  et  al., 
1981;  Smith  and  Ftitts,  1983;  Vincent  and  Reid,  1983;  Meek  et  al. ,  1985). 
The  small-scale  fluctuations  hen/e  senetiaes  been  interpreted  as 
instabilities  (Hodges,  1967;  Than  et  al.  1979;  Schoeberl  et  al.,  1983; 
KLostermeyer  and  Ruster,  1981,  1984). 

In  this  paper,  insteeK)  of  searching  for  different  possible  meefetnisns 
that  produce  these  scall-scade  fluctuations,  we  shall  make  a  systematic 
investigatiian  of  short-periad  gravity  waves  %#ith  periods  in  the  range  of 
the  yzuriation  of  the  Brunt  period  vrith  height,  a  variation  «hidh  is  in  the 
zfpnadjnata  range  of  the  ^^ically  observed  periods  just  mentioned.  We 
will  focus  cur  attention  primarily  on  the  small-period  gravity  waves  with 
higher  phase  velocities. 

The  purpose  for  maJdng  a  study  of  the  short-period  linear  gravity 
waves  are  as  follows; 

(1)  Pswer  !^»ctrum  analyses  of  simultauiecus  photometric  observations 
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of  mesospheric  auxi  lower  thermospheric  optical  esnissions  such  as  the  557.7 
nm  01,  the  730  rm  and  790  m  OH  band  as  well  as  the  1.53  m  OH  radiaaetric 
data  (taken  by  Peterson  and  Pendleton  in  the  1984  MAPSTAR  Campaign)  serins 
to  shew  a  surprisingly  large  number  of  power  peaks  with  short  periods. 

This  would  seem  to  suggest  that  the  higher  Fourier  components  as  well  as 


the  short  paeriod  gravity  waves  play  an  important  role  in  saall-scale 
at3no^]heric  fluctuations. 

(2)  Sijtultaneous  observations  of  557.7  m  and  630  m  01  fay  Weinberg 
et  al.  (1973)  would  seem  to  show  that  considerably  more  ^K>rt  period 


structures  occur  in  the  557.7  m  at  lower  altitudes  rath^  than  the  630  m 
in  the  iono^heric  altitudes.  Vliilst  other  eigilanations  are  possible,  ve 
will  show  that  this  is  consistent  with  the  behavior  of  the  sijiple  short- 
period  lineau:  gravity  wave  theory, 

(3)  The  theoretical  justification  for  our  investigations  is  based  on 
the  fact  that  so  faur  the  treatiaent  of  gravity  wave  ducting  has  been 
confined  to  three  paurticular  mechanians:  the  variation  in  atmospheric 
structure  (e.^.  Fjrancis,  1975  etc.);  the  vauriation  in  dissipation 
(Richmond,  1978;  Yu  et  al.,  1980;  lUan  and  Tadic,  2982);  and  the  variation 
in  faackgreund  winds  (Qiiacnas  and  Hines,  1986) .  Ha  have  fiawd  there 
exists  anothej:  type  of  ducting  machanuai  valid  only  for  gravity  waves,  or 
higher  Fourie^  ocaponents  of  gravity  waves  with  periods  in  the  range  of.  the 
Brunt  period  (about  4.5  min  at  low  altitudes  to  15  min  at  highar 
altitudes) .  '[his  variatior  in  the  Brunt  period  with  height  allows  a  low 
altitude  gravilty  vave  with  initial  Fourier  period  oceponent  greater  than 
its  local  Bru]{it  period  to  propagate  upwards  to  sane  height  level,  Zq,  where 
it  becomes  egiiial  to  the  local  Brunt  period  (Fig.  1) .  Ihe  resulting 
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(downward  reflecticjr.  together  with  the  up«/ard  reflection  frcm  the  rigid 
ground  would  produce  a  series  of  g»iided  modes  for  each  such  Fourier 
frequaicy  ocnponent,  leaving  only  the  much  Icwer  frequency  conponents  to 
continue  their  upward  propagation.  Thus,  the  variation  in  the  Brunt 
frequency  curve  could  serve  as  a  ducting-filtering  mechardsm  (Wan^,  and 
TUan,  1985,  1986)  which  filters  out  and  restricts  the  hi^ier  frequency 
ccnpcnents  to  only  the  lower  altitudes  (from  new  on  to  be  referred  to  as 
Brunt  ducting) . ' 

Whilst  the  eKxustic  iiiave  branch  can  readily  propagate  at  higher 
frequencies  above  2  “  zq,  there  is  no  way,  of  course,  that  a  gravity  yave 
C2in  turn  into  an  acoustic  %«ave,  since  the  latter  requires  a  horizontal 
phase  velocity  greater  than  the  local  speed  of  sound.  In  a  horizontally 
stratified  ataosphere,  the  horizontal  phase  velocity  cam  not  vauy  with 
heic^.  Thus,  for  z  >  zq  the  short-period  gravity  waves  beocne  evanescent. 
Tha  further  rt  guirement  of  finite  energy  source  would  limit  the  hi^ 
altitude  solutions  to  those  which  are  "exponentially  decaying"  emd  which  in 
turn  limit  the  propagation  modes  to  only  the  fully  guided  inodes. 

DacaMHW  the  horizontal  winds  will  zULways  have  an  inportant  effect 
above  80  Jon,  va  %dLll  in  this  p^ier  ocobire  the  Doppler  ducting  of  Chinonas 
and  Hines  (1986)  and  the  Brunt  ducting  (Wang  and  TUan,  1985  aid  1986)  to 
produce  a  simultaneous  ducting  nachaniaB  which  can  produce  both  the  fully 
and  the  partially  guided  modes  and- which  will  be  referred  to  as  die  Brtmt- 
Dcpplar  ducting  mechanism  for  ahort  (See  also  wang  and  TUan,  1987) .  The 
zonal  and  meridioiwd  wind  profiles  of  lUrlses  and  Gillette  (1982)  based  on 
Lie  beKdaground  mean  circulation  wind  profiles  of  Lindzen  and  Hong  (1974) , 
and  Roble  et  al.  (1977)  «dll  be  used  for  die  horizontal  wind  profiles.  The 
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OOSPAR  nodel  will  be  used  for  general  background  atstosphere . . 

In  Section  2  we  shall  develop  cind  describe  the  fundamental  equations, 
the  boundary  conditions  and  the  methods  of  solution.  In  Section  3  we  shall 
present  the  solutions  for  the  case  when  the  wind  is  absent  and  the  guidance 
is  provided  by  Brunt  ducting  alone.  The  dispersion  curves  for  the  fully 
guided  modes  at  different  gravity  wave  periods  and  the  kinetic  energy 
density  will  be  given  as  functions  of  2,  the  altitude,  (Fig.  4-a  and  b)  for 
the  various  discrete  fully  guided  modes. 

In  Section  4,  we  include  the  wind  profiles  (Fig.  5-a  amd  b)  and  will 
disojss  the  effects  of  the  zonal  auid  meridional  winds  on  the  different 
directions  of  wave  propagation.  Dispersion  curves  as  a  function  of  periods 
for  different  propagation  directions  will  be  given  and  oonpared  with  the 
di.«yersion  carves  for  the  windless  case  (Fig.  8-a,b  and  Fig.  9-:a,b) .  Vte 
also  shew  the  di^iersion  curves  as  a  function  of  propeigatim  direction  at 
fixed  period  (Fig.  10-a,b) . 

In  Section  5,  we  consider  the  effects  of  stability  and  dissipation  of 
the  different  propagation  modes,  the  range  of  vertical  and  horizontal  wave 
lengths  for  vhich  neither  effect  is  significant.  In  Section  6  we  present 
the  conclusion. 


2.  RUnftKEJnM.  HQUATICKS 


The  differential  equations  gcsveminq  the  propagation  of  gravity  waves 
in  a  horizontally  stratified  atmosphere  whose  tenperature  and  horizontal 
wind  velocity  vary  in  an  artoitrary  manner  with  height  have  teen  derived  by 
Pitteway  and  Hines  (1965) .  Choosing  the  pressure  perturbation  ap  cuid  the 
vertical  velocity  perturbation  aw  eis  dependent  variables,  the  ^stera  of 
equations  can*  be  transformed  into  two  coupled  first-order  equations  given 
by 

“  i  _!L  (wb^  -  o2)  (AW) 
dz  0 

(i: 

t  («»  +  !^  1  (AW)  »  i  JL  ( -  -i- 

dz  QdZ  w  ^ 

0 


where  o  is  the  OappleiMahifted  frequency  that  would  be  naasured  in  a 
coordinate  system  noving  with  the  atiaosphere.  It  is  a  ftmcticn  of  height 
through  the  ptreeenoe  of  horlzcntal  wind  ooBponent,  Vq^Cz)  .  w  is  the  wave 
ftacpanq^  in  a  fixed  obbxdinate  systen.  They  are  related  by 


0  -  «  -  k^Vox 


(2) 


Itie  veuricibles  tv  ani  f)  defined  by 


(3) 

(4) 

(5) 

(6) 


Ihe  other  quantities  have  their  usual  meaning. 

Eq.  (1)  is  subject  to  the  following  jassinptians  and  limitations:  (1) 
the  effects  of  viscosity,  themal  csnduction  auid  souroe  terns  are 
neglected;  (2)  the  aabient  atnosphere  is  horizontall.v  ’'^ratified,  all 
unperturbed  quantities  have  z-dependenoe  only;  ant'  ['■}  xhe  perturbations 
are  small  so  that  the  nonlirteeu:  terns  can  be  emitted,  and  all  perturbation 
quantities  can  be  Fourier-deoenposed  into  ncnochronatic  plane  %«vee  given 
by  - 


f  (X,  z,  t)  «  f  (2)  exp  f  i  («t  -  JtxX)  1  (7) 

The  helg^-varying  beebground  %finds  are  aneinwd  to  be  horizontal, 
in(;d.uding  zonal  and  meridional  oonponents,  Vqj  and  Vom,  respectively. 

Since  the  vertical  wind  is  generally  much  less  than  the  horizontal  winds, 
it  vaxild  seem  reasonable  to  sinplify  the  discussion  by  neglecting  the 


vertical  wind.  The  fact  that  the  wi^  perpcndiculcu:  to  the  wave  vector  has 
no  effects  at  all  on  the  wave  propagation  allows  us  to  consider  only  the 
wind  along  ky,  the  direction  of  wave  propagation.  Thus,  from  Fig.  2,  we 
may  write 

Vox  =  V02  oosor  +  Von,  sin*  (8) 

where  a  is  the'  angle  between  ky  and  the  East  dirbcticn. 

There  nay  be  a  critical  height  level  v4iere  the  horizontal  phase  speed 
equals  the  wind  speed,  and  (l  >=  0.  The  solutions  of  the  Bq.  (1)  are 
singular  (Lin,  1955)  due  to  the  fact  that,  for  harmonic  wave  trains,  the 
inviscid  equations  intrxsducse  infinite  ^tears  at  the  critical  level.  This 
inplies  that  in  the  neighboihood  of  such  a  level  viscous  focoe  cannot  be 
ignored.  The  singularity  may  also  be  circuavented  by  dealing  with 
transient  wove  packets  (Tblstoy,  1973) .  Such  a  theory  renains  to  be  fully 
developed,  and  attention  at  the  present  time  is  confined  to  those  modes 
which  do  ix3t  enocunter  a  singular  level  at  any  heicht. 

Tb  solve  ths  equations  appropriate  bounary  conditions  are  required. 

At  the  rigid  growid  bouidazy  the  vertlcid  velocity  muet  vanish.  In  the 
upper  half  apaoa  taapazature  and  wind  can  ba  viawad  as  oca»tant  (Fig.  1  and 
5-a,b) ,  the  energy  propagation  has  to  ba  only  iqwezd  at  Thus,  the 
radiation  condition  must  be  satisfied  at  infinity  (Rsddy,  1969;  Fansim, 
1974) .  With  the  help  of  Bq.  (1)  these  conditions  can  ba  eaqaesaad  as 

[  “  »h^3z-o  ■  0 

dz 
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0(2^)  o  exp(i3c„z) 


(10) 


wheni  ka  is  the  vertical  wave  nunber  for  the  upper  half  space  with  constant 
tenperature  and  wind. 

For  a  given  atzaaspherlc  model  Dq.  (1)  can  be  runerically  solved  by 
full  wave  oonputations  under  the  boanii  -y  conditions  of  Bcjs.  (9)  and  (10) 
for  the  range  of  horizontal  phzise  velocities  we  consider  in  this  paper. 


3.  WINIX£SS  HXXS  IN  'HIE  OOfiPAR  AlMCSniERE 

In  order  to  understand  tha  Brunt  duc±ing  mechcUiism  for  small-period 
gravity  waves  in  a  realistic  atmo^^here,  we  first  consider  an  atmospheric 
itcdel  without  wind.  The  1972  OOSPAR  teirperature  profile  with  an  exospheric 
tenperature  of  1000  is  used  in  the  present  calculation.  The  upper 
boundcuy,  above  vhich  the  atmosphere  is  assumed  to  be  uniform,  is  fixed  at 
400  kn. 

The  phase  velocity  di^sersion  curves  for  the  atmosphere  aure  shown  in 
Fig.  3.  Each  curve  represents  a  mode  labelled  by  the  node  number  in 
parenthesis  on  the  right-hand  side  of  Fig.  3.  The  discrete  inodes  indicated 
by  short  dash  lines  can  be  viewed  as  bound  states.  Sinoe,  in  general,  the 
eigenvalues  are  proportiensa  to  inverse  horizontal  phase  velocity  (Yu  et 
ad.,  1980) ,  we  plot  these  discrete  mextes  against  decreasing  velocity.  The 
noite  with  the  same  nodes  are  connected  by  solid  lines,  and  the  lanb  modes 
as  pseudomodes  (Press  and  Hauicrider,  19o2;  and  Pfeffer  and  Zarichny,  19(53) 
are  connected  by  dotted  lines.  The  Lanb  nodes  have  horizontal  phase 
velocity  of  about  316  ly^sec,  and  are  nearly  non-dispersive.  The  nodes  with 
la«r  phaM  velixi.ties  aure  ducted  between  the  ground  and  the  lower 
themd^diere,  and  can  beexmw  highly  dLqpezsive  at  the  low  velod^  end. 

Fig.  4-a  gives  the  first  few  solutions  of  wave  function,  for  wove 

period  T  >■  6  ein.  The  nods  nuaber  assigned  to  each  mla  lo  to  node 
ruber  of  the  node,  including  the  infinity  as  a  nods.  It  can  be  seen  that 
decays  to  zero  very  quickly  above  1^  km,  and  the  weves  are  guided  by  the 
ground  and  the  thermospheric  tenperature  gradient. 


The  time-averaged  kinetic  energy  per  unit  volume,  associated  %ri.th  the 
wave  percurbarion,  is  given  by 


|au|2  +  |aw|2  ] 


(ll> 


v4iere  au  is  the  horizontal  velocity  perturbation.  The  energy  densities  for 
all  ducted  modes  with  periods  of  4.5  -  14  min  are  concentrated  at  low 
altitudes  either  neeur  the  ground  or  near  the  lower  thermosphere.  Fig.  4-b 
shows  the  results  for  T  *  6  min. 


4.  FlOEftavnON  MXES  IN  WINOT  AMOSFHERE 


The  zcwil  and  meridional  wind  profiles  used  in  the  present  stufy  are 
shown  by  the  bold  lines  in  Fig.  5-a  <ind  b.  They  cure  given  for  the  winter 
solstice  at  mid-latitude  (42*  N)  and  at  local  midni^t.  The  East  is  taken 
as  positive  for  the  zoral  wind  conpsnent,  while  the  North  is  taJcen  as 
positive  for  the  meridional  wind  ocnfxxiait.  These  wind  ocnponents  cune 
tcOcen  from  Porbers  and  Gillette  (1982)  v^vsse  calculations  are  in  turn  based 
on  the  mean  circulaticn  conponents  (giv^  by  the  thin  line  in  Fig.  5-a  and 
b)  prcviided  by  Lindzen  and  Hong  (1974)  and  Roble  et  al.  (1977).  These  wind 
ccnporients  may  be  taken  as  time-independent  due  to  the  rather  short-period 
%aves  we  ocnsider  in  ccoisarisan  with  the  tidal  periods.  It  is  scaieMhat 
unfortunate  that  the  region  of  greatest  iaportanoe  (80  -  120  loa)  in  the 
present  paper  is  edso  a  region  of  some  laioertainly  in  the  wind  profiles  due 
to  boundary  effects.  Some  ?qodification  of  the  winds  was  performed  to  allow 
for  a  smooth  merging  with  profile  below  100  km.  The  new  "dRA"  model  (MAP 
Haoidbook  (1985) )  would  help.  The  temperature  profile  is  the  same  as  that 
described  in  Section  3. 

For  the  windy  ateosphere,  model  solutions  have  bean  oceputed  fbr  wave 
periods  at  T  ■  4.5,  5.5,  6,  8,  10,  127  14  sin  at  propagation  direction 
angle  a  -  0* ,  30* ,  45* ,  60* ,  90* ,  120* ,  135* ,  150* ,  180» ,  210* ,  225* ,  240* , 
270*,  300*,  315*  and  330*.  The  preecnoc  of  wind  and  the  resulting  Doppler 
frequency  shift  (Bq.  2)  can  sufficiently  alter  the  beciagraund  atmoephere  to 
allow  seme  fUlly  guided  modes  to  change  into  pertially  guided  modes.  We 
will  specifically  consider  two  cases  to  show  how  the  effect  of  the  Winds 
can  change  fully  ducted  into  partially  ducted  modes.  The  waximLW  wind 
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speeds  ii  the  wind  profiles  of  Fig.  5-a  and  b  would  restrict  our  presott 
treatment  to  gravity  waves  with  horizontal  phase  velocity  greater  than  150 
ny'sec  to  avoid  the  singulaurities  generated  by  critical  layers  already 
mentic»ied  in  Section  2. 

Fig.  6-a,b  and  Fig.  7-a,b  show  the  wave  function  ^  and  the  kinetic 
energy  density  for  the  6  min  modes  at  a  =  0*  (Eastward  prcpagaticn)  and  90* 
(Northweird  propagation) ,  respectively.  As  cam  be  seen,  the  presencje  of  the 
wind  system  greatly  modifies  the  properties  of  wave  guidance.  Some  inodes 
with  lower  horizontal  pha^  velocities  become  peurtiadly  ducted  (Figs.  6-a 
and  b) .  Ihe  energy  densities  of  the  fully  ducted  modes  are  concentrated  in 
the  region  below  120  km,  much  lower  than  that  for  the  landless  nodes.  Hiis 
is  caused  by  the  rapid  increase  in  the  wind  system  with  height  iduch  can 
substantially  reflect  wave  energy  downwards  to  lower  adtitudes.  The 
horizontal  phase  ^leeds  of  the  wave  modes  (as  measured  in  a  fixed 
ooorcUnate  system)  are  adso  much  altered  by  the  presence  of  the  wind 
profiles.  The  effect  of  wind  structure  is  quite  significant  even  for  the 
non-dispersive  Lairb  modes  with  a  high  horizontal  phase  velocity  Vptgf. 

The  dispersion  curves  for  the  variation  in  the  modal  horizontal  phase 
velocity  with  different  wnve  periods  at  each  propagation  dizectlan,  and  the 
curves  for  the  variation  in  modal  phase  velocity  with  different  propegatich 
directions  at  each  fixed  period  are  ooitpiled  but  will  not  all  be  shown  due 
to  lack  of  spaoe.  Here  we  show  the  dispersion  curves  for  a  *  0* ,  180* 
(Figs.  8'*a,b),  90*,  270*  (Figs.  9^,b),  and  the  directional  variation 
curves  of  phase  velocities  for  T  »  6  euid  14  min  (Figs.  10<-a,b) .  The  other 
curves  have  chauracteristics  siinil^u:  to  these  more  extreme  cases.  In 
eKldition  to  the  notations  and  symbols  defined  for  Fig.  3  in  Section  3,  a 
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sign  is  used  to  indicate  partially  ducted  iwades. 

The  dispersion  curves  in  Fig.  8-a  and  b  shew  the  ducted  nodes  when  the 
propagation  direction  is  alcsTg  the  direction  of  the  zonal  wind  (q=o*  ) ,  and 
against  the  zcxial  wind  (q  =  180*).  Coipared  to  Fig.  3,  the  calLy  possible 
similarity  is  the  b^avior  of  the  Lamb  node  v^ch  continues  to  be 
relatively  non-diqpersive,  although  the  phase  speed  is  shifted  fran  316 
To/sec  to  337  ii^sec  for  a  =  0*  and  296  i^sec  for  a  =  180* .  This  is  to  be 
expected,  since  the  zonal  wind  has  a  more  or  less  const^mt  magnitude  of 
about  20  xa/s^  from  the  ground  ip  to  about  30  km,  a  height  range  where  a 
significant  portion  of  the  Lamb  mode  energy  is  located.  In  other  re^)ects, 
Fig.  8-a  looks  quite  different  from  Fig.  3.  CVing  to  Doppler  wind  shift, 
no  fully  ducted  modes  can  be  found  for  waves  with  periods  greater  than  8 
min.  Fig.  8-b  bears  some  resemblance  to  Fig.  3,  and  shews  the  existence  of 
fully  ducted  modes  for  waves  with  periods  ip  to  14  min.  The  physics  b^rind 
it  is  that  for  a  »  0’  the  waves,  propegating  ed.ways  in  the  direertion  of  the 
background  wind,  can  be  Doppler  shifted  by  the  wind  to  longer  periods  which 
leads  to  reducing  the  reflection,  and  allowing  greater  energy  leakage 
ipwards  into  higher  altitudes.  Thus,  the  waves  with  periods  greater  than  8 
min  are  no  longer  ftilly  ducted.  At  the  ^xsrt  period  end  of  wove  spectrum, 
however,  sene  waves  with  periods  less  than  4.5  min  (the  lower  limit  of  the 
Brunt  period)  nay  be  shifted  to  form  fully  ducted  gravity  wave  modes.  For 
a  »  180*  the  wave  propagation  is  opposite  to  the  wind.  The  wave  period  is 
shifted  to  the  shorter  end,  and  the  wave  reflection  and  ducting  become 
»ih2moed.  In  this  case,  even  waves  with  period  greater  than  14  min  (the 
ipper  l^t  of  the  Brunt  period)  caut  be  turned  into  fully  guided  modes. 

The  dispersion  curves  in  Fig.  9-a  and  b  zure  due  to  pur^y  meridional  , 
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wind.  They  bear  more  simileority  to  Fiq.  3.  Ihe  phase  apoods  for  the  Iamb 
modes  eure  nearly  the  same  as  that  in  the  windless  atmo^here,  ed.thauc^ 
quite  significant  differences  can  appear  for  other  modes  with  Icwer  phase 
velocities.  We  point  out  that  vhile  the  meridional  wind  in  the  region 
curound  100  km  has  anplitudes  as  leurge  as  zoned  wind  (Fig.  5-a,b),  its 
effect  ai  the  wave  guidance  arB  much  weaker.  This  can  be  understood  in 
terms  of  the  general  pattern  of  meridional  wind  thich  reverses  itself  at 
different  height  levels  in  a  synmetric  way.  The  wind  reversals  cancel  cut 
much  of  the  effects  of  Doppler  shifting  that  the  more  mcno-directional 
zonal  wind  can  prtxiuoe.  Nevertheless,  the  ceuioellation  is  not  ccoplete, 
and  the  effects  of  meridional  wind  can  not  be  ignored  at  all,  particularly 
vhere  lower  phase-velocity  waves  are  concerned. 

Fig.  10-a  and  b  show  the  directional  dispersicn  of  the  guided  modes  at 
fixed  periods  (6  min  cind  14  min,  respectively) .  Tlie  relative  ^mnetry  of 
the  modes  about  a  =  180*  is  due  again  to  the  meridiunal  wind  pattern.  The 
rise  of  the  modes  eind  hence  a  lowering  of  the  phase  velocity  as  a  increases 
from  zero  in  Fig.  10-a  can  again  be  understood  in  terms  of  the  dispersicn 
relaticxi  which  provides  for  lower  horizontal  phase  velocity  as  Q  increeises 
throuch  Doppler  shift,  bnee  again,  the  chzuige  in  phase  velocity  is  the 
least  for  the  Lamb  mode  which  is  sufficiently  oaperitrated  at  lower 
altitudes  where  the  veuriation  with  hei^t  of  either  the  zonal  or  the 
meridionad  wind  is  relatively  small.  The  preponderance  of  modes  near  a  > 
180*  for  the  14  min  wave  (Fig.  10-b)  and  near  a  «  0*  or  360*  for  the  6  min 
wave  can  again  be  tinderstood  in  terms  of  Doppler  wind  shift.  The  longer 
period  modes  cam  only  be  Doppler  wind  shift.  The  longer  period  modes  can 
cxdy  be  Doppler-shifted  to  the  higher  frequency  end  auid  hence  require  winds 
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at  a  =  180* ,  whilst  the  shorter-period  nodes  Cein  be  more  readily  Doppler- 
shifted  to  the  lower  frequency  end  and  hence  require  a  =  0*  or  360* . 
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5.  siABnxnr  ce  nm  hkve  modes 


All  dissipations  and  nonlinear  effects  have  been  ignored  so  far.  To 
determine  the  ocnditions  for  which  a  linear,  inviscid  wave  is  a  good 
approximation,  and  the  cxmditians  for  which  nonlinearity  and  dissipations 
would  have  to  be  included,  the  inertial,  nonlinear,  and  viscous  terns  have 
been  calculated  in  terns  of  the  linear,  inviscid  modal  soluticns.  Ihe  data 
of  viscosity  coefficient  cure  taken  from  the  model  used  by  Midgley  and 
Licmohn  (1966) .  The  Richardson's  nutber  given  by 

(12) 


has  been  ccrputed  to  investigate  the  possibility  for  instabilities  to 
occur.  In  all  our  calculations  the  gratvity  waves  have  been  normalized  to 
density  fluctuation  of  20%,  15%,  and  10%  at  a  height  of  100  km. 

In  the  windless  atmosphere  for  guided  modes  with  period  less  than  9 
min  the  nonlinear  term  is  much  less  than  the  inertial  tern  at  lower 
ailtitud^  as  expected.  The  same  is  true  at  higher  altitudes,  where  the 
waves  have  purely  imaginary  vertical  wave  numbers  with  values  greater  than 
1/2H  (H  is  sc2LLe  height) ,  and  decay  very  quickly  so  that  their  amplitudes 
are  limited  to  small  values.  The  nonlinear  term,  however,  can  be 
corparable  to  the  inertial  term  cuound  120  -  140  km.  The  viscous  term 
below  240  km  is  much  less  than  the  nonlinear  term,  and  can  be  neglected. 
Viscosity  beccmes  uiportant  at  very  high  altitudes,  but  does  not  affect  the 
guided  modes  too  nwch,  since  their  wave  energy  is  confined  to  lower 
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altitudes.  Ihe  Richeuxlsan's  numbers  cf  these  wave  inodes  are  cilways  much 
greater  than  one  qu£u±er  at  edl  altitudes.  Ihus  there  is  no  cheuxe  for 
instability  to  ocxair.  Fig.  11-a  and  b  show  the  three  terms  and  the 
Richaurdson's  number  for  the  5.5  min  windless  mode  with  phase  velocity  of  72 
nv'sec  at  density  variation  of  20%. 

For  the  inodes  with  periods  greater  than  9  min  in  the  windfree  case, 
the  nonlinear  term  becomes  much  larger  than  inertial  term  above  100  ]an, 
quite  different  frcm  the  previous  case,  since  the  wave  my  not  decay  fast 
enouc^  to  balance  the  cagxmmntial  growth  due  to  the  decrease  of  the 
background  densi^.  The  viscous  term  which  is  £d.ways  less  than  the 
nonlinear  term  on  still  be  ignored,  but  becxxnes  conparable  to  the  inertial 
term  at  hi^ier  <iltitudes.  The  Richardson’s  numbers  at  some  height  levels 
are  less  than  one  quarter,  and  instability  can  arise  at  these  levels. 

The  presence  of  wind  will  change  the  pheuse  velcsity  of  the  modes  and, 
in  consequence,  greatly  affect  the  stability  of  such  modes.  Calculations 
show  that  only  the  Iamb  inodes  of  4.5  -  14  min  period  propeigating  in  any 
direction  can  still  be  viewed  as  stable  in  the  presence  of  winds. 

Stability  of  other  inodes  with  lower  phase  speeds  strongly  depends  on  wave 
period,  horizontal  phase  v<elacd.ty,  and  propagation  direction.  For  a  >  0* 
only  waves  %d.th  periods  less  than  6  min  can  have  one  or  mote  other  stable 
inodes  with  lower  phase  speed.  The  nusbers  of  stable  modes  increase  as  wave 
period  decreases.  The  same  thing  happens  at  periods  less  than  12  min  for 
a  »  180* .  In  the  case  of  purely  meridional  wind,  each  period  within 
.4.5  -  14  min  can  have  one  or  more  stable  inodes  other  than  the  Iamb  mode. 

The  lower  limit  of  the  horizonted  phase  velocity  for  stable  modes,  roughly 
speaking,  is  around  180  -  200  m/sec.  Figs.  12-*a,b  and  13-a#b  display  the 


three  terns  and  the  Rict^iaidson's  nuntexs  foL  5.5  nin  wave  with  pha^  speed 
of  260  nv^sec  (a  =  0* ) ,  amd  12  rain  wave  with  pheise  speed  210  Vsec 
(a  =  180*)  reflectively.  Those  ncdes  aure  stable  even  for  the  relative 
density  fluctuation  of  20%  at  100  Icra. 

In  ooi:>clusicn,  in  the  real  ataofihere  vrith  arbitrarily  vaurying  wind 
system,  cne  would  expect  that  only  waves  with  periods  less  than  6  min  and 
horizontal  phase  velocity  greater  than  200  ny'sec  can  occur  more  frequently 
as  stable,  fully  ducted  nodes.  These  are  consistent  with  the  experinentail 
results  mentioned  in  Introduction  of  this  paper.  More  detailed  discussion 
about  isplicaticn  of  the  theory  to  experiments  will  be  presented  in  other 


6.  OCNdOSION 


We  haive  found  and  investigated  a  duc±ing  medwiisah,  the  Brunt  IXxrting, 
valid  only  for  gravity  waves  with  periods  in  the  range  of  the  Brunt  period. 
UhliJce  structural  ducting  or  dissipative  ducting,  the  Brunt  ducting  can  by 
itself  produce  fUllv  guided  mooes,  ihe  only  other  mechanism  that  can  do  so 
is  the  Doppler  ducting  when  the  wind  velocity  along  the  direction  of  v^ve 
propeigatian  exceeds  the  pheise  velocity.  The  combined  Brunt-Doppler  ducting 
can  produce  both  the  fully  and  the  partially  guided  modes.  We  haive 
runericedly  oonputed  both  these  modes  with  a  CDSPPR  atmospheric  model  and 
the  horizontal  wind  profiles  of  Forbes  and  Gillette  (1982) ,  Lindzen  and 
Hong  (1974),  and  Itable  et  al.  (1977)  .  The  dispersion  of  these  modes  with 
wave  period  at  f imad  propagation  direction  and  with  propagation  direction 
at  fixed  period  have  both  bem  investigated.  We  have  also  investigated  the 
effects  of  viscosity,  non-linearity  and  dissipation  for  all  the  guided 
modes.  In  general,  we  have  found  the  following  results: 

(1)  In  the  absenoe  of  horizontal  winds  the  guided  modes  with  periods 
^xnrter  than  8  min  are  confined  to  the  lower  atanosphere  below  160  )aB. 

(2)  In  the  presence  of  winds  everything  depends  on  the  direction  of 
horlzcntzd  wave  propeigatlon.  If  the  wave  propagates  along  the  direction  of 
the  zonal  %/ind  (due  East) ,  then  the  wave  period  is  Doppler-shifted  towards 
the  longer  periods  and  many  of  the  higher  modes  (small  horizontal  phase 
velocity)  %buld  have  their  wave  energy  spread  over  a  much  greater  height 
range.  If  the  wave  propagates  along  the  direction  of  the  meridional  wind 
(due  North) ,  then  the  North-South  symnetry  of  the  meridicnal  wind  pattern 
has  oonsiderably  less  effect  on  the  guided  modes  «4i06e  wove  energy  remains 
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in  lore  or  less  same  location  as  the  windless  case  (i.e.  belcw  160  kn 
for  periods  less  than  8  min) .  If  the  wave  propagates  agednst  the  directiori 
of  the  zonal  wind  (due  Vtest) ,  then  the  wave  periods  are  shifted  to  the 
shorter  period  end  and  only  one  or  two  lower  modes  (hig^  horizontal  phase 
velocity)  can  exist  with  most  of  their  energy  located  below  80  km.  Dius, 
the  net  effect  of  the  winds  is  to  increase  the  spread  of  the  wave  energy 
over  a  greater  height  range  when  the  propeigaticn  is  adong  the  wind  and 
decrease  it  when  the  prcpagation  is  against  the  wind. 

(3)  The  %rinds  also  have  an  effect  on  the  ffeguenq^  dispersion  of  the 
guided  modes.  If  the  wave  {xrpagates  due  East,  only  the  short-period  modes 
survive;  the  longer  period  ones  have  been  Doppler-shifted  out  of  eadstanoe. 
If  it  is  due  North,  then  the  dispersion  curves  would  zesenble  the  windless 
case  due  agzdn  to  the  synnetry  of  the  meridional  wind.  The  dispersion  with 
respect  to  propagation  direction  at  a  fixed  but  short  period  (6  min.  Fig. 
lO-a)  shows  that  a  nmter  of  modes  vanishes  if  the  propagation  is  due  Hast 
when  most  of  the  higher  modes  are  again  Doppler-shifted  cut  of  existence. 
The  apposite  is  true  for  the  dispersion  curve  at  a  fixed,  but  relatively 
long  period  (14  min.  Fig.  10-b) .  In  this  case  the  higher  modes  tend  to 
survive  only  if  the  propagation  is  due  Hast  when  ttmrm  is  room  for  Dcppla^ 
shifting  to  ehortar  periods. 

(4)  The  oonputations  for  the  effects  of  visooeity,  non-linearl^  and 
instability  seen  to  show  that  the  visooeity  is  unisportant  for  most  heicht 
ranges  of  interest.  In  general,  within  a  height  range  between  110  km  and 
180  km,  the  non-linearity  and  instability  can  be  important  for  many  of  the 
higher  modes.  The  ^jrportanoe  decreases,  with  small  wave  ai^litudes  as 
expected.  For  reasonable  wave  anplitudes  (say  waves  which  produce  15-20% 
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flurtuaticns  in  the  major  species  at  ,100  km)  the  ren-linearity  and 
instability  can  be  totally  neglected  only  for  the  Iamb  mode  with  a 
horizcnt2d  phase  velocity  of  almost  315  m  sec“^.  For  other  modes  they  can 
be  neglected  for  seme  horizontal  phase  velocities  when  the  wave  anplitudes 
happen  to  be  moderate  in  the  height  ranges  mentioned.  So  far  as 
es^ierimental  observations  are  ooncemed,  for  most  of  our  calculations  we 
hcT'e  found  that  the  mesospheric  edrglaw  peaks  at  height  levels  well  below 
the  region  where  viscosity,  veve  instability  and  non-linearity  become 
important  for  gravity  waves,  (we  deal  with  horizontal  phase  velocity  higher 
than  100  Bv^sec  in  this  paper  to  an/oid  the  singulaur  solutions  produced  by 
critical  layers  mentioned  earlier  in  the  paper) . 

Vie  vnuld  lilce  to  suggest  doing  simultaneous  airglow  observations  at 
mesospheric  and  Icwer  thermospheric  height  levels  (e.g.  557.7  no  01,  the  GH 
emission  bands,  etc.)  and  the  ionospheric  height  level  (e.g.  the  630  ran 
01) .  The  data  should  be  taken  at  short  time  intervals  (<  min)  to 
aocounodate  the  shott-period  structure.  This  can  provide  us  with 
information  on  the  energy  distzibution  of  the  short-period  ducted  waves 
vAiich,  as  we  have  shown,  %auld  oonoentrate  »t  lower  edtitudes.  They  can 
also  provide  verification  on  the  long-perioi  and  laujge  scale  TID's  produoed 
by  vmves  whose  energy  distribution  is  oonoehtrated  well  above  100  km 
(Theme,  1968;  Tfti  et  al.,  1980).  So  far  as  jwe  know  Theme's,  (1968)  redzur 
oheervation  eit  Axecibo  wes  the  only  oonfimiation  of  such  an  energy 
distribution  for  l2u:ge  seeds  TTD's  and  that  no  airglow  observations  have 
ever  provided  such  evidence. 
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FIGURE 

1.  Brunt  period  curve  for  the  OOSPAR  atmosphere 

2.  Ooccdinate  system  in  the  horizontal  plane 

3.  Dispersion  curves  in  1±e  OQGEAR  model 

4.  Ihe  windless  modes  at  T  «  6  min 
-a  Wave  Functions 

-b  Kinetic  Eher^  Densities 

5.  Hind  Etofiles 

-a  Zonal  Oonponents 

Meridional  Ocnponents 
Thin  line  for  meeui  circulahicn 
Bold  line  for  mean  plus  tidal  ocnponents 

6.  The  %dnd/  modes  at  T  »  6  min,  and  a  -  Q* 

Have  Functions 

*b  Kinetic  Diergy  Densities 

7.  The  tdndy  modes  at  T  »  6  min,  and  a  «  90* 

-a  .  Have  Functions 
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-b  Kijietic  Eiiergy  Densities 


8.  Dispexsicn  curves  in  the  zonal  idnd 

-a  a  =  0*  (wave  propagation  towards  the  East) 

-ii  a  =  180*  (wave  propagation  towaurds  the  West) 

9.  Dispersion  curves  in  Idle  neridicnal  %dnd 

-a  a  =  ,90*  (wave  propagation  toweuTds  the  North) 

^  o  =  270*  (wave  propagation  towards  the  South) 

10.  Directional  variation  of  Vpi^^  in  the  winds 
-a  for  T  =  6  min 

-b  for  T  »  14  min 

11.  Ihe  three  terns  and  the  Richardson's  nntezs 

Windless  spde,  T  =*  S.Smin,  ii/sec 

-a  Inertial,  viscous,  nonlinear  terns 
a  Inertial  tern 
+  Nonlinear  tern 
*  Viscous  tern 
-4}  Richardson's  numbers 

Bold  line  gives  the  wave  function  in  arbitrary  scale 
Relative  variation  of  density  at  100  )an  is  20% 
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12.  The  three  terns  and  the  Richardson's  numbers 

Hincty  BDde,  T  =  5.5  min,  =  260  Vsec,  o  =  0* 

-a  Inertial,  visocus,  nonlinear  terms 
A  Inertial  term 
+  Nonlinear  tern 

*  Viscous  term 

-b  Richazxiscn's  numbers 

13.  The  three  terns  and  the  Richardson's  rutbers 

Wind/  node,  T  =  12  nin,  =  210  x^soc,  a  =  180 
-a  Inertizd,  viscous,  nonlineaur  terms 
A  Inertial  tern 
-f  Nonlineeu:  term 

*  Viscous  tern 
Richaurxlson's  numbers 
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Height 
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Profiles 


The  windy  inodes  at  T  «  6  min,  and 


Bold  line  gives  the  wave  function  in  arbitrary  scale 
Relative  variation  of  density  at  100  km  is  20% 


The  three  terms  and  the  Richardson's  numbers 
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Figure  13-a.  Inertial,  viscous,  nonlinear  terma  Figure  13-b.  Richardson's  numbers 
A  Inertial  term 


PART  B: 


ON  THE  IMPORTANCE  OF  THE  PURELY  GRAVITATIONALLY  INDUCED 
DENSITY,  PRESSURE  AND  TEMPERATURE  VARIATIONS  IN  GRAVITY 
WAVES;  THEIR  APPLICATION  TO  AIRGLOW  OBSERVATIONS. 


U.  Makhlouf*,  E.  Dewan*,  J.  Isler'  and  T.  F.  Tuan' 


'Physics  Department 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 


'Optical  Physics  Division 
Air  Force  Geophysics  Laboratory 
Hanscom  AFB 

Bedford,  Massachusetts  01731 


55 


^tetcagt 


A  quantitative  study  is  made  on  the  relative  inportanoe  of  the  purely 
gravitationally  induced  ocnpression  (G.I.C. )  due  to  fluid  particle  altitude 
change  and  the  actucil  "wave  oonpiession''  vAiicii  can  occur  at  a  fixed 
altitude  in  a  gravity  wave.  The  results  for  density,  pressure  and 
tegoperature  variations  show  the  following:  (1)  the  G.I.C.  effects 
predaninate  (>95%)  for  v/c  <  20%  where  v  is  the  horizontal  phase  velocity 
euxi  very  sisple  fomulas  can  be  obtained:  (2)  the  relative  i]sport2uxae 
depends  strongly  on  frequency  for  wave  periods  less  than  10  min.,  but 
bftcmifts  totally  independent  of  frequency  for  periods  greater  than  20  min. : 
(3)  the  teaperature  measuranents  can  be  quidcly  converted  to  heig^ 
variations  wherever  the  G.I.C.  effect  predcninates;  in  generza  the 
conversion  is  equivalent  to  the  adiabatic  lapse  rate,  i.e.  a  10* 
tenpexature  variation  oorre^xvids  to  a  height  change  of  1  Im. 

In  ctddition,  the  total  kinetic  energy  density  can  be  sinply  expressed 
in  terms  of  height  vzuriation  and,  vhenever  the  G.I.C.  effects  predaninate, 
can  be  very  easily  obtained  from  tenperature  maasurements.  An  interesting 
by-product  has  been  that  for  waves  of  snail  horizontzd.  phase  speed,  the 
total  wave  )dnetic  energy  at  any  frequency  is  equal  to  the  kinetic  energy 
of  the  natural  (Brunt)  oscillation  of  an  air  parcel  with  the  Mma  vertical 
di^lzx^ement. 


56 


c  Effect 


Ihe  purpose  of  the  present  note  is  to  determine  the  degree  to  vAiich 
Gravity  Waves  (to  be  abbreviated  G.W.)  are  subject  to  purely 
Gravitaticnally  Induced  ejqpansicns  and  ocnpressicns  (G.i.C.)  from  edtitude 
cheuiges  of  a  given  parcel  of  2dr  and  the  degree  to  which  they  are  caused  by 
purely  '*wave  oonpressian".  If,  for  certain  types  of  G.W. ,  the  G.I.C. 
effects  can  cause  most  (say  90%)  of  the  G.W.  variations  so  that  we  can 
neglect  'hove  ccBpression”,  then  the  resulting  siagplification  (Dewan  et  al 
(1988))  can  allow  us  to  xise  very  sinple  fonailas  for  the  G.W.  parameters, 
(density,  pressure  and  temperature  fluctuations)  as  well  as  to  make  quick 
estimates  of  other  different  G.W.  parameters.  For  instance,  a  quick 
oonversion  trm  tenperature  fluctuations  to  height  fluctuations  can  be 
readily  affected  emd  the  latter  cm  be  used  for  a  quick  estimate  of  the 
kinetic  energy  density,  ihis  is  also  very  useful  in  alrglow  maasuramants 

I  .  ,  , 

where  interf ercmeter  measurements  give  tenperature  directly  and  hence 
airglow  height  changes  can  be  estimated.,  Vhils  soma  of  our  results  are 
kncwn  queditatively  to  G.  W.  specjalists,  to  our  ]aiawledge  there  has  never 
been  a  systematic  quantitative  tr^tment  of  the  subject. 


We  will  begin  by  shewing  that  the  density,  pressure  and  temperature 


fluctuations  can  be  explicitly 


as  a  sin  of  two  terms;  one 


describing  the  G.I.C.  oonpression,  while  the  other  the  'Snve  motion".  Ihe 
former  vanishes  in  the  absenoe  of  gravity  (g*6)  while  the  latter  remains  in 
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atltered  fora  and  prcn^ides  for  purely  longitudinal  aocustic  '*wave 
ocnpressicn". 

For  q^Q,  the  relative  infxsrtanoe  of  these  two  effects  depends  on.  the 

negnitude  of  the  horizontal  phase  velocity  of  the  G.Vf.  as  ocnpared  to  the 

speed  of  sound.  For  ssall  scale  G.W. ,  uhere  by  definition  ue  neem  the 

phase  velocity  is  such  less  tham  the  speed  of  sound  and  the  wavelength  is 

also  snail,  we  have  found  that  indeed  the  G.I.C.  effects  predcminate. 

Since  much  of  the  oosetved  ndd-latitude  airglow  fluctuations  (incltriing  CH) 

involve  relatively  snail-scale  G.W.  (Witt  (1962) ,  aairaidl  et  al  (1985) , 

Taylor  &  Hopgood  (1968)  with  low  horlzbntal  phase  velocity  v  <  70  m  sec~^, 

we  say,  for  most  cases,  assume  that  the  density,  pressure  and  tsaparature 

variations  axe  G.I.C.  and  a  qaldc  oonversion  to  height  variations  nay  be 

made.  Very  sisple  fomulas  for  the  kinetic  energy  density  can  be  derived, 

2 

where  >0,  both  for  the  genered  case  and  even  note  so  in  the  G.I.C. 
approxisation. 

We  should  point  out  that  there  is  considerable  inberferanca  between 
the  two  effects  and  that  only  if  either  is  greater  in  aegnltuds  than  the 
other  by  a  large  factor  can  we  oanfortably  neglect  the  ssallar  terms.  To 
test  the  limits  of  this  criterion  we  will  consider  three  specific 
experimental  obsexvatlons;  the  first  being  the  airglow  dbeexvation  made  by 
Taylor  and  Hapgood  (1988)  (in  association  with  the  MKFSIUt  program)  which 
easily  meets  our  criterion;  the  second  being  the  wave  ofasexvatione  by 
dairemidi  et  al  (1985)  which  also  meets  cur  criterion;  third  beinr  j 
observations  on  noctilucent  clouds  made  by  Witt  (1962)  and  is  mrginally 
edxive  the  limits  of  our  criterion.  The  results  show  that  for  all  three 


58 


sets  of  cteervations,  the  sinple  C.I.C.  nodel  is  surprisingly  accurate  <uid 
that  cne  nay  safely  assume  that  such  of  mid-altitude  density,  pressure  and 
teoperature  fluctuations  is  governed  by  G.I.C.s.  In  general,  ve  have 
proved  that  in  the  region  v«C  where  the  G.I.C.  esqpansion  or  contraction 
effects  dominate,  a  10* K  vzuriation  in  tenperature  would  indicate  an 
edtitude  change  of  1  km  of  a  parcel  of  edr.  He  have  edso  found  that  in  the 
sane  region  the  kinetic  energy  density  is  equivzdent  to  tiiat  of  the  free 
vertical  oscillation  of  an  air  parcel  with  the  same  verticad  displaced 
anplitude  and  is  independent  of  the  wave  frequency.  Ihus,  ttw  above 
conversion  czm  be  ianediately  used  to  calculate  the  kinetic  energy  density. 


Hhen  a  fluid  element  is  displaced  by  a  vertical  distance,  h,  in  a 
hydrostatic  flviid,  its  density  is  ocnpressed  by  an  amdunt  This 

daisity  change,  Apa»  is  produced  by  a  change  in  background  pressure,  ap, 
(vAiere  ap  ••  (dp^^/dz)  h  «  and  Aps  is  the  density  variation  due  purely 

to  adiabatic  ocnpressicn)  minus  the  change  in  background  density,  a^q  • 
(dpf/dz)h,  one  can  show  that  in  an  uniform  isothermal  bacJo^xund 
atmo^here  with  scale  height  H, 

^£a  a  .  ^£a  (1) 

Po  Po  Po 
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liiere  we  have  used  the  f2kct  that  c^PoTrPo 


and 


i£o  .  . 1.  iBd 


Pq 


Po  32 


III.  A  ocnoarison  of  the  G.I.C.  and  wave  ocnpressioo 


In  a  G.W.  with  frecpency  «,  h  ■  w/i»  (to  first  order)  vtere  h  and  w 
aure  the  vertipal  di^laoement  and,  velocity  fields  respectively.  At  the 
Brunt  frequency,  m  «■  the  atmosphere  %iould  just  oscillate  vertically  at 
its  natural  frequency  %d.th  the  density  tluctuations  qiven  ty  (1)  •  In 
general,  by  using  the  Hines*  (1960)  fonnila  for  Ap/Po»  h  •* 
horizontal  velocity  field  u  wo  can  show  that  to  the  first  order  (see 
Appendix  (I)): 


l)h  ^ 

•  H 

(-)^  (-) 

'c*  W 

Apa 

+  ^ 

(2) 

Po 

Po 

where  u  ■  horizontal  particle  velocity 

V  -  Vx  ■*  horizontal  phase  velocity 

Apc/po  ■  (v/c)2  (Vv)  -  density  fluctuation  due  to  %nve  oo^pcession 
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Fran  (2)  it  is  cleaur  that  -*  0  as  g  -*  0  sinoe  H  *  C?/yq,  and  the 

terms  tpc/Po  will  be  c^propriately  changed  to  provide  for  the  purely 
longitudincil  "wave  oonpression"  at  Jd.1  frequencies  for  sound  wave 
propeigation  in  the  absence  of  gravity.  As  g  -*  0  ,  H  ■*  «  and  using  (A22)  we 
obtain 

lA£i  - )  LAf»d.  ,(2a) 

emd  iuj_  V  kv 
|w|  ^-*o  kg 

where  is  the  horizontal  wave  vector.  (2a}.  clearly  shows  a  purely 
longitudinal  wave  ocnpression  when  g  »  0.  Frtaa  now  on  we  shall  only  be 
concerned  with  gp^O.  In  that  case,  from  (2)  the  relative  iiqportance  of 
^Pa/Po  and  Apc/po  depends  on  (v/c)  and,  in  order  to  remain  well  below  the 
Lindzen  (1981)  liiait  for  G.W.  saturation,  iu|  <  v  (see  Ftitts,  1984) . 
deeurly  for  anall-scade  G.W.  where  v«C,  the  densi^  veuriation  ccnes  nainly 
frcn  G.I.C.  conpr^ion,  tpa/Po-  ^  the  other  h2md  for  v  -»  C  and  close  to 
the  limit  for  G.  W.  saturation  the  wave  ocnpression  Upc/po  ^  large  and  we 
need  to  include  both  effects. 

In  general,  since  both  Ap^pQ  and  tp^p^  are  cxmplex  quantities  with 
different  phases,  considerable  interference  occurs  between  the  two  terms 
and  only  if  one  has  a  magnitude  of  at  least  4  to  5  times  greater  than  the 
other  can  we  neglect  the  other  term. 

Tb  make  a  quantitative  conparison  between  Ap^  and  Ap^,  we  plot 
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jAPal/|A^|  and  IapcI/I^pI  ^  a  function  of  v/c.  It  can  be  shown  for  ]^  >  0, 
(see  Appendix  U) , 


lA£ai 
|Ap  I 


^(7-1)  -  (u./u,b)^(v/c)2  J 


lAp  I 


=  [ 


(7-I)  -  {«/wb)^(v/c) 


's' 


Ihus,  in  general  |apal/|Ap|  begins  %d.th  the  value  unity  (i.e.  densi'^ 
variations  are  due  entirely  to  G.I.C.  ocnpression)  at  v/c  ■  0,  while. 
I^PcI/i^^l  t^egins  at  zero  (i.e.  no  "wave  ocDpression") .  As  can  be  seen 
frcm  (3)  and  (4),  when  v/c  increases,  the  behavior  of 
IVci/l^^l  (impend  entirely  on  w/w^  for  any  value  of  v/c. 

For  long  period  gravity  waves  where  w  «  we  any  neglect  the  term 
involving  w/w],  in  the  denominator  for  both  |apa[/|Ap|  and  |apcl/l^^i* 
such  cases 


and 


lA£ai 
|Ap  I 


2 


|Apf;|  , 

|Ap  I 


(5) 

(6) 


and  |Apal/|Ap|  drops  from  unity  as  v/c  Increases  while  {ApcI/I^aI  increaaes 
lineeurly  with  v/c. 

Fig.  1  shows  this  behavior  for  gravity  waves  with  periods  greater  than 
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20  min.  Actually,  there  is  little  vauriation  with  period  fran  20  min.  c>rt 
and  a  plot  at  r  »  120  min.  falls  exactly  on  top  of  the  oorresponding  oirves 
at  20  minutes.  Both  curves  are  plotted  frcm  (3)  amd  (4),  but  there  would 
have  been  no  difference  if  we  had  used  (5)  and  (6) .  At  r  =  1C  min. ,  (Fig. 
2),  |Apcl/|Ap|  begins  to  deviate  frcn  a  straight  line  but  only  slightly. 
|Apa|/|Ap|  drops  scncuhat  less  as  v/c  increases  showing  aui  irtcrease  in 
infxsrtance  relative  to  '\jave  ocnfsession"  for  a  given  phase  velocity. 
Furthermore,  the  cross  ever  betwenu  the  G.I.C.  oenpression  and  the  "wave 
oonpression"  mains  fairly  constant  at  between  v/c  *  55  to  60%.  In 
general,  tor  G.W.  with  freqi>ft.-yjies  in  this  range  (i.e.  r  >  10  min.),  if  we 
use  the  approxiaate  criteria  tlat  |Apal/|Apcl  ^  ^  greater  than  a 
factor  of  5,  then  v/c  t)aa  to  be  less  than  15%  or  the  horizontal  phase 
velocity  has  to  be  less  than  45  m  sec**^.  Diis  covers  a  fairly  large  range 
and  include  most  of  the  observed  midhaltitude  small  scale  G.W.  fxon  aiiglow 
data.  ■ 

Ihe  verticaO.  green  line  ii^cates  the  position  where  •  0.  It  is 
easy  to  show  frcn  the  Hines'  di^»rsicn  fonaila  that  )c§  -  0  corresponds  to 


to  the  left  of  the  grean  .  kj  >  0  and  G.W.  may  propagate  freely.  To  its 
right  kj  <  0  and  the  G.W.  is  evanescent  .in  the  vertical  direction,  although 
it  still  propagates  along  k^.  Ihe  dotted  lines  represent  the  evanescent 


region.  At  v/c  >  1  we  have  the  Laab  node  for  liiich  there  is  no  G.I.C. 
ocDpression  since  there  is  no  verticzd  notion  euid  variations  in  density  are 
caused  entirely  by  "Vave  oon^uression",  an  eKZKrt  reversal  of  vAiat  twippens  at 
v/c  >  0.  The  scune  reverscil  occurs  in  all  the  Figures. 

As  r  -  r^,  the  Brunt  period,  both  |A^cl/l^^l  incr®*®® 

with  v/c,  (see  Fig.  3  &  4)  in  such  a  way  that  the  G.I.C.  ocnpzessicn 
continues  to  predoninate  over  a  ULXh  greater  range  of  horizontal  phase 
velocity.  In  feet,  in  Fig.  *4  where  r  «  5.2  odn.,  ti)e  G.I.C.  oespressien 
oentinues  to  deninate  all  the  way  to  the  green  line  vhere  v/c  *  42% 
corresponding  to  a  horizontal  phase  velocity  of  well  ewer  120  n  aec.'^;  a 
strictly  neditarsczLle  G.W.  (Fiancis  1975) .  Thus,  as  the  G.H.  period 
approaches  the  Brunt  period,  the  G.I.C.  cenprassion  doodnatas  cK/er  a  larger 
range  of  v  and  hence  a  larger  ranc^  of  wavelengths. 


The  pressure  veuriation  can  be  easily  shown  to  be  the  sane  as  ^pq. 
Using  the  linearized  horizontal  ncmentun  conservation  equation,  it  can  be 
^heur.  that  for  xnonochrcBatic  waves. 


or 


so 


AP  -  VpfJH 

4E  .  ivu 

Po 


MeJ—  .  l^ri 

C^|A/»|  |A^| 


(8) 


Frcn  Fig.  l->4  we  see  that  there  is  little  pressure  variation  produoad 
Ijy  "wave  oesprassion"  for  snail-scale  G.W.  Only  as  the  scale-size  of  the 


G.W.  incxeeises  does  it  becxme  inportant.  Ihus,  for  snail-scale  G.W.  (<45  ra 
sec“^) ,  the  pressure  vciriations  are  relatively  insignificant. 


The  tenperature  variations  may  be  very  eeisily  obtained  from  the 
linearized  perfect  gas  law  given  by 


To  Po  Po 


Frcta  (1) ,  (9)  ar^  (2)  /  we  obtain. 


To 


(7  -  1)  (P) 


4£a 

Po 


(9) 


(10) 


11x03,  for  V  «  C  (snail-scale  G.W.), 

^  *  -  4£a  =  -  ,  (7  -  1)  h 

To  Po  7  H 


(11) 


and  the  tenp.  variation  is  apposite  in  phase  to  the  G.I.C.  density 
ocnpression.  We  will  define 

^T  _  alia  ATb 


(12) 


where 


AT. 


-  4£a  .  ^Ib  ^ 

Pa  To  ^c2 


^o  Po  *o 

Ag2dn,  as  in  (3)  &  (4)  it  is  possible  to  show  that  for  >  0,  (see 
^pendix  III), 


|AT| 


[ — 7-1 V2 

'■  1  »  (7-l)(w/«-b)^(V/C)^  -I 


(13) 


-  r-  -  (Trim  2)  1  V 

I  ATI  ^  1  -  (7-I)  (w/whPiv/C)^  ■» 


(7-I)  (w/wb)^(V/C)- 


(14) 


Ihus,  once  again,  for  v/c  »  0  |ATal/|ATl  *»  1  and  |aTcI/|AT|  »  0  and  we  see 
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that  in  the  lew  pheise  velocity  limit  (very  small-scede  G.W.)  the 
tenperatuxe  veuriaticn  is  caused  entirely  by  the  adiabatic  o^ansion  and 
ccppression  of  an  adr  paroel  oscillating  in  the  vertical  direction.  For 
long  period  G.W.  (<.>«u)b) ,  (13)  &  (K)  siitplify  to; 

-  t  1  -  (|)2  ]V2  (15) 

.  (7-1)V2(3)  (16) 

Conparing  (15)  and  (16)  with  (5)  and  (6)  we  see  that  for  long  period  G.W. 
|aTal/|AT|  behaves  the  same  way  as  |A^a/I^A>l  in  (5)  whilst  |aT^|/|aT|  hais  a 
far  more  gradual  slope  than  |Apcl/l^/’l>  "Oius,  for  tespexature  variations, 
the  G.I.e.  effect  predoninates  over  the  "wave  effect"  through  2^  greater 
range  of  phase  velocities. 

All  this  can  be  seen  from  Fig.  5  which  shows  Ihe  behavior  of 
|ATal/|AT|  and  |aT(.|/|aT|  for  G.W.  with  periods  fren  30  min.  on  up.  Again, 
a  plot  of  both  curv^  at  r  »  120  min.  can  be  exactly  superiaposed  over  the 
two  curves  for  r  ■  30  min.  Ihis  is  similar  to  tite  curve  of  the  density 
Veuriations  (Fig.  1)  which  has  a  lower  limit  at  r  >■  20  min.  Again,  for  r  > 
30  min.  there  is  no  significant  difference  between  (13)  and  (15)  or  (14) 
and  (16) . 

Fig.  6,  7  and  8  show  the  behavior  of  tenperature  variations  for 
f  =  10,  6.4  eud  5.2  min.  respectively.  The  big  difference  between  these 
tenperature  veuriations  and  the  density  veuriations  given  by  Fig.  2,  3  and  4 


is  that  for  the  former  the  effect  of  G.I.C.  ocnpression  is  even  much  more 
inportant  than  the  "wave  ocnpression".  For  instance,  at  r  =6.4  min.,  the 
wave  ocnpression  effect  overtakes  the  G.I.C.  effect  at  v/c  =  0.72,  whereas 
the  G.I.C.  effect  remains  predominant  in  the  corresponding  tenperature 
veuriation  all  the  way  to  the  green  line.  For  both,  however,  the  G.I.C. 
effect  predominates  cill  the  way  to  rhe  green  line  when  r  =  5.2  min., 
approaching  the  Brvint  period. 

Ihe  dotted  curves  in  adl  the  Figures  on  the  right-hard  side  of  the 
green  line  oorre^xxd  to  evanescent  waves  with  kl  <  0  and  thus  a  purely 
imaginary  In  thi^  region  one  can  no  longer  use  the  relatively  sinple 
ejgnessions  given  by  (3) ,  (4) ,  (13)  and  (14) ,  which  were  derived  for  ]^  >  0 
and  real  Insteztd  one  has  to  work  with  the  general  polarization 
relations  given  by  Hines  (1960) . 

From  the  above  discussions,  it  is  cleeu:  that  ocnparatively  snail-scale  G.H. 
with  horizontad  phcise  velocities  less  than  40-50  m  sec'^  may  be  considered 
to  be  purely  under  the  influence  of  G.I.C.  expansion  and  ocnpression  (> 

901) .  Ihis  will  allow  us  to  make  sinple  conversions  from  the  measurement  . 
of  tenperature  fluctuations  to  hei^t  fluctuations.  Ihus, 

im  .  iaiai  .  llZil  M  (17) 

To  To  7  H  '  ^ 

eud  using  the  usual  values  for  7,  H  and  To.  we  obtain  the  appraximate 
oortversion  that  a  10*  variation  in  tenperature  oorre^jcnds  approocimately  to 
a  height  vzuriation  of  i  km.  at  90-100  km  altitude. 


Dlls  GUI  be  isaiediately  used  for  cxxixiting  the  time  averaged  kinetuLc 
energy  density,  K.E. ,  given  by 

K.E.  =  1/4  po  (Iu|^  +  |w|2)  (18) 

vhere  w  =  vertical  particle  velocity.  Using  the  continuity  equation  and 
the  Hines'  di^iersicn  relations  one  can  show  that  (Appendix  IV) , 

K.E.  «  i  p<^g|h|2  t^-TJy/^)_^v/c)2  ^ 

For  large-scale  G.W.  .vhere  v  -»  c,  the  dencninator  can  be  much  less 
than  the  numerator.  Diis  is  e^iecially  true  for  long  period  G.W.  Du:^, 
the  K.E.  can  then  be  much  greater  than  1/4  />d-i|h|2.  ^  small-scale  G.W. 
(v«c) ,  on  the  other  hand,  the  K.E.  is  very  accurately  given  by: 

K.E.  =  1/4  Pc^g|h|2  (20) 

It  is  interesting  to  note  that  (20)  is  true  irrespective  of  the 
frequency  of  the  small-scale  G.W.  Since  (20)  is  also  the  K.E.  for  the 
natural  Brunt  oscillation  of  a  parcel  of  air,  we  may  ocnclude  that  the 
total  K.E.  of  small-scade  G.W.  at  any  frequency  is  equivalent  to  natural 
vertical  atmospheric  oscillations  with  the  same  vertical  displacement. 

HGre,  the  G.I.C.  approximatian  for  the  density,  pressure  and  taiperature 
variations  ate  also  good.  Using  (17),  (20)  acquires  the  form 
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(20a) 


K.E. 

Po 


i  (JL_ji4Ta,2 
4  S-r  't„  ' 


i  /_1_\i4Ti2 

4  S-l^'Tn’ 


0.8751—1 


2 


vAiere  pto  =  badoground  pressure  energy.  Thus,  a  mecisuresnent  of  the 
rotational  temperature  Vciriation  can  unnediately  yield  the  total  kinetic 
energy  density  at  any  frequency  for  sjnall-scale  G.W.  In  fact,  (19)  shows 
that  the  K.E.  depends  strongly  on  frequency  only  for  large-scade  G.W.  where 
V  -*  c. 

To  sutmarize:  (1)  for  small-scale  G.W.  (v  <  45-50  n  sec"^) ,  the  K.E. 
doisity  is  the  same  as  that  for  natural  atmospheric  oscillations  and  is 
ind^aendent  of  wave  frequency;  (2)  for  ary  G.W.  (eudpitzary  v/c) ,  the  K.E. 
eigfciin  beccmes  equal  to  natural  atmo^sh^nx:  oscillations  when  w  -*  (.>b  ^ 
expected:  (3)  for  large-scede  G.W.  (v/c  .*  1),  the  K.E.  for  the  sane 
verticad  di^laoeinent  become  very  large  euid,  unli3oe  the  two  previous  cases, 
can  also  become  very  strongly  frecpiency  dependent. 


(XV)  ArpJjgatiaig 

To  ^Jecifically  apply  the  edxve  anedyses,  we  apply  than  to  three 
^»cific  sets  of  observed  data:  two  sets  involve  G.W.  peurameters  well 
within  the  above  mentioned  limits  for  the  validity  of  the  purely  G.l.C. 
model;  the  other  with  parameters  above  the  borderline.  The  first  set  were 
data  taken  by  Taylor  ctnd  Kapgood  (1988).  The  observed  data  include  a  G.W. 
with  a  horizant2d  phase  velocity  of  v  >■  21.7  n  sec'*^,  a  horizontal  wave 
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length  of  ^  ^  period  r  >  20  min.  Using  g  *■  9.8  n  sec**^  H  >■ 

6.56  ]an,  7  »  1.4,  rj^  «  5.1  min.  (the  Brunt  period)  and  the  ^>eed  of  sound  c 
=  300  m  /sec”^  ,  we  find  that  if  we  enploy  the  full  expression  given  by 
(10)  with  help  fzxxD  the  Hines'  polarization  and  dispersion  relations  (Hines 
1960).  we  obtain 

=  4.364  X  10"2  (h|  (21) 

io 

where  h  is  expressed  in  )an.  Using  the  C.I.C.  2ipprtadmation,  we  obtain 

*  ,(2^  M  «  4.355  X  10“2  |h|  (22) 

ihe  error  from  using  (22)  is  only  0.2». 

secxand  set  were  taken  by  Clairenidi,  Hersh  and  Moceels  (CLairavidi  et 
ail  (1985) ) .  They  have  measured  v^ives  with  velocity  v  <■  15.6  m  sec”^  and  a 
period  of  48  min.  The  full  tenperature  fluctuation  without  aipproximation 
is  given  by 

^  -  4.36  X  10-2|h|  (23) 

*0 

Ccrpared  with  the  G.I.C.  approximation  given  by  (22),  we  have  an  error  of 
the  order  0.1%. 

Ihe  third  set  of  data  were  taken  by  Witt  (1962) .  The  G.W.  had  v  ■  75 
m  sec~^,  Ax  *  ^  of  the  parameters  eu»  the  >«*■»«»  as  the 

first  data  set.  We  obtain 
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(24) 


»  4.5  X  10"2  |h| 

Ttie  G.I.C.  £q:proxuatian  gives  the  same  result  cis  (22).  Ihus,  the  error 
for  this  exanple  is  3%.  The  saraa4)at  bigger  error  is  understandable 
because  at  75  m  sec~^  the  horizontal  phase  velocity  is  scmsMhat  above  the 
45-50  a  sec**^  upper  limit  that  we  mentioned  eeurlier  as  the  safe  limit. 

(i.e.  v/c  *  15-20t) .  Actually  for  this  exanple  there  is  considerable 
interference  between  the  two  terms  on  the  R.H.S.  of  (10) . 

(V)  ft)rglMsi<?n 

The  present  note  seems  to  show  that  for  small-scale  G.W.  (say  v/c  < 
18%),  the  G.I.C.  approximation  is  very  good  to  above  95%  in  tenperature  and 
density  flxictuations.  The  appraxisation  appears  to  be  applicable  to  nuch 
of  the  observed  mid-altitude  airglow  data.  Quick  estimates  of  euid  simple 
formulas  for  G.W.  parameters  such  as  height  veuriations  and  kinetic  energy 
density  may  be  immediately  obtained  from  meeisurement  of  tenperatuze 
V2uriaticns. 
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Figure  captions 


Fig.  1  a  plot  of  the  Induced  density  variation  due  to  the  gravitational 
(red)  and  the  euxustic  ocnpressions  (bl2Kdc}  as  a  function  of  v/c'  for  wave 
periods  >  20  min. 

Fig.  2,3,4  plots  of  the  induoed  density  variation  due  to  the  grztvitationzd 
(red)  cuid  the  acoustic  oonptessions  (black)  as  a  function  of  v/c  for  wave 
periods  »10  ,  6.4  ,  5.2  min.  respectively. 

Fig.  5  a  plot  of  the  induoed  tenperature  variation  due  to  the 
gracvitatlonal  (red)  and  the  acoustic  ocnpressions  (black)  as  a  function  of 
v/c  for  wave  periods  >  30  min. 

Fig.  6,7,8  plots  of  the  induoed  tenperature  variation  due  to  the 
gravitationed.  (red)  and  the  acoustic  ocnpressions  (black)  as  a  function  of 
v/c  for  wave  periods  »  10,  6.4,  5.2  min.  respectively. 

Fig.  9  a  plot  of  the  ratio  of  the  kinetic  energy  density  to  the  backgrcud 
pressure  as  a  function  of  v/c  for  wave  periods  of.  20,  10,  6.4  &  5.2  min. 
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Fiquro  I  -  a  plot  of,  the  imluccd  variation  due 
to  tho  gravitational  (red)  and  the  acoustic 
compressions  (blac)()  as  a  function  of  v/c  for 
wave  periods  20  min. 
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Figures  2,  3,  4  -  plots  of  the  induced  density 
variation  due  to  the  gravitational  (red)  and 
the  acoustic  compressions  (blacic)  as  a  function 
of  v/c  for  wave  |>eriods  •  10,  .6.4,  S.2  min. 
respectively. 


-JOrn 


CHtiM 


uccc 


HK  79JiiA 


Fiijuris  3,  4  -  plots  of  the  induced  density 
variation  duo  to  the  gravitational  (rod)  and  the 


acoustic  O'.)!*!!  rt  jsions  (black)  as  a  function  of  v/c 
for  wave  peri<^ds  “  10,  f>.4,  5.2  min.  respective  It. 
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Fi<jures  2,  3,  4  -  plots  of  tho  induced  density  variation 
due  to  the  gravitational  (red)  .ind  the  acoustic  compressions 
(black)  as  a  function  of  v/c  for  wave  periods  *  10, 

6.4,  5.2  min  resf)octively . 
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i'lrjuii  S  -  .1  (>l()t  <jf  the  induced  tomperature 
variation  duo  to  the  qravitational  (red)  and 
th(;  acoustic  compressions  (black)  as  a  function 
of  v/c  for  wqvo  periods  30  min. 
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Figures  6,  7,  8  -  plots  of  the  induced  temperature 
variation  due  to  the  gravitational  (red)  and  the 
acoustic  compressions  (blacit)  as  a  function  of 
v/c  for  wave  [x^riods  *.  10,  6.4,  5.2  min.  respecti 
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7,  8  -  plots  of  the  induced  temperature 
due  to  the  gravitational  (red)  and  the 
acoustic  dompressions  (black)  as  a  function  of 
v/c  for  wajve  periods  ■  10,  6.4,  5.2  min.  respectively 
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Fi<^ir.!  J  -  a  I'lot  of  the  ratio  of  the  kinetic 
eii<  rqy  dt.MsUy  to  the  backqround  pressure  as  a 
f.inctn^ti  of  v/c  for  wave  jHiriods  of  20,  10, 

0,4  and  0,2  tniti. 


Acpendix  I 


We  can  derive  equation  (2)  by  using  the  usual  expression  for  tp/po, 

(Hines  (1960) ) ,  given  by 

^  =  {  w2k2  +  i[(7  -  l)gl<x  - 

where  i(wt-k.r) 

We  elininate  w  by  using  w.s  kxV  2uxl  substitute  in  for  the. scale. hei^t 
H  “  c^/7g  to  obtain, 

^  -  {k^kzv2  +  ik^[(2^)i  -  2^])f  (A2) 

we  obtain  by  adding  2md  subtracting  (^)^  ^ 


^  -  (4i=zv2  +  u^td  -  (|)2)(^)  -  "  1 

)J)f 

(A3) 

Introducing  the  vertical  disylaoanent  field, 

a  -  35i  -  id's'll  - 

• 

(M) 

WB  get  ^  -  (I^)g  +  y^»z  -  id  - 

(A5) 

Finally,  wa  use  u  ■  )«^c2v  [k^  -  i(l  -  7/2)g/c2]f  ^  obtain 

(A6) 

which  is  exactly  equation  (2). 

ftnffPTrily  TT 

■ 

Fron  now  on  we  consider  only  kj>0  (i.e.  kz  real) . 
(3) ,  we  note  (using  (A4) ) , 

To  obtain  equation 

-  (7-l))c|g(l  -  (|)2]  |f| 

-  1)  w2  |f| 

(A7) 
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so 


(A14) 


JAfici  =  ( _ _ _ ,)i/2 

|ap1  ‘ki  +  [((7-1)/7)(1/H)V  + 


k9  + 


where  we  have  used  (AS)  amd  (A13) . 
Using  (A9)  we  obtain: 


so 


^  ‘(^  -^2,  *  (1^1 
.  /■  ,  _ _ — - )V2 

i^Sci  .  f - - — ^ — _)V2(v/c5  (A15) 

-lap  I  ^  (7-1)  -  («/«-b)^  (v/c)^  ^ 


which  is  equation  (4) . 


Frtn  (10)  we  write 


42 » .  tzm 

To  ”  tH 


-  (7-l)(v/c2[„,t^(j;^  +  iij))  -  ixig(v/c)2(c2/v2  -  l))f  (A16) 
vhere  we  have  used  (A4)  &  (A13) .  Thus, 


AT/To  -  (7-l))^{)Cz  +  Ht  “  g/c2(c2/v2  -  l)J)f 


since 


4la  «  .  4£a 


lATal 

|AT| 


(see  (12)  ),  we  obtain  by  using  (A7) 


(7-l)()^v2g/c2), 


(7-l))cgv2  ()^  +  (,  -  (q^c2)o2jV2 


(A17) 


Using  agaih  the  dispersion  relation  for  )c^,  we  obtain 
lAT^l  -  ■  ^  ^  (q/c2)« 

I^Ti  ji^!7TlprTpr2jV2 


.  — jV2 

*  i  -  («^)/g^  * 

so  1^  =  ,_ _ 1  -  (v/c)^  ■  ■■■  --.V2 

|AT|  ‘  1  -  (7-l)(u/ub)^(V/C)^  ’ 

whic*  is  equation  (13) . 


(A18) 


Equation  (14)  say  be  similarly  obtained.  The  denoainator  is  the  sans 

as  (13)  as  expected,  but  the  runerator  is  given  by  t(7-l)uv/c^t  vAiicit  is 
identiczd  to  (4)  except  for  the  c^dditioned  (7-I)  factor. 


jy 

Ihe  K.E.  density  is  given  by: 

K.E.  -  V4  Po  (|ul2  +  lw|2)  (A19) 


Trap  Hines  (1960)  ^^ld  the  definition  of  7,  we  nay  write 
u  -  w)Cj{c2  [kz  +  i»?]f 
w  ■  -  (1  -  (v/c)2jf 


so 


|U| 


2  . 


ki  +  f»- 


kid  -  (v/c)2)2 


|W|‘ 


Ftora  the  disperiscn  relation 

.2  . 


£md 


-  ,X,2 


Hmoe,  |u|^  .  |w|^  -  |w|»(l  . 


(A20) 

(A21) 


(A22) 


Ihe  total  K.E.  then  booones: 


Obviously  (A23)  is  eqiiation  (19). 


89 


III. 


NONLINEAR  RESPONSES 


PART  A: 


On  the  Nonlinear  Response  of  Airgiow 
to  Linear  Gravity  Waves 


.  R.  Isler’*,  T.  F.  Tuan',  R.  H.  Picard'  and  U.  Makhlouf'* 


*  Physics  Department 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 


#  Optical  Environment  Division  (OPS) 
Geophysics  Directorate,  Phillips  Laboratory 
Hanscom  AFB,  Bedford,  MA  01731 


+  Present  Address:  Geophysical  Institute 
University  of  Alaska,  Fairbanks,  AK  99775 


S  Present  Address:  Optical  Environment  Division  (OPS) 
Geophysics  Directorate,  Phillips  Laboratory 
Hanscom  AFB,  Bedford,  MA  01731 


90 


Abstract 


The  question  of  whether  a  linear  gravity  wave  will  give  rise  to 
nonlinear  effects  in  ground-based  airglow  observations  is  important 
for  the  proper  interpretation  of  gravity  wave  dynamics.  In  this 
paper,  we  obtain  a  closed  form  solution  for  the  integrated  airglow 
response  to  a  linear  gravity  wave,  containing  all  the  higher-order 
nonlinear  response  terms.  A  comparison  is  made  to  the  linear 
response,  and  the  higher  orders  are  seen  to  be  significant.  In 
addition,  the  wave-induced  airglow  intensity  fluctuations  are  shown 
to  be  much  greater  than  the  corresponding  major  species  density 
fluctuations. 
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1.  Introduction 


Ground-based  airglow  observations  have  often  been  used  for 
studying  middle  and  upper  atmospheric  dynamics.  The  common 
assumption  has  been  that  the  airglow  intensity  (or  "brightness") 
fluctuations  correspond  directly  to  atmospheric  disturbances. 
However,  if  the  airglow  response  is  nonlinear,  many  of  the  airglow 
power  spectrum  peaks  may  be  due  to  the  nonlinearity  of  the  airglow 
response,  rather  than  an  indication  of  the  actual  dynamic  behavior  of 
the  atmosphere  as  a  whole.  This  paper  deals  primarily  with  the 
purely  dynamical  effects  of  a  linear  gravity  wave  on  airglow  as 
observed  by  ground-based  equipment.  We  ignore,  as  have  other 
authors,  photochemistry,  quenching,  and  molecular/eddy  diffusion 
on  the  time  scale  of  the  wave  period. 

Minor  species  involved  in  airglow  reactions  often  have  a  layered 
structure,  i.e.  a  sharply  peaked  unperturbed  density  profile,  with  a 
very  steep  vertical  gradient  just  below  the  peak.  As  a  consequence, 
it  has  been  realized  for  some  time  that  the  fluctuation  of  the  minor 
constituent's  density,  induced  by  a  linear  gravity  wave,  can  be  very 
large  at  a  point  where  its  unperturbed  gradient  is  steep.  For 
example,  a  linear  gravity  wave  which  produces  a  10%  density 
fluctuation  in  the  major  species  may  induce  a  50-100%  density 
fluctuation  in  the  minor  species  at  the  same  point.  Thus,  the  local 
response  of  the  minor  species  at  particular  height  levels  may  not  be 
amenable  to  a  linear  treatment. 

Weinstock  (1978)  computed  the  first-order,  linear  airglow 
response  and  the  unperturbed  background  airglow  for  02( *£■*■£) 
emissions.  He  also  attempted  to  calculate  the  nonlinear  response. 
Hines  and  Tarasick  (1987)  (hereafter,  HT87)  have  recently 
questioned  the  necessity  of  invoking  the  higher-order  nonlinear 
response  terms  produced  by  the  steep  gradient.  They  have  reasoned 
that,  while  the  local  response  rnay  be  nonlinear,  the  integrated 
response  as  measured  from  the  ground  should  be  linear,  since  the 
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effect  of  the  large  density  gradient  can  be  removed  by  the 
integration.  They  have  demonstrated  this  by  means  oT  a 
transformation  from  the  Eulerian  to  the  semi-Lagrangian  '"rame,  in 
which  the  air  parcels  in  motion  due  to  the  gravity  wave  are  mapped 
b.ack  to  their  original  resting  positions  (before  the  onset  of  the 
gravity  wave).  This  is  a  very  important  result,  since  it  permits,  if 
true,  a  great  simplification  in  calculating  tiie  column-integrated 
response  of  airglow  to  a  gravity  wave.  Radar  or  resonance  lidar 
which  measures  the  local  response  of  a  minor  specie:.-  must  contend 
with  nonlinearities  associated  with  steep  graJ'  ts,  leading  to 
harmonics  not  present  in  the  original  wa\e  pe'turbaiioii;  on  the 
other  hand,  passive  instruments  such  rs  ph-‘Jo.neters  and 
radiometers,  which  measure  total  column  brightness,  effectively 
filter  out  such  nonlinear  responses. 

In  this  paper  we  reexamine  the  order  nonlinear  terms  in 

the  semi-Lagrangian  formulation  of  HZZl  (for  convenience  "semi" 
will  be  dropped  for  the  remainder  of  the  paper).  We  propose  to 
show  that  whilst  liic  steep  ve.tirs:  tensity  gradient  of  the  minor 
species  no  tOuy,ex  shows  up  cv;  Lagrangian  system,  the  higher- 
order  response  terms  *ma;i  important.  This  is  because  a  new 
nonlinearity  is  introduerd  k/  th;,  Lagmagian  transformation  itself. 

In  general,  an  atnospheric  Held  variable  (for  example,  pressure 
or  density)  which  is  linear  in  it?  deper4-nce  on  velocity  at  a  fixed 
point  (in  Eulerian  coordinates)  m*y  Te  nonlinear  along  a  Lagrangian 
trajectory  passing  through  that  pc.ri  For  example,  while  the  density 
may  fluctuate  with  a  small  fmplitu^e  at  a  fixed  point,  along  a 
trajectory  passing  through  that  po;nt  it  may  fluctuate  with  a  much 
larger  amplitude,  due  in  part  lO  the  large  variation  in  background 
density  and  pressure  experienced  during  vertical  motion.  By  a 
simple  extension  of  the  procedure  in  HT87,  it  is  possible  to  obtain  the 
total  airglow  response  to  a  general  gravity  wave  in  closed  form, 
which  can  then  be  used  to  quantitatively  demonstrate  the 
importance  of  the  higher-order  nonlinear  response  terms. 


It  is  tempting  to  think  that  the  large  nonlinear  response  in  the 
Lagrangian  frame  is  due  implicitly  to  the  steep  density  gradient  in 
the  Eulerian  frame.  However,  we  have  not  been  able  to  show  this 
explicitly.  Nevertheless,  we  believe  it  is  worthwhile  to  point  out  that 
the  nonlinear  response  can  be  shown  to  be  large,  whatever  frame  is 
chosen  for  the  calculation. 


2.  Airglow  Intensity  Fluctuations 

Since  we  focus  on  purely  dynamical  effects  which  are  largely 
independent  cf  the  particular  radiator,  we  will  discuss  solely  the 
02(*2l'*‘g)  emission  considered  by  Weinstock  and  by  HT87,  We  ignore 
quenching  and  all  other  losses  except  radiative  loss.  We  assume,  for 
the  sake  of  discussion,  that  this  state  is  produced  and  destroyed 
solely  through  the  reactions 

O  +  O  +  M  — O  )  +  M 

o,('r-,)-»Oj+hv. 

where  M  denotes  the  major  species,  O  is  a  minor  species  with  a 
layered  structure  in  the  atmospheric  region  of  interest  (lower 
thermosphere  and  upper  mesosphere)  and  k  is  the  reaction  rate.  O  is 
assumed  to  be  a  passive  tracer  of  the  atmospheric  motions.  The  local 
intensity  will  be  given  by  , 

I  -  i  n2  N  (1). 

Here  n  (N)  is  the  minor  (major)  species  number  density.  The  reactioft 
rate  k  is  known  to  be  temperature-dependent,  but  again,  for  the  sake 
of  clarity  and  tneviiy,  we  will  take  it  to  be  constant. 

If  we  separate  the  number  densities  into  an  unpenurbed  pan 
and  a  wave  induced  fluctuation 

n(x,y,z,t)  a  no(z)  n’(x,yA0 
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N(x,y,z,t)  =  No(z)  +  N’(x,y,z,t) , 

the  volume  emission  rate  profile,  I,  upon  linearization,  is  given  by 

I  =  *nX|l  +  2^  +  ^]  (2). 

However,  as  mentioned  in  section  1,  the  minor  species  number 
density,  n,  may  require  a  nonlinear  treatment  even  though  the  wave 
itself,  carried  by  N,  is  lineaf.  One  method  (Isler  et  al.,  1988)  is  to 
solve  for  n  order  by  order  via  a  perturbation  expansion  in  the 
velocity  field,  and  so  obtain  ni,  n2,  n3  and  so  forth.  This  allows  one  to 
calculate  1  to  any  order  in  the  velocity  field  and  thus  obtain  the 
column-integrated  brightness 

B=Jldz  (3) 

to  any  order  in  the  velocity  field.  The  results  show  that  the  higher- 
order  terms  are  important  and  do  not  cancel. 


HT87  maintained  that  the  dynamical  nonlinearities  in  I,  which  are 
local,  are  not  present  in  B,  which  is  height-integrated.  To  show  this, 
they  made  two  assumptions:  (i)  that  n/N  is  conserved  exactly  along 
a  parcel's  trajectory  and  (it)  that  the  essentially  nonlinear 
transformation  from  Eulerian  to  semi-Lagrangian  coordinates  can  be 
adequately  represented  by  a  first-order  Taylor  expansion,  given  the 
linearity  of  the  gravity  wave  carried  by  the  major  species.  As  to  the 
first  assumption,  along  a  Lagrangian  trajectory  the  continuity 
equations  can  be  written 

IdN^- 

N  dt  (4) 


(5). 


where  v  (v.^)  is  the  velocity  of  the  major  (minor)  species  and  P 
and  L  are  the  production  rale  and  loss  frequency  for  the  minor 
species.  This  implies 
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(6). 


Thus,  the  first  assumption  is  correct  if  production  and  loss  processes 
are  negligible  over  gravity  wave  time  scales  and,  in  addition,  the  two 
species  have  the  same  velocity  fields  at  all  times.  But  tlie  former  is 
only  true  for  some  airglow  reactions  (see.  for  example,  Walterscheid 
et  al.,  1987).  while  the  latter  ignores  the  diffusion  velocity  of  the 
minor  species,  inherent  in  any  hydrostatically  unstable  density 
profile,  such  as  a  layered  structure.  The  effects  of  such  diffusion  ar; 
not  negligible  during  the  transient  onset  of  the  gravity  wave,  nor 
during  some  fraction  of  each  gravity  wave  period  when  the  wave 
amplitude  is  small.  However,  for  the  remainder  of  this  paper,  we 
will  ignore  these  effects  of  diffusion  and  chemistry  and  like  HT87 
allow  that  (he  mixing  ratio  is  conserved  along  a  trajectory.  We  turn 
now  to  the  second  assumption,  namely  the  linearization  of  the 
Lagrangian  transformation. 


3.  Lagrangji3Ji.J:Qrmulau.Qn 


A.  The  Exact  Solution 


We  begin  by  writing  the  coupled  set  of  equations  expressir.g  mass 
and  momentum  conservation  and  the  adiabatic  equation  of  state  for 
the  major  species  in  the’  Eulerian  system  (using 
D  /  Dt  =  <?/ A  +  V  •  V): 


1  DN 
N  Dt 


N 


Dv 

Dt 


=  -V.v 


I 


(7) 

(8) 


1  Dp  iDN 
m  Dt  ^  Dt 

where  m  is  the  mean  molecular  mass,  g  is  the  acceleration  due  to 
gravity  and  c  is  the  speed  of  sound.  Here,  the  only  assumptions 
made  arc  that  m  is  constant  and  that  the  dissipative  loss  terms 

representing  kinetic  viscosity  and  heat  conduction  are  negligible. 
These  assumptions  are  generally  valid  in  the  mesosphere  and  lower 
thermosphere,  where  the  02(*£'*'g)  airglow  and  other  important 
airglow  emissions  originate.  Since  the  momentum  transferred  to  the 
minor  species  is  overwhelmingly  dominated  by  collisions  with  the 

major  species  and  since  we  neglect  chemistry  and  diffusion,  the 

minor  species  may  be  assumed  to  have  the  same  velocity  and 

acceleration  Helds  as  the  major  species,  so  that  v„^  =  v,  and  the 
continuity  equation  for  the  minor  species  reduces  to 

i^  =  -Vv  (10). 

n  Dt 


With  no  further  assumptions,  we  can  in  pnnciple  solve  equations 
(7)-<9)  for  the  exact  v  in  the  original  Eulerian  system.  This  allows  us 
to  transform  the  continuity  equations,  (7)  and  (10),  using  this  v  into 


the  Lagrangian  system.  Thus, 

N  dt 

1  dn 
n  dt 


Vv  (11) 

Vv  (12). 


We  next  integrate  these  equations  along  the  trajectory  of  a  fluid 
parcel,  from  time  t  s  0,  before  the  gravity  wave  arrives  when  the 
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parcel  is  at  rest  at  z  =  to  time  t,  when  the  parcel  is  moving  with 
velocity  v.  We  obtain 


(13) 

(14) . 


Equations  (13)  and  (14)  are  a  direct  consequence  of  the  mass 
conservation  equations,  (11)  and  (12).  Obviously, 


n  n,(C) 

N  N.(f)  ■ 


and  the  mixing  ratio  is  conserved  along  the  trajectory.  The 
exponential  term  in  these  expressions  is  the  Jacobian  of  the 
transformation  from  Lagrangian  to  Eulerian  coordinates  [see,  e.g., 
Dutton  (1986)1,  in  agreement  with  the  usual  form  in  which  the 
Lagrangian  continuity  equation  is  written. 


N  = 


I  9{%,y,z) 


N. 


Clearly,  the  Jacobian  in  equations  (13)  and  (14)  is  nonlinear  in  vin 
general. 


The  Eulerian  luminosity  profile,  Ie=  k  ne^  Ne  (where,  in  the  future, 
the  subscripts  E  and  L  will  denote  the  Eule<ian  and  Lagrangian 
systems,  respectively),  can  then  be  transfonhed  to  the  Lagrangian 
system  as  follows: 


lE=ifcn|NE->lL  = 


* 

k 

M'f 

Ln.(oJ 

k 

Ln,(Oj 

(15  a) 


(15  b) 
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where  the  same  expression  for  the  conservation  of  the  mixing  ratio 
yields  the  two  alternate  forms  (15a)  and  (15bV  Using  either  form  for 
II  ,  the  integrated  brightness,  Bl  ,  is  given  by 

Bl  =  I  ^  (^)N„  (C)exp^-3  J  V  •  vdt' Jj^l  ( 1 6). 

The  factor  after  the  exponential  term  is  the  Eulerian  dz  transformed 
to  Lagrangian  coordinates,  using  z  =  ^  -k  h,  where  h  is  the  vertical 
component  of  the  parcel  trajectory.  In  equation  (16)  and  from  now 
on,  we  assume  that  the  horizontal  fluid  displacement  is  small 
compared  to  the  horizontal  wavelength.  We  are  then  justified  in 
considering  only  the  vertical  trajectoiy.  This  was  implicitly  assumed 
by  HT87. 

This  is  the  total  integrated  airglow  response  to  a  general  gravity 
wave  velocity  field  for  the  three-body  reaction  considered, 
neglecting  only  the  effects  of  diffusion  and  chemistry.  Note  that,  of 
course,  the  column-integrated  brightness,  being  a  real  physical 
quantity,  should  be  independent  of  the  coordinate  system  used,  so 
that  Bl  =  Be  .  where  Be  is  the  Eulerian  brightness  integral.  However, 
since  the  transformation  from  Eulerian  to  Lagrangian  coordinates  is 
nonlinear  in  v,  one  cannot  assume  that,  in  a  series  expansion  in  v,  Bl 
will  be  equal  to  Be  term  by  term. 

‘  * 

B.  Approximation  Schemes 

We  begin  by  making  two  observations.  First,  from  equations 
(13)  and  (14)  U  is  obvious  that,  along  a  trajectory,  n  and  N  vary  as 
complicated  functions  of  the  trajectory.  For  a  monochromatic  simple 
harmonic  trajectory  with  frequency  o>,  n  and  N  are  periodic  with 
fundamental  frequency  to  and  all  the  higher  harmonics.  The  question 
is:  Under  what  conditions  are  the  higher-order  harmonics  negligible? 
Second,  we  note  that  any  approximation  scheme  may  violate  some 
conservation  principle.  Actually,  in  this  context,  there  are  only  two 


basic  conservation  laws  involved,  namely  the  mass  conservation  of 
the  minor  and  of  the  major  species.  From  them,  one  can  deduce  the 
conservation  of  the  mixing  ratio.  However,  the  converse  is  not  true; 
the  conservation  of  the  mixing  ratio  does  not  imply  that  the  species 
densities  are  each  individually  conserved.  In  HT87  the  conservation 
of  the  mixing  ratio  was  initially  assumed.  As  a  result,  their 
subsequent  Lagrangian  linearization  of  the  major  species 
automatically  forced  the  linearization  of  the  minor  species. 


Before  we  consider  the  effect  of  linear  gravity  waves  we  consider 
a  hypothetical  wind  or  v'ave  velocity,  v,  which  is  exactly  known  and 
is  caiiied  by  the  major  species.  Once  again,  we  neglect  chemistry 
and  diffusion  and  assume  that  v^i,K>r  =  v  so  that  the  behavior  of  the 
minor  species  is  determined  completely  by  (10). 

We  introduce  a  dimensionless  parameter  X  and  let  v-»Av  so  that 
the  magnitude  of  the  velocity  can  be  adjusted  by  varying  X. 
Substituting  Av  in  (10),  transforming  into  the  Lagrangian  system  and 
integrating  along  the  trajectory  of  a  fluid  parcel,  we  obtain 


(  '  ^ 
n  =  n,(^)exp|^-A  J  V  •  v  dt' 


(17). 


where  ^  is  the  vertical  coordinate  of  the  parcel  before  the  wind  or 
wave  arrives  at  time  t  =  0. 


Since  V  is  known  exactly,  n  from  (17)  is  a  power  series  in  X. 
Clearly  if  X  is  chosen  small  enough,  we  need  only  keep  the  Hrst-order 
term  in  X.  If  X  is  large,  on  the  other  hand,  we  need  to  use  many 
terms,  or  even  the  entire  scries  which  is  summable  in  closed  form. 
Thus  for  a  given  wind  or  wave  velocity,  Av,  the  Lagrangian  density 
response  is  highly  nonlinear  in  X. 
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In  the  case  of  gravity  waves  we  expand  the  dynamical  variables 
of  equations  (7)-(9)  in  a  power  series  in  a  dimensionless  parameter  X 
(X  here  does  not  mean  quite  the  same  thing  as  before).  Thus, 


p  =  p„+Ap,  +  A*pj+.. 
N  =  N„  +  ANi  +  A*Ni+. 


v  =  Av, +A^V2+- 


,  (18). 


where  we  have  assumed  (hat  the  velocity  is  entirely  v/ave-induced 
(no  prevailing  wir>.d),  and.  for  convenience,  we  choose  the  coordinate 
system  so  that  v=:(u,0,w).  Substituting  this  expansion  into  the 
dynamical  equations  and  solving  them  to  first  order  in  X  giver  the 
well-known  expressions  of  Hines  (1960)  for  v„  pi,  and  Ni.  Solving 
the  equations  to  X2  gives  Vj,  p2,  and  N2  in  terms  of  products  of  first- 
order  terms.  For  example,  V2  would  be  a  linear  combination  of  Nipi, 
NiVi,  Vi*  etc..  For  Iu|«Vpte,  where  Vphx  is  the  horizontal  phase  trace 
speed,  0)/kx,  one  finds  (Fritts,  1984), 


A*|vjJ«A|v,l 

A*|p2j«A|p,j 

A^|N,j«AjN,l 


(19). 


Here  we  can  see  that  the  power  series  expansion  in  X  given  by  the 
third  line  of  (18)  is  an  expansion  for  the  major  species  velocity  which 
has  to  satisfy  the  hydrodynamic  equations,  with  major  species 
background  density,  pressure,  etc..  •'  The  expansion  in  X  given  by  (17) 
is  an  expansion  in  velocity  for  the  lAinor  species  density  which  needs 
to  satisfy  only  the  continuity  equation  for  the  minor  species,  with  a 
minor  species  background.  The  two  are  very  different  expansions 
and  would  not  in  general  converge  at  the  same  rate. 

Through  (19)  and  the  continuity  equation  for  the  major  species  it 
is  easy  to  show  that 

A*|Vv2|«^jV*Vi|  (20). 

Thus,  using  v^Av, +A*V2+ . »  equation  (1 7)  becomes 
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In  the  second  line  of  (21),  we  see  that  for  each  order  of  velocity 
perturbation,  (v„  Vj,  etc.)t  the  response  of  the  minor  species  is  in  the 
form  of  an  exponential  function  of  that  particular  order:  a  very 
nonlinear  response.  In  this  paper  we  shall  discard  all  the  terms  of 
0(X2)  and  higher  in  the  exponent  of  the  second  line  of  (21),  leaving 
only  the  exponential  response  to  a  strictly  linearized  gravity  wave, 
i.e.  v  =  Av,.  The  same  approximation  was  made  by  Gardner  and 
Shelton  (1985)  when  they  also  considered  the  nonlinear  response  to 
a  linear  gravity  wave. 


To  see  that  there  can  be  no  cancellation  between  the  higher- 
order  teiiTts  in  the  nonlinear  response  to  Xv|,  and  a  corresponding 
higher-order  term  in  the  expansion  of  v  such  as  we  consider 

the  third  lino  in  (21).  Of  the  two  X?  terms  on  the  right-hand  side  of 
(21)  the  first  belongs  to  a  series  expansion  of  (17),  the  second 
belongs  to  the  velocity  expansion  in  (18),  and,  as  already  mentioned, 
the  two  series  do  not  converge  at  the  same  rate.  Furthermore,  v,  and 
Vj  are  linearly  independent  so  that  there  is  no  way  (for  arbitrary 
time  t)  that  the  two  X2  terms  can  cancel.  Indeed,  Owing  to  the 
oscillatory  nature  of  Vj  the  fact  that  the  term  involving  v,  is 
always  positive,  there  are  times  when  the  two  terms  would  reinforce 
rather  than  partially  cancel.  Whilst  for  the  larger  gravity  waves  V2 
may  well  be  important,  for  this  paper  we  shall  only  be  concerned 
with  a  iincit  gravity  wave,  i.e.  v  =  Av,  and  as  in  the  case  of  small 
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magnitltde  winds  we  shall  consider  the  nonlinear  response  to  a  linear 
velocity  perturbation  (bearing  in  mind  that  the  nonlinear  response 
cannot  in  general  be  cancelled  by  higher-order  terms  such  as  Vj). 
Our  numerical  results  seem  to  show  that  the  higher-order  corrections 
in  the  expansion  of  (17),  even  for  a  gravity  wave  with  a  very  modest 
amplitude  (say  a  4%  fluctuation  in  the  major  species  density,  which 
should  allow  (2G)  to  be  valid),  the  nonlinear  response  of  the  airglow 
is  important  (see  Table  1).  From  now  on  we  shall  omit  writing  X. 

It  is  customary  in  the  Hines  linear  gravity  wave  model  to  replace 
the  real  Helds  v,  ,  and  N]  by  complex  Helds  to  take  advantage  of 
the  simplicity  of  the  algebra  in  the  complex  notation.  This  is 
permissible  because  only  linear  transformations  are  carried  out  on 
the  variables.  It  is  tacitly  understood  that  the  real  part  is  to  be 

taken  at  the  end  of  the  calculation.  On  the  other  hand,  when 

nonlinear  transformations  are  carried  out,  as  in  the  present  case, 

greater  care  is  required  to  ensure  that  only  real  quantities  are  used 
in  nonlinear  transformations.  As  an  example,  we  note  that,  in 

squaring  a  sinusoidal  Held, 

cos*  oit  Re  e**®’ 

The  right  member  correctly  yields  the  second  harmonic  term,  cos2a)t, 
in  cos^,  but  completely  misses  the  dc  term.  In  our  case,  while,  we 
choose  to  keep  the  linear  Held  varijables  v,,  pj  and  Ni,  complex,  so 
that  the  usual  linear  theory  may  be  applied,  we  also,  make  sure  to 
take  the  real  part  of  the  field  variables  that  appear  either  in  an 
exponent  or  in  a  product  of  terms. 

In  the  linear  theory,  we  have  v,  =^d,/^t  instead  of  v,  =Dd,/Dt, 
where  dj  is  the  Eulerian  displacement  H:ld  (we  keep  only  the  linear 
terms  in  the  Eulerian  system;  signiHcant  effects  such  as  the  Stokes 
drift  (Coy  ct  al.,  1986)  have  been  neglected).  Thus,  when  the  velocity 


103 


div»*.rgence  is  evaluated  in  the  Eulerian  frame  prior  to  the  Lagratigian 
transformation,  we  may  set 


Vv, 


<9(Va,) 

<?t 


(22). 


Since,  unlike  the  major  species,  the  minor  species  obeys  only  the 
continuity  equation  00).  with  v  =  v,,  we  transform  this  equation  to 
the  Lagrangian  system  exactly,  keeping  all  orders  in  Vj.  Later,  we 
will  investigate  the  magnitudes  of  the  higher-order  terms  introduced 
by  this  transformation.  Thus,  using  (14)  and  (22),  after  transforming 
to  the  Lagrangian  system,  we  may  write 

Hl  =  n.(Ocxp^-Re£^^~^dt'j 

=  n,(Oexp[Re(-Vd,{,  +  Vd,|,^)] 

=  n,(C)cxplRci(Vv,)]  (23), 


where  in  the  second  step  we  have  used  the  fact  that  all  integrations 
are  partial  integrations,  and  in  the  last  step  we  have,  as  stated, 
neglected  any  effects  of  the  initial  diffusion  velocity  so  that  at  t  =  0, 
taken  to  be  before  the  onset  of  the  gravity  wave,  d,  *  0  for  all  x,  y 
and  z.  At  time  t,  we  have  assunr.ed  that  the  gravity  wave  has 
reached  a  steady  asymptotic  state  with  a  monochromatic  velocity 
Held  so  that  v,  siojd,.  The  exact  conservation  of  the  mixing  ratio  then 
implies 


NL  =  No(C)exp[Re(^Vv,)]  (24). 

With  these  expressiovis,  we  have  the  exact  conservation  of  each 
species  and  of  the  mixing  ratio.  The  column-integrated  brightness  of 
equation  (16)  becomes 
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BL=J*n;(f)N,(f)exp[Re(3iV-v,)]|l-Re(i^))df  (25'. 

0 

The  factor  after  the  exponential  reflects  the  fact  that  h(z,t),  in  steady 
state,  is  a  simple  harmonic  vertical  trajectory. 

It  is  important  to  stress  again  a:  this  point  that  N2L  n>ay  be 
significant  even  though  N2E  is  negligible,  and  this  is  due  to  the 
nonlinearity  of  the  transformation  relating  them.  In  fact,  while  Nl 
and  Ne  noth  refer  to  the  major  species,  they  are  actually  quite 
different  physically;  being  the  density  of  one  air  parcel  along  its 
trajectory  is  affected  by  background  variations  over  a  cycle,  while  Ne 
is  the  density  of  a  succession  of  air  pan  els  passing  through  a  fixed 
point  where  the  background  is  fixed.  So  there  is  no  reason  to  expect 
Nl  and  Ne  to  be  the  same,  or  even  to  converge  at  the  same  rate  in  a 
perturbation  expansion.  As  evident  in  equation  (24),  following  a 
closed,  simple  harmonic  trajectory,  NL=NL(d,,t)ssNL(-iv„t)  is  in 
general  nonlinear  in  v,,  even  though  at  ly  fixed  point  Ne  may  be 
linear  in  v,.  Again,  the  reason  is  the  transformation  from  the 
Eulerian  to  the  Lxigrangian  system,  that  is,  the  Jacobian  given  by  the 
exponential  in  equation  (24)  is  nonlinear  in  Vj. 

This  transfonnation  wa.<i  linearized  by  HT87  who  transformed  Nie 
(which  is  already  linear  in  v,)  to  the  Lagrangian  system  via  a 
linearized  transformation,  or  a  first-order  Taylor  expansion.  Thus, 
according  to  HT87. 
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(26) 


=  N.(C  +  h,)  +  N^E(C  +  h,.0 

dN 

-N.(0  +  ^h,  +  N,E(C)  +  O(V) 
dC 

-  N.(0  +  ^h,  +  i-N.V-  V,  +  -^w,  +0(V) 

<a  u> 

-N.(0(l>-^^  v,^(  ‘N.tN,, 

where  the  first  two  lines  are  nonlinear  in  V|,  white  the  last  three  are 
linearized.  In  going  ftom  the  third  line  to  the  fourth  line,  the  first- 
ojder  Eulerian  continuity  equation  for  a  sinusoidal  wave, 

ia>N,(  =-N.Vv, -WidN.fdC  (27) 

has  been  used. 

So  far  as  we  can  see,  the  rationale  behind  this  linearizadoM  is  that 
if  Ne  is  linear,  then  Hl.  being  still  the  major  species,  should  also  be 
used  in  its  linearized  form.  However,  Ne  at  a  fixed  point  is  a 
different  function  of  v  from  Nl  which  follows  a  trajectory  through 
the  same  fixed  point.  In  general  if  g^v)  and  f(v)  are  two  arbitrary 
functions  related  by  a  transformation  nonlinear  in  7,  the  fact  that 
g(v)  can  be  accurately  linearized  does  n^I  ensure  that  f(^)  can  also 
be  accurately  linearized.  Hence,  in  our  case,  the  linearity  of  Nl  dees 
not  follow  from  the  linearity  of  Ne. 

It  is  clear  that  equation  (26)  is  the  same  as  the  first-order  term  of 
equation  (24),  as  it  should  be.  If  we  adopt  this  form  for  Nl, 
conservation  of  the  mixing  ratio  would  force  ol  to  acquire  the  same 
form.  We  then  have  the  peculiar  case  of  an  exact  conservation  in  the 
mixing  laiio,  wi.ile  each  species  is  conserved  only  to  Hrst  order;  a 
point  we  have  mentioned  earlier.  In  addition,  if  we  now  write 
equation  (25)  to  first  order  in  7,,  we  obtain 

>0f» 


(28). 


BL  =  |*nl«:)N.(C)(l+3Re(iV.v,)|5-Re(i^)|5]dC 

This  is  just  the  expression  for  B  given  by  equation  (20)  of  HT87  (with 
a=l,  A=0).  (As  remarked  in  HT87,  this  expression  may  also  be 
obtained  directly  from  the  linearized  Eulenan  equation  (2).  after  an 
integration  by  parts).  The  linearization  of  the  coordinate 
transformation,  called  a  Lagrangian  linearization  by  H787,  has 
guaranteed  a  linear  response  for  B. 

To  summarize  this  section,  HT87  assumed  that  Nte*  being  a 
solution  to  the  linearized  Eulerian  hydrodynamic  equations,  could  be 
transformed  linearly  to  the  Lagrangian  system  with  no  loss  of 
accuracy.  Then,  through  the  consuncy  of  the  mixing  ratio,  (15s)  was 
employed  to  express  II  in  terms  of  Nil.  which  upon  further 
iinearization  gave  equation  (28).  We,  in  contrast,  use  (ISb),  which  is 
just  as  valid  as  (15a),  to  express  II  in  terms  of  nt.  Since  ng  satisfies 
a  single  equation  (the  continuity  equation)  exadtly,  in  the  Lagrangian 
system  the  solution  nt  must  satisfy  equation  (23).  With  the 
constancy  of  the  mixing  ratio  we  obtain  equation  (25).  If  the 
transformation  can  be  accurately  linearized,  as  HT87  assumed,  there 
should  be  no  significant  difference  between  nsing  either  equation 
(25)  or  (28).  We  now  explicitly  calculate  the  higher-order 
Lagrangian  terms,  to  determine  if  this  is  true. 


flJ  imB  t  wBI  iTiiiM 


To  first  order  in  V|  in  the  Eulerian  system,  using  the  polarization 
relations  of  Hines  (1960),  one  can  show  (Makhlouf  et  al.,  1989) 


N.»,(y-i)i». 

N.  y  H  c* 


(29). 


lb? 


Here,  y  is  the  ratio  of  specific  heats  and  H  is  the  atmGfpheric  density 
scale  height.  Combining  this  with  the  first-order  continuity  equation, 
equation  (27),  we  have 


Thus,  to  first  order  in  v,  in  the  Lagrangian  system,  using  equation 
(24). 


Nil 

N. 


hi 


(30). 


For  the  small  scale  waves  (vpha  «  c),  often  found  in  middle 
atmosphere  airglow,  we  may  neglect  the  second  term  on  the  right- 
hand  sides  of  (29)  and  (30)  apd  obtain 


N.  r-1  N.  N, 

where  we  have  used  y=  1.4.  This  demonstrates  that  the  Lagrangian 
density  fluctuation  is  bJT.er  than  the  Eulerian  by  a  factor  of  2.5..  For 
larger  scale  waves,  the  relation  is  not  sc  simple,  but  the  Lagrangian 
density  fluctuation  remains  on  average  significantly  greater  than  the 
Eulerian  density  fluctuation.  So,  in  an  obvious  expansion  of  It.  the 
integrand  of  equation  (25),  we  write 
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=  knJ{f)N,(f)x 


ju3Re(^V.v,]  +  i[Re(iV.v,] 


+  0(v,’) 


(32). 


Then,  wc  obiain  kom  equations  (30)  and  (31) 


(33) 


or,  dividing  equation  (34)  by  equation  (33) 

y  =  (35, 

First  of  all,  we  notice  from  (33)  that  the  Hrst-order  luminosity 
fluctuation  in  the  Lagrangian  system  is  7.5  times  the  Eulerian  major 
species  density  fluctuation;  a  10%  major  species  density  fluctuation 
can  produce  a  75%  brightness  fluctuaitioft.  This  is  very  similar  to  the 
purely  Eulerian  computations  in  which  the  first-order  brightness 
fluctuation  can  be  very  large  whilst  the  major  species  density 
fluctuation  is  only  modest.  Notice  that  this  factor  7.5  is  a  product  Of 
two  factors:  a  factor  of  2.5  from  equation  (31).  which  shows  that  the 
first-order  density  fluctuation  at  a  flxed  point  is  only  40%  of  the 
fluctuation  along  a  trajectory  passing  through  that  point;  and  a  factor, 
3,  which  reflects  the  conservation  of  the  mixing  ratio  along  a 
trajectory,  coupled  with  *he  assumption  that  the  airglovr  in  question 
is  produced  by  a  three-body  reaction. 
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Secondly,  we  now  that  for  a  10%  Euierian  density  fluctuation  the 
second-order  intensity  fluctuation  is  ‘31.5%  of  the  first  order,  and 
28.1%  of  the  background  --  clearly  not  insignificant.  This  agrees  witit 
the  purely  Eulerisn  calculations  (Islcr  et  al.,  1988),  in  which  the 
higher  orders  are  also  seen  to  be  important.  Table  I  shows  the  first- 
and  second-order  intensity  fluctuations  in  the  Lagrangiac  system 
against  given  Euierian  major  species  den.sity  fluctuaticrs.  Wc  see  in 
all  cases  that  the  second -order  intensity  fluctuation  is  greater  than 
the  first-order  density  fluctuation. 


labt^L 


j  i!y 

N. 

ly 

1. 

!kJ 

ivl 

y 

I. 

4% 

30% 

15% 

4.5% 

6% 

45% 

22.5% 

10.1% 

8% 

60% 

30% 

18% 

10% 

75% 

37.5% 

28.1% 

Finally,  while  the  two  points  mentioned  above  deal  with  local 
intensities  and  densities,  the  integrated  brightness  wM..h  is  being 
measured  will  have  the  same  numerical  factors,  and  would  hence  be 
expected  to  remain  important.  We  can  convince  ourselves  of  the 
importance  of  the  nonlinear  response  terms  by  looking  at  the  column 
brightness.  In  Fig.  1-4,  wc  have  used  the  Hines  gravity  wave  model 
which  satisfies  the  rigid  surface  ground  boundary  condition.  The 
unperturbed  density  profiles  (O,  N2.  etc.)  were  taken  from  the  U.  5- 
Standard  Atmosphere  (1976).  Two  types  of  diagrams  have  been 
plotted.  One  (Figures  1  and  4)  deals  with  the  variation  of  brightness 
as  a  function  of  horizontal  phase  trace  speed,  Vphx«  at  a  fixed  time 
and  period.  Oscillations  occur  because  of  the  varying  wavelengths 
produced  by  varying  vphx  and  the  fact  that  the  ground  boundary 
condition  ensures  that  all  waves  begin  with  the  same  phase.  The 
second  type  of  figure  (Figures  2  and  3)  illustrates  the  variation  in 
brightness  versus  time  at  fixed  Vphx 
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Here  we  should  mention  that  we  assume  the  background  to  be 
approximately  uniform  so  far  as  the  gravity  wave  is  concerned  all 
the  way  up  to  and  including  the  airglow  region  under  discussion. 
The  general  gravity  wave  in  the  airglow  region  can  then  be  a  linear 
combination  of  an  upwardly  travelling  Hines  gravity  wave  and  a 
vertically  stationary  non-evanescent  Hines  gravity  wave.  The  latter 
provides  for  the  large  downward  reflection  from  the  sudden 
temperature  rise  at  the  base  of  the  thermosphere  for  larger  scale 
waves.  For  smaller  scale  waves  there  can  sometimes  be  significant 
downward  reflection  from  critical  layer  regions  in  the  actual 
atmosphere  (He  et  a!.,  1990).  Since  the  primary  purpose  of  the 
present  paper  is  to  demonstrate  the  importance  of  nonlinear  airglow 
response,  even  for  relatively  small  amplitude  gravity  waves,  we  use 
only  the  Hines  gravity  wave  satisfying  the  rigid  surface  boundary 
condition;  in  the  Eulerian  system  a  perturbation  expansion  for  this 
particular  type  of  wave  does  not  require  renormalization  and  is 
therefore  simpler  to  use. 

In  Fig.l,  we  show  a  plot  of  the  first-order  Lagrangian  and  first- 
order  Eulerian  brightness  terms  (which  are  the  same,  of  course, 
because  of  the  relation  via  integration  by  parts  referred  to  above)  at 
a  fixed  time  but  for  different  horizontal  phase  trace  speeds,  for  a 
two-hour  period  gravity  wave  that  induces  a  5%  major  species 
density  fluctuation  at  97  km.  In  the  same  figure  we  h^ve  also 
plotted  the  second-order  Lagrangian  (dashed  line)  and  second-order 
Eulerian  (dotted  line)  terms.  Clearly  both  are  significant  for  ali 
horizontal  pha::e  trace  speeds.  In  Fig.2  we  plot  the  first  and  second . 
order  terms  for  a  fixed  Vphx  -  ISO  m/s  as  a  function  of  time,  where 
the  maximum  density  fluctuation  of  the  major  species  is  again  set  to 
5%  at  97  km. 

Fig.  3a-e  show  a  comparison  of  the  total  Lagrangian  brightness 
(equation  (25),  crosses),  as  a  function  of  time,  and  the  sum  of  the 
first  three  perturbation  terms  in  t.he  Eulerian  brightness  (doited 
lines)  well  as  the  first-order  brightness  (equation  (28),  solid  lines), 
using  the  same  parameters  as  Fig.  2.  Each  figure  i$  drawn  for  a 


different  fractional  rr^ijor  species  density  fluctuation  at  97  km, 
beginning  with  7.5%  in  Figs.  3a  through  5%,  3J5%,  2.5%  and  ending 
with  0.5%  for  Fig.  3e.  As  can  be  seen,  the  nonlinear  response, 
producing  a  cusp  at  the  waveform  crest  and  a  flattening  Of  the  wave 
trough,  is  very  important  for  the  7.5%  case,  and  almost  completely 
absent  for  0.5%.  Figs.  4a-e  show  the  brightness  response  at  a  fixed 
time  but  for  different  horizontal  phase  trace  speeds;  tfie  importance 
of  the  higher  orders  is  clear. 


5.  Conclusion 


The  question  of  whether  a  linear  gravity  wave  will  give  rise  to 
higher-order  harmonics  of  the  basic  gravity  wave  frequency  in 
ground-based  observations  of  airglow  brightness  fluctuations  is 
important  for  the  investigation  of  atmospheric  dynamics,  especially 
with  regard  to  the  proper  interpretation  of  the  power  spectra  of 
airglow  data.  That  such  harmonics  are  significant  has  been  shown  to 
be  the  case  when  they  are  calculated  in  the  Eulerian  and  in  the 
Lagrangian  frames. 

The  net  consequence  of  all  this  is  that  one  would  expect  the 
frequent  occurrence  of  higher-order  harmonics  in  the  power  spectra 
of  airglow  brightness  fluctuations.  The  results  presented  here  may 
be  tested  by  the  simultaneous  measurement  of  ground-based 
brightness  fluctuations  and  local  major  species  density  fluctuations 
via  radar  or  lidar  or  in  situ  instruments.  Such  an  experiment  could 
determine  whether  higher  harmonics  present  in  the  nonlinear 
airglow  response  are  in  fact  absent  from  the  actual  gravity  wave 
fields. 
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Figure.  .Captigns 


Figure  1.  Comparison  of  the  second-order  Eulerian  and  Lagrangian 
integrated  brightness  (with  a  second-order  dc  offset  omitted)  to  the 
first-order  integrated  brightness  (eqn.  (28)),  at  a  fixed  time  for  a  two 
hour  period  wave  over  a  range  of  horizontal  phase  trace  speeds 
(vphx)*  The  fractional  major  species  density  fluctuation  (Nje/No)  at 
97  km  has  been  set  to  5%.  The  vertical  scale  is  given  in  arbitrary 
units. 

Figure  2.  The  same  as  in  Figure  i,  but  for  a  fixed  Vphx  of  150  m/s, 
plotted  as  a  funciion  of  time.  The  fractional  major  species  density 
fluctuation  at  97  km  has  been  set  to  5%.  The  vertical  scale  is  given  in 
arbitrary  units. 

Figure  3.  Comparison  of  the  sum  of  the  first  three  orders  (including 
the  dc  offset)  of  the  integrated  brightness  in  the  Eulerian  system  and 
the  total  Lagrangian  integrated  brightness  (eqn.  (25))  to  the  first- 
order  integrated  brightness  (eqn.  (28)),  for  a  two-hour  period  wave 
with  a  horizontal  phase  trace  speed  of  150  m/s.  The  fractional  major 
species  density  fluctuation  at  97  km  is  7.5%  in  Fig.  3a,  5%  in  Fig.  3b, 
3.75%  in  Fig.  3c,  2.5%  in  Fig.  3d  and  0.5%  in  Rg.  3e.  The  vertical  scale 
is  given  in  arbitrary  units. 

Figure  4.  The  same  as  in  Figure  3,  but  at  a  fixed  time  and  ever  a 
range  of  horizontal  phase  trace  speeds.  The  vertical  scale  is  given  in 
arbitrary  units. 
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Figure 


speeds  (Vphx) •  The  fractional  iiajor  species  density  fluctua 
((jjg/No)  at  97  km  has  been  set  to  5%.  The  vertical  scale  is 
in  arbitrary  units. 


Figure  3a 
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Figure  3  -  ^omparisoa  of  the  sum  threeordors  (includiiKj  the  dc  offset)  of  the 

integrated  brightness  in  the  Eulcrian  system  and  the  total  Lagrangian  integrated  brightness 
(eqn.  25))  to  the  first-order  integrated  brightness  (eqn.  (2B)),  fi.r  a  two-hour  i>et  iud  wave 
with  a  ))orizuntal  phase  trace  speed  of  150  m/s.  The  fractional  maior  species  density  fluctuation 
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Figute  4c 


vphx  (M/t*c) 

Fi9ure  4  -  The  sisme  as  in  Figure  3,  but  at  a  fixed  tijne  and  o</er  a  range 
of  horizontal  phase  trace  speeds.  The  vertical  scale  is  given  in 
arbitrary  units. 
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Abstract 


The  local  fluctuations  of  a  minor  atmospheric  species  roncentiation  caused  by 
linear  gravity  wave  motions  can  be  nonlinear  if  the  minor  species  density  gradient 
is  sufficently  steep.  To  treat  such  a  case,  a  method  is  outlined  by  which  it  is 
possible  to  calculate  the  minor  species  response  to  any  order  in  the  linear  gravity 
wave.  Calculations  to  third  order  over  a  wide  range  of  wave  parameters  show  that 
the  nonlinear  effects  can  be  substantial.  As  a  result,  caic  must  be  taken  when 
analysing  data  from  minor  species  fluctuations,  so  that  frequencies  due  solely  to 
the  nonlinear  nature  of  the  minor  species  respond  are  not  attributed  to  gravity 
waves. 
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1.  Introduction 


The  study  of  gravity  waves  in  the  middle  and  upper  atmosphere  has  relied 
from  its  beginning  upon  the  detection  of  the  motion  of  minor  atmospheric  con¬ 
stituents.  Early  examples  include  deformations  of  meteor  trails  and  radio  echoes 
from  ionospheric  layers.,  while  more  recently  many  observations  have  been  made  of 
airglow  intensity  fluctuations  and  resonance  lidar  backscatter, fluctuations. 

Many  of  the  minor  constituents  of  the  atmosphere  are  distributed  in  layers 
which  are  sharply  peaked.  The  atomic  sodium  layer,  for  example,  is  concentrated 
around  90  km  in  a  layer  which  has  on  average  a  half-width  much  less  than  the 
atmospheric  scule  height.  Because  of  such  background  structures,  small-amplitude 
(i.e.  linear)  gravity  waves  can  induce  local  fluctuations  of  much  larger  amplitude  in 
a  minor  species  concentration.  Modelling  the  response  of  minor  atmospheric  species 
to  gravity  wave  perturbations  has  been  carried  out  by  many  investigators  (Thome, 
1968;  Porter  and  Tuan,  1974;  Chiu  and  Ching,  1978;  Weinstock,  1978;  Gardner  and 
Shelton,  1985;  Molina  et  al.,  1985;  Walterscheid,  et  al.,  1987;  Hines  and  Tarasick, 
1987;  Hickey,  1988;  Tarasick  and  Shepard,  1989).  With  the  exception  of  Weinstock 
(1978)  and  Gardner  and  Shelton  (1985),  previous  studies  have  focused  on  the  linear 
response.  The  present  work  develops  a  perturbation  expansion  approach  that  can 
be  used  to  model  the  nonlinear  response  of  neutral  minor  species  to  linear  gravity 
waves.  The  approach  taken  is  Eulerian,  so  that  the  method  has  direct  application 
.  to  measurements  made  at  a  fixed  height.  It  is  apparent  that  with  the  inherent 
reliance  of  so  many  remote  sensing  technologies  (ionosondes,  airglow  photometers, 
radars,  lidars)  on  the  detection  of  minor  species  motions,  a  clear  understanding 
of  the  response  of  minor  atmospheric  layers  to  linear  gravity  waves  is  essential  for 
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the  proper  interpretation  of  such  data,  particularly  when  the  data  are  spectrally 
cinalysed. 

The  present  paper  confines  itself  to  the  treatment  of  the  local  minor  species 
response,  as  measured  at  a  given  height  by  active  instruments  such  as  radars  or  res¬ 
onance  iidars.  A  companion  paper  addresses  the  miner  species  response  integrated 
along  a  iine-of-sight,  as  measured  by  passive  instruc’ents  such  as  ground-based 
photometeri 

Section  2  displays  the  origin  of  the  nonlinear  dynamical  response  of  minor 
species  with  sharply  peaked  profiles.  Section  3  develops  the  perturbation  expansion 
method  in  general,  and  in  section  4  it  Is  used  to  calculate  higher-order  terms  in  the 
minor  species  dynamical  response  over  a  broad  range  of  wave  parameters  as  well 
as  for  an  illustrative  specific  case.  Conclusions  are  presented  in  section  5. 


2.  I'^*^amical  Nonlinearity 

I.in'rar  gruvitv  wave  tlu’ory  has  on  the  wImU*  h<vii  drvelcj)od  for  the  special 
case  of  a  homogeneous  fluid,  where  the  mean  molecular  mass  is  constant  and  mixing 
processes  overwhelm  difFiisive  separation.  This  is  in  general  a  good  approximation 
below  the  mesopause;  near  and  above  100  km,  however,  the  mean  molecular  mass 
becomes  strongly  variable  with  increasing  altitude,  and  at  the  same  time  the  mean 
free  path  is  large  enough  for  species  to  separate  diffusively.  In  the  thermosphere, 
therefore,  even  for  neutral  minor  species,  a  proper,  treatment  of  the  response  to 
gravity  wave  perturbations  mu.st  incorporate  multi-component  fluid  effects  and 
have  separate  moment ;-ni  iqtiations  for  each  component  (Volland,  1969;  Mayr  et 
al.,  1984).  In  Use  piesent  paper,  attention  is  restricted  to  neutral  minor  species  in 
the  stratosphere  and  mesosphere.  In  those  regions,  where  minor  spedes  concen¬ 
trations  never  exceed  1%  of  the  total  atmospheric  concentration,  one  is  justified 
in  considering  only  the  continuity  equation  for  the  investigation  of  neutral  minor 
species.  This  is  because  the  overwhelming  majority  of  minor  species  collisions  are 
with  major- species  particles,  and  one  can  assume  that  they  have  the  same  velocity 
as  the  major  species,  except  possibly  for  diffusion,  which  we  will  consider  next. 

The  Eulerian  continuity  equation  for  any  minor  species  is 

^  =  -V.(nf?„)  +  /^-n£  (1) 

where  n{i,y,z,t)  is  the  minor  species  number  density,  Vn«(x,  y,  r,t)  is  the  minor 
species  velocity  field,  P  is  the  minor  species  production  rate  per  unit  volume,  and 
L  is  the  minor  species  loss  frequency.  Since  we  want  our  tlicory  to  apply  to  airglqw 
emission  as  well  as  active  probing,  we  will  take  the  term  minor  species,  in  such 
cases,  to  mean  the  particular  electronic  or  vibrational  state  of  the  species  which  is 
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the  source  of  the  airglow.  We  wish  to  use  this  equation  to  examine  the  response  of  a 
minor  species  to  dynamical  perturbations  induced  by  a  gravity  wave,  in  particular 
the  nonlinear  response  to  such  perturbations. 

Other  processes  which  must  be  considered,  in  addition  to  gravity-wave-induced 
djTiamical  perturbations,  are  diffusive  processes  and  chemical  (Including  collisional 
and  radiative)  production  and  loss.  Even  in  the  absence  of  any  motion  of  the 
atmosphere  as  a  whole,  it  may  be  that  Un,  ^  0.  This  is  because  any  minor  cp^’ies 
with  a  layered  structure  is  in  diffusive,  not  hydrostatic,  equilibrium,  and  is  therefore 
continually  diffusing  from  regions  of  net  production  to  regions  of  net  loss.  The 
vertical  diffusive  velocity  for  an  isothermal  atmosphere  is  given  by  (Chamberlain 


and  Hunten,  1S87) 


1  dno  . 


where  D  is  the  diffusion  coefficient,  m,  the  molecular  mass  of  the  minor  species,  g 


is  the  acceleration  due  to  gravity,  k  is  Boltzmann’s  constant,  and  T  is  the  temper¬ 


ature. 


If  there  is  ptir<*  gravity  wave  motion  (i.e.  no  prevailing  wind)  with  parcel 
velocity  v,  so  long  as  >  |u'ei//|t  one  can  s«*t  w,,  =  v.  However,  for  small  scale 
waves  it  may  be  that  |t^  ss  |w4,//|,  and  ever  for  somewhat  largetr  scale  waves  the 
diffusive  velocity  may  be  comparable  to  the  gravity  wave  velocity  for  an  appreciable 
part  of  each  gravity  wave  period,  as  v  passes  through  zero  twice  per  cycle. 

In  addition  to  this  diffusive  velocity,  ixrturbations  to  the  production  and  lorn 
terms  induced  by  the  gravity  wave  motion  may  be  comparable  to  the  other  per¬ 
turbation  terms.  If  this  is  (he  case,  the  effects  of  chemistry  will  be  as  important  as 
dyi'amic.J  effc-'t.s  in  drterinining  minor  s|>ecies  density  fluctuation^.  One  measure 
of  whether  tips  is  so  ( Waiterscheid  el  al.,  1987;  FVitts  and  Thrane,  1989)  is  if  the 
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diffusion  time  (or  the  lifetime  of  an  excited  state)  is  much  greater  than  the  v/ave 
period,  Tjiff{cktm)  ^  T.  The  range  of  diffusion  times  (lifetimes)  of  minor  con¬ 
stituents,  from  seconds  to  days  (Winick,  1983),  encompasses  the  gamut  of  gravity 
wave  periods,  which  are  from  5, or  6  minutes  up  to  several  hours.  Thus,  whether 
or  not  perturbations  to  the  production  chemistry  should  be  incorporated  in  the 
solution  to  equation  (1)  must  be'  decided  on  a  case  by  case  basis.  For  example,  for 
the  OH(i/')  levels,  which  are  the  source  of  the  Meirel  bands,  Tektm  ~  I  sec,  and 
chemistry  must  be  included.  For  atomic  sodium,  on  the  other  hand,  r^i  fy  >  1  day, 
so  for  this  spsdes  chemistry  may  be  safely  neglected. 

Since  the  focus  of  this  paper  is  on  dynamical  effects,  we  will  arbitrarily  assume 
i)  v,„  =  V  to  a  good  approximation,  that  is  vV,//  =  0,  and  ii)  perturbations  to  the 
production  and  loss  terms  are  negligible.  Under  these  assumptions,  equation  (1) 
becomes 

~  =  -V.(ntr),  (3) 

where  v  is  a  known  quantity,  and  (3)  is  therefore  an  uncoupfed  equation  for  n  only. 
To  simpU^  the  analysis,  for  the  remainder  of  this  paper  we  shall  take  this  v  to  rep¬ 
resent  a  monochromatic,,  linear  gravity  wave  srith  angular  frequency  u,  horixontal 
wavenumber  k^  .  vertical  wavemunber  k,,  propagating  in  a  srindless,  isothermal, 
horizontally  stratified  atmospiiere  with  components  (u,O,t0)  in  suitably  chosen 
Cartesian  coordinates,  with  z  and  t  dependence  given  by  exp  i{ut  —  kgz).  Such 
a  wave  is  a  solution  to  the  linearized  hydrodynamic  equations  of  the  major  atmo¬ 
spheric  species,  implying  that  the  major  species  number  density  (hereafter  denoted 
by  N),  is  fluctuating  with  a  small  amplitude  about  an  unperturbed,  background 
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value.  Given  this,  we  can  write  the  relative  fluctuation  in  the  major  specie*  ntunber 
density  directly  from  the  linearized  continuity  equation  as 


N’  i,  i  dw  t  1  dNo 

-  :=  - u  4- - 4-  —  —  W 

No  tAt  dz  ut  No  dz 


(4) 


where  N'{x,  z,  t)  =  perturljed  part  of  N  and  Noiz)  =  unperturbed  part  of  N. 

This  wave  will  induce  a  fluctuation  in  the  minor  species  density  (hereafter 
denoted  by  n).  However,  it  is  important  to  note  that  the  minor  species  fluctuation  is 
not  necessarily  a  “wave”  in  the  same  sense  as  the  major  species  fluctuation  is  a  wave. 
While  the  latter  is  a  solution  to  the  entire  set  of  linearized,  coupled  hydrodynamic 
equations  for  a  fluid,  the  former  is  the  response  of  a  minor  constituent  embedded  in 
the  fluid  to  wave-induced  perturbations  and  is  the  solution  to  just  the  uncoupled 
continuity  equation  given  by  (3)  for  the  minor  species  density  n.  This  is  true 
because  the  minor  species  momentum  equation  is  overwhelmingly  dominated  by 
the  collision  term  representing  momentum  transfer  from  particles  of  the  major 
species.  In  a  fashion  quite  analogous  to  the  derivation  of  equation  (4)  for  the 
major  species,  linearizing  equation  (3)  gives  the  first-order  approximation  to  the 
minor  species  response 


TXo 


t  dw  i  1  drto 

—  —  —  4-  _  _ — Iff 

W  Oz  W  Tlo  dz 


(5) 


where  n'(j,  z,  t)  =  perturbed  part  of  n  and  n,(r)  —  unperturbed  part  of  n. 

Two  points  to  consider  about  the  minor  species  response,  evident  in  the  first- 
order  equation  (5).  are  the  following.  First,  for  minor  species  with  a  layered  struc¬ 
ture,  the  unperturbed  density  gradient  is  positive  in  the  region  below  the  peak 
and  negative  in  the  region  above  the  peak,  while  the  unperturbed  major  species 
gradient  is  always  negative.  This  implies  a  shift  in  the  phase  difference  between  the 


136 


major  and  minor  species  density  fluctuations  at  the  point  of  the  minor  species  peak. 
Indeed,  for  a  wide  range  of  gravity  waves,  equations  (4)  and  (5)  are  dominated  by 
their  third  terms  (Makhlouf,  ct  al.,  1087),  so  that  below  the  peak  the  first-order 
minor  species  response  is  180  degrees  out  of  phase  with  the  major  species  fluctua¬ 
tion,  while  above  the  peak,  it  is  in  phase  with  the  major  species.  This  phase  reversal 
in  the  region  of  the  profile  peak  has  been  known  for  some  time  on  both  theoretical 
(Thome,  196S;  Porter  and  Tuan,  1974;  Weinstock,  1978)  and  observational  grounds 
(Thome,  1968;  Noxon,  1978;  Shelton  et  al.,  1980). 

A  second  point  to  consider  about  ,  equation  (5)  is  its  magnitude  relative  to 
equation  (4)  For  the  case  of  pure  linear  gravity  wave  motion  which  we  are  ton- 
sidcring,  the  magnitude  of  the  major  species  fluctuation,  given  by  equation  (4),  is 
by  assxunption  small,  say  10%  or  less.  However,  when  the  same  velocity  field  is 
used  in  equation  (5),  the  magnitude  of  the  first-order  mirxx^  species  response  can 
be  much  greater.  This  n»ult  ocevirs  because  the  advective  part  of  the  Eulerian 
time  derivative,  which  produces  the  third  terms  of  equations  (4)  and  (5),  involves 
gradients  of  backgroimd  densities.  As  mentioned  previously,  suck  gradients  can 
be  much  larger  for  minor  species  than  the  corresponding  major  species  gradient. 
Consequently,  though  equation  (4)  is  by  assumption  valid  for  the  major  species, 
retaining  only  first-order  tcrm.s  for  the  minor  species,  as  in  equation  (5),  is  not  in 
general  justified.  As  an  explicit  example,  if  one  uses  the  U.S.  Standard  Atmosphere 
(USSA)  1976  model  of  atomic  oxygen,  at  90  km  (just  below  the  O  peak  at  97  km) 
a  gravity  wave  which  ccuresponds  to  a  10%  density  fluctuation  in  the  major  species 
will  induce  al  40%  first-order  fluctuation  in  the. atomic  oxygen  density.  Thus  dy¬ 
namical  nonlinearities  enter  naturally  into  the  minor  species  response.  In  the  next 
section,  a  method  is  developed  to  treat  these  dynamical  nonlinearities. 


3.  Perturbation  Expansion  Method 

To  calculate  higher-order  solutions  to  equation  (3),  we  begin  by  rewriting  it 
(Isler,  et  al.,  19S8)  as 


where 


^  -h  +  AS0  +  A5 


xj)  =  n'/n,, 


(6) 


R  =  -w, 


/  du  dw  1  djio  \ 


and  we  have  simply  used 


1  tfh'  n'  drta 


dz  n«  dz  n\  dz 

in  rewriting  equation  (3).  Note  that  we  have  rewritten  the  exact  equation  (3), 
retaining  all  orders  in  the  wave  velocity  iT.  A  dimensionless  parameter,  A,  has  bien 
introduced  as  a  bookkeeping  device  to  keep  track  of  orders  in  tT,  or  equivalently  of 
orders  in  N'/N„  since  u  and  w  arc  linearly  related  to  Next,  we  write 


0  =  A0I  +  +  A’03  +  . . . 


(7) 


where  xpi  is  linear  in  N'/No,  ^7  i«  quadratic  in  etc..  Inserting  this  expansion 

into  equation  (6),  we  can  proceed  to  solve  for  the  minor  species  response  to  a  linear 
gravity  wave  to  any  order  in  A,  that  is  in  N'/N,. 

To  first  order  in  A,  equation  (6)  is  Otj^i/dt  =  5,  and  so 


(*  ( du  dw  1  dn„  \  , 


(8) 
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Since  we  are  interested  in  fluctuations  at  some  time  t  Ion/;  after  the  gravity  wave 
has  settled  into  its  monochromatic  state,  we  take  the  initial  time  in  equation  (8), 
t  —  to,  to  be  a  time  before  gravity  wave  motion  begins,  when  tT  =  0.  We  assume 
that  the  transient  state  develops  smoothly  and  continuously  into  a  monochromatic 
state.  As  it  should  be,  0i  obtained  in  this  way  is  equivalent  to  equation  (5).  Using 
the  above  notation,  tj^i  =  Re  —iS/ui, 

When  calculating  higlier-order  solutions  to  equation  (6),  it  must  be  remem¬ 
bered  that  the  real  parts  of  complex  quantities  must  be  taken  prior  to  any  multi¬ 
plication;  otherwise,  difference-frequency  effects  can  be  missed  altogether.  Thus, 
to  second  order  in  A,  equation  (6)  is  given  by 

^  =  ReQ  Re^  +  ReR  Re^  +  ReS  ReV-i  (9) 

at  dx  dz 

which  can  be  solved  since  is  known.  Once  agaiii,  the  integration  is  chosen  to  be 
from  a  time  prior  to  the  gravity  wave  motion. 

In  general,  at  the  nth  order,  we  have  the  recursion  relation 

^  =  ReQRc^^-»-ReRRe^^^-|-ReSReV'»-,.  (10) 

In  this  way,  one  can  obtain  0,  the  minor  species  response,  to  any  order  in  the  linear 
gravity  wave  parameters. 

.  Gardner  and  Shelton  (1985)  devised  a  similar  perturbation  expansion  ap¬ 
proach;  they  too  neglected  production  and  loss  perturbations  and  diffusibn,  and 
began,  m  we  have  done,  from  equation  (3).  However,  their  approach  assumed  a 
form  for  the  solution  in  advance,  and  they  considered  only  large  peiiod  gravity 
waves,  while  the  approach  taken  here  is  more  general.  , 


It  may  be  suggested  that  to  be  consistent,  one  would  heve  to  retain  all  nth 
order  terms,  and  so  in  equation  (6)  put  in  for  v  the  expansion 

V  =  At?i  +  A^t72  +  A^vs  +  —  (11) 

These  higher  orders  in  v  would  be  the  result  of  retaining  higher  orders  in  the 
solutions  to  the  coupled  hydrodynamic  equations  of  the  major  species,  i,e.  incor¬ 
porating  nonlinear  gravity  wave  effects.  However,  we  have  assumed  that  the  minor 
species  fluctuations  are  due  to  perturbation  by  a  linear  gravity  wave,  that  is  v  =  vi 
to  a  good  approximation,  so  we  further  assume  thsd  any  contributions  to  the  minor 
species  response  due  to  higher^order  corrections  to  v  are  negligible  relative  to  those 
due  to  the  linear  v  alone.  This  was  also  implicitly  assumed  by  Gardner  and  Shelton 
(1985). 
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4.  Results 


Wr  havr  «:alculat(^<l  tin*  minor  sjMvic's  n*sj)ons<*  to  tliird  order  in  the  linear 
gravity  wave  for  several  different  background  profiles  and  over  a  wide  range  of 
gravity  wave  parameters.  As  one  would  expect,  the  importance  of  the  higher-order 
terms  in  the  response  depends  strongly  upon  the  steepness  of  the  vertical  density 
gradient  of  the  minor  species.  Specifically,  we  have  tised  various  Chapman  function 
and  Gaussian  analytic  models,  and  the  USSA  197C  model  of  atomic  oxygen  as  well. 
Figures  la-lc  show  the  profiles  used  in  obtaining  the  subsequent  plots:  namely,  a 
Chapman  function,  a  Gaussian  distribution,  and  the  USSA  1976  atomic  oxygen 
model.  For  ca.sy  comparison,  the  two  analytic  profiles  have  each  been  set  to  peak 
at  92  km. 

We  have  considered  two  gravity  wave  models-  The  first  is  a  free  travelling 
wave  in  an  infinite  isothermal  atmosphere  (Hines,  1960),  with  velocity  components 

u  =  (12) 


«;  =  Au;(u;*-Jb*c*)e*'2"c'^‘-‘-*-*'*>  (13) 

where  H  is  the  atmospheric  scale  height,  c  is  the  speed  of  sound,  A  is  an  amplitude 
factor,  and  q  =  — (1  —  7/2)g/c*,  with  7  the  ratio  of  specific  heats.  The  second  is  a 
vertically  stationary  wave  in  a  semi-infinite  atmosphere  with  a  rigid  surface  ground 
boundary,  havine  velocity  comiamcnts 


^cos  k,z  +  ^  sin  ‘ 

=  ia sin  k,if  ^  7  , 


u»  =  losin — 5 — lL_ 

VwJ-wV 

where  ut  is  the  Brunt- Vaisala  frequency  and  a  is  another  amplitude  factor.  Such 
waves  come  from  either  a  high  altitude  source,  or  a  once7reflected  wave  from  a  low 


altitude  source,  or  a  guided  mode;  in  general,  an  arbitrary  gravity  wave  is  a  linear 
combination  of  (12),  (13)  and  (14),  (15).  The  detailed  algebra  leading  to  xjj2  and 
03  for  each  case  will  not  be  presented  here. 

The  stability  of  the  perturbation  method  outlined  above  is  dependent  upon 
the  type  of  gravity  wave  model  used;  in  particular,  secularity  (divergence  in  time) 
is  present  in  the  response  to  the  travelling  wave,  but  is  absent  in  the  response  to 
the  vertically  stationary  wave.  The  absence  of  secularity  in  the  latter  case  is  due 
to  the  fact  that  the  horizontal  and  vertical  velocity  fields  are  90"  out  of  phase, 
so  that  at  the  second  order  the  time  integral  contains  a  puK  elliptically  polarized 
sinusoid.  For  the  travelling  wave,  on  the  other  hand,  the  time  integral  contains  an 
offset  sinusoid,  and  secularity  result;..  It  is  interesting  to  note  that  the  growth  rate 
of  the  instability  is  proportional  to 


+ 


1  <fno\ 
rio  dz  }' 


(16) 


Thus,  i)  the  instability  disappears  if  the  minor  species  density  falls  off  with  the 
atmospheric  scale  height,  and  ii)  the  instability  is  more  pronounced .  for  small- 
scale  waves  (w  C  fc,c)  than  for  larger-scale  waves,  as  one  would  expect,  since 
the  relative  importance  of  vertical  motion  is  greater  for  the  former  than  foi;  the 
latter,  and  it  is  vertical  motion  which  produces  the  nonlinear  resjponsc.  To  correct 
for  the  secularity,  we  are  investigating  the  use  of  renormalization  schemes  which 
would  ultimately  lead  to  a  solution  containing  non-integer  multiples  the  gravity 
wave  frequency.  To  avoid  the  secularity,  for  the  present  paper  we  have  restricted 
attention  to  the  vortically  stationary  wave  model  (it  is  this  model. which  has  been 
used  in  the  Figures  2-5);  however,  the  relative  magnitudes  of  the  higher-order  terms 
arc  comparable  in  the  two  cases,  so  that  the  results  would  not  be  greatly  altered 
had  we  used  a  renormalized  solution  for  the  vertically  travelling  wave  model. 
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To  ensure  the  linearity  of  the  gravity  wave,  the  envelope  of  the  amplitude  of 
the  vertically  s^tionary  velocity  field  has  been  normalized  to  induce  at  most  a  15% 
major  species  density  fluctuation  at  the  profile  peak.  In  the  case  of  the  wave  used 
in  Figures  2-4,  this  results  in  an  actual  major  species  density  fluctuation  of  only 
5%  in  the  major  species  density  at  the  profile  peak. 

The  gravity  wave  used  in  Figures  2-4  has  a  horizontal  phase  trace  speed  of 
150  m/s,  an  intrinsic  period  of  2  hrs,  a  horizontal  wavelength  of  iOSO  km,  and  a 
vertical  wavelength  of  55  km.  The  results  for  this  particular  gravif  y  wave  eure  not 
imique,  and  we  have  used  it  solely  for  illustrative  purposes.  In  ail  three  figures, 
we  have  plotted  the  major  species  density  fluctuation,  the  first -oider  minor  species 
response,  and  the  sum  of  the  first  three  orders  in  the  minor  species  response. 

In  Figure  2,  the  background  minor  species  profile  is  that  of  Figure  la.  Figure  2a 
shows  the  response  at  87  km,  where  the  background  density  is  orders  of  magnitude 
below  its  peak  value.  The  first-order  response  here  is  greater  than  400%  (the 
major  species  fluctuation  is  barely  visible  on  this  scale).  As  a  result,  the  nonlinear 
response  is  overwhelming.  In  Figure  2b,  the  response  at  89  km  is  shown.  The 
nonlinear  response  at  this  height  remains  considerable,  and  is  responsible  for  the 
cusp  effect  of  sharpening  the  waveform  peak  emd  flattening  the  wavefrom  trough. 
In  Figure  2c,  the  response  at  95  km  is  shown.  At  this  height,  the  response  is  for 
the  most  part  linear. 

In  Figure  3,  the  background  minor  species  profile  is  that  of  Figure  lb.  Figures 
3a,  3b  and  3c  chow  the  response  at  87  kiii,  89  km  and  95  km,  respectively.  The 
nonlinear  effects  for  this  profile  arc  not  as  pronounced  as  they  were  in  the  previous 
figure,  but  they  are  noticeable  nonetheless.  The  cusping  pf  the  waveform  is  evident 
at  87  km.  Also  noticeable  is  the  phase  reversal  of  the  response  at  the  profile  peak. 
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The  response  is  180  degn-es  out  of  phase  with  the  major  species  fluctuation  below 
the  profile  peak,  and  in  phase  above  the  peak. 

In  Figure  4,  the  background  profile  is  that  of  Figure  Ic.  Figures  4a,  4b  and 
4c  show  the  response  at  89  km,  92  km  and  102  km,  respectively.  Even  for  this 
relatively  broad  profile,  there  is  a  noticeable  cusp  effect  at  certain  height  levels. 

Figure  5  serves  to  show  the  importance  of  higher-order  effects  over  a  broad 
range  of  wave  parameters  for  the  response  at  89  km,  using  the  profile  of  Figure 
la.  In  Figure  5a,  we  have  plotted  the  first,  second  and  third-orde**  response  as 
a  function  of  time  for  a  gravity  wave  with  a  horizontal  phase  trace  speed  of  15C 
m/s  and  an  intrinsic  period  of  2  hrs..  It  is  the  sum  of  these  which  appears  in  large 
dashes  in  Figure  2b.  In  Figure  5b  we  have  plotted  the  first,  second  and  third-order 
response  at  a  fixed  time  as  a  function  of  horizontal  phase  trace  speed  (vpht)  front 
35  m/s  up  to  250  m/s  for  an  intrinsic  period  of  2  hrs.  Oscillations  occur  over  the 
range  of  phase  speeds  because  of  the  vertically  stationary  gravity  wave  model  we 
employ.  Figure  5c  is  the  same  as  Figure  5b,  except  the  intrinsic  period  is  now 
20  min.  We  should  point  out  that  it  is  unlikely  that  smaller-scale  waves  (phase 
speeds  less  than  50  m/s,  say)  will  reach  the  amplitude  necessary  to  induce  even  a 
10%  major  species  fluctuation  before  becoming  unstable  or  encoimtering  a  critic^d 
level,  and  so  smaller-scale  waves  should  not  be  expected  to  produce  a  significant 
nonlinear  response.  However,  it  is  clear  that  the  significance  of  the  higher  order 
terms  is  not  sensitive  to  the  particular  gravity  wave  parameters. 

The  above  results  demonstrate  the  significance  of  the  nonlinear  response  for 
minor  atmospltcric  constituents  which  are  distributed  in  narrow  layers.  For  exam¬ 
ple,  there  cjui  be  a  substiuitial  cusp  effect  in  the  response  at  certain  height  levels. 
As  a  consequence,  it  is  to  be  expected  that  the  spectra  obtained  from  minor  species 


144 


fluctuations  at  particular  height  levels  can  contain  peaks  not  only  at  frequencies 
associated  with  the  actual  gravity  wave  field,  but  at  integer-multiple  harmonics  of 
the  gravity  wave  frequencies  eis  well. 
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5.  Conclusions 


We  have  outlined  a  method  which  can  be  used  to  calculate  higher-order,  non¬ 
linear  terms  in  the  local  response  of  minor  atmospheric  species  to ,  linear  gravity 
waves.  The  method  results  in  an  expansion  of  the  minor  species  response  to  any 
order  in  the  linear  gr  'vity  wave.  Dvsc  to  sectilar  terms  arising  when  the  method,  is 
applied  to  vertically  travelling  waves,  for  the  present  we  have  therefore  only  consid¬ 
ered  the  case  of  vertically  stationary  waves.  Higber-order  effects  have  beeii  shown 
to  be  important  over  a  broad  range  of  wave  parameters.  Oinsequently,  particuleir 
care  must  be  taken  in  inferring  from  the  power  spectra  of  minor  species  flucuation.s 
the  wave  field  of  the  atmosphere  as  a  whole.  In  future  work  we  intend  to  extend  the 
method  to  treat  char  ged  minor  species  and  incorporate  the  effects  of  perturbations 
to  the  production  and  loss  chemistry. 
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1  igar*^  Captf  .inr 

I.  a.  Ch^i)mc-i  function  used  in  Figures  2  and  5.  (b)  Gaussian  distri¬ 
bution  useu  f*-)  b^SSA  atomic  ox>g.i.  inodel,  used  in  Figure  4. 

r  gi>re  2.  (a)  ?  order  major  species  fluctiia'/^r  (solid  line),  first-order  minor 
spfj-.K::  response  (iar;:  ;  u  if  bis),  and  -'ru  o’ t  e  fit?  three  orders  of  the  m’r,or 
sp-r'es  response  (sm-Jl  da-'oesi,  at  c  /  )■  .  fo-  \  gravity  wave  with  an  intrinsic 
peri  J  of  2  hrs  and  a  liorizontal  phusr:  f;ace  speed  of  150  m/s.  The  minor  species 
profile  is  that  of  Figure  la.  (b)  Sn'ue  as  (a),  hut  at  80  km.  (c)  Same  as  (a),  but  at 
95  km. 

Figure  3.  (a)  Fi  st-  .irdci  jna/  »•  -  pecics  fluctuation  (solid  line),  first-order  minor 
specie?  response  and  the  sum  of  the  first  three  orders  of  the  minor 

cpccie.,  response  ( 5mall  dashes),  at  87  km,  for  a  gravity  wave  with  an  intrinsic 
period  of  7  hrs  arid  a  horizontal  phase  trace  speed  of  150  m/s.  The  minor  species 
profile  is  that  of  Figure  lb.  (b)  Same  as  (a),  but  at  89  km.  (c)  Same  as  (a),  but  at 
95  km. 

Figure  4.  (a)  First-order  major  species  fluctuation  (solid  line),  first-order  minor 
species  response  (large  dashes),  and  the  sum  of  the  first  three  orders  of  the  minor 
species  response  (small  dashes),  at  89  km,  for  a  gravity  wave  with  an  intrinsic 
period  of  2  hrs  and  a  horizontal  phase  trace  speed  of  150  m/s.  The  minor  species 
profile  is  that  of  Figure  1-'.  (b)  Same  as  (a),  but  at  92  km.  (c)  Same  as  (a),  but  at 
102  km. 

Figure  5.  (a)  First-order  minor  species  i-esponse  (solid  line),  second-order  mi¬ 
nor  sfiecies  resjxm.M?  (large  dushirs),  luid  third-order  minor  species  response  (small 
dashes),  at  , 89  km,  for  a  gravity  wave  with  lui  intrinsic  period  of  2  hrs  and  a  hori¬ 
zontal  phase  trace  speed  of  150  in/s.  The  mi.<or  species  profile  is  that  of  Figure  la. 
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(b)  Same  as  (a),  but  at  a  fixed  time  as  a  function  of  horizontal  phase  trace  speed 
(Vfh,),  for  a  gravity  wave  with  an  intrinsic  period  of  2  hrs.  (c)  Same  as  (b),  but  for 
a  gravity  wave  with  a  period  of  20  min. 


Figure  1.  (a)  -  Chapman  function  used  in  Figures 
and  and  5.  (b)  Gaussian  distribution  used  in 

Figure  3.  (c)  USSA  atomic  oxygen  model,  used 

in  Figure  4. 
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Figure  1,  (a)  -  Chapman  function  used  in  Figures 
and  and  5.  (b)  Gaussian  distribution  used  in 

Figure  3.  (c)  USSA  atomic  oxygen  model,  used 

in  Figure*  4. 
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Figure  2.  (a)  First-order  major  species  fluctuation 
(solid  line),  first-order  minor  species  response 
(large  dashes) ,  and  the  sum  of  the  first  three  orders 
of  the  minor  species  response  (simill  dashes) ,  at 
87  )cn,  for  a  gravity  wave  with  an  intrinsic  period  of 
2  hours  and  a  horizontal  phase  trace  speed  of  ISO  m/s. 
The  minor  species  profile  is  that  of  Figure  la. 

(b)  Same  as  (a),  but  at  89  Icm.  (c)  Same  as  (a) ,  but 
at  95  km. 
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N‘/N.<i.P»i.1  .p«l.l*2*3 


r i Qure  2b 


Fiquro  2.  (j)  t’irsfc-order  major  species  fluctuation 
(solid  line),  first-order  minor  species  response 
(large  dashes) ,  and  the  sum  of  the  first  three  orders 
of  the  minor  species  response  (small  dashes),  at 
37  )tm,  for  a  gravity  wave  with  an  intrinsic  period  of 
2  hours  and  a  horizontal  phase  trace,  speed  of  150  m/s. 
The  minor  sp«;cios  profile  is  that  of,  Figure  la- 
(b)  Same  as  (a)  ,  but  at  8V  )im.  (c)  Same  as  (a)  ,  but 
at  y5  km. 
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Figure  2.  (,a)  First-order  major  species  fluctuation 
(solid  line) ,  first-order  minor  species  response 
(large  dashes) ,  and  the  sum  of  the  first  three 
orders  of  the  minor  species  respionse  (small  dashes), 
at  87  km,  for  a  gravity  wave  with  ah  intrinsic 
period  of  2  hrs  and  a  horizontal  phase  trace  speed  of 
150  m/s.  The  minor  species  profile  is  that  of  Figure 
la.  (b)  Same  as  (a),  but  at  89  km.  (c)  Same  as  (a), 
but  at  95  km. 
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F i gure  3c 


Fi-iuru  3.  (a)  First-order  major  species  fluctuation 
(solid  lirio),  first-order  minor  species  response 
(largo  dastes)  ,  and  the.  sum  of  the  first  three 
orders  of  the  minor  species  response  (small  dashes) , 
at  87  km,  for  a  gravity  wave  with  an  intrinsic 
period  of  2  hrs  and  a  horizontal  phase  trace  speed  of 
150  m/s.  The  minor  species  profile  is  that  of  Figsire 
lb.  (b)  Same  as  (a),  but  at  89  km.  (c)  Same, as  (a) , 
but  at  95  km. 


r • gur  e  4 a 
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riqurp-4.  (a)  First-ordor  maj'jr  sj^cies  fluctuation 
(solid  lin«),  first-order  minor  species  response 
(larcje  dashes),  .ind  the  sum  of  the*  first  three 
orders  of  the  r>irH>r  species  response  (small  dashes), 
at  8'»  )(Ni,  for  a  'jr.ivity  wave  with  an  intrinsic  p«‘riod 
of  2  hrs  .ind  <i  horizontal  those  tr.ice  speed  of  150  m/s. 
Th»>  minor  ai>ecirs  profile  j«  that' of  Kiqure  Ic. 

(b)  S.ime  IS  (a)  ,  but  at  92  Km.  (c)  Same  4s  (a)  ,  but 
at  1C2  km. 
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Fi  ,ut»'  4.  (j)  FxrRt-ordor  Mjor  specius  fluctuation 
(ijlid  liJn),  first-order  Minor  spociea  response 
(l>ircjc  dishes),  and  the  sum  of  the  first  three 
orders  of  the  minor  species  response  (sMail  dashes). 

It  «9  )tm,  for  a  gravity  wave  with  an  intrinsic  period 
of  2  hen  and  i  horizontal.  j)haBO  trace  speed  of  150  m/s. 
The  Mir, or  s}*,cita  profile  is  that  of  Figure  Ic. 

(b)  Sitn*'  ,is  hut  at  92  )tM,  (c)  Same  as  (a),  but 
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Fiquro  4.  (a)  Fir8t-ord«rr  major  sjx'cies  fluctuation 
(solid  line),  first-order  minor  species  response 
(large  dashes) ,  and  the  sum  of  the  first  three 
orders  of  the  minor  species  response  (small  dashes) , 
at  89  Xm,  for  a  gravity  wave  with  an  intrinsic  period 
of  2  hrs  and  a  horizontal  phase  trace  speed  of  150  m/s. 
The  minor  species  profHe  is  that  of  Figure  Ic. 

(b)  Same  as  (a) ,  but  at  92  km.  (c)  Same  as  (a) ,  but 
at  102  km. 
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Fiqurc  5.  (a)  First-order  minor  species  response 
(solid  liiic"),  second-order  minor  species  response 
(large  dashes),  and  third-order  minor  species 
response  (small  dashes),  at  89  km,  for  a  gravity 
wave  with  an  intrinsic  period  of  2  hrs  and  a 
horizont'il  phase  trace  s|^»ed  of  150  m/s.  The 
minor  s] ecies  prof ilo  is  that  of  Figure  la. 

(b)  Same  as  (a),  but  at  a  fixed  time  as  a  function 
of  horizi  iit.il  phase  trace  speed  (Vp^^^) ,  for  a  gravity 
wave  with  an  intrinsic  period  of  2  hrs.  (c)  Same  as 
(h) ,  but  for  a  gravity  wave  with  a  period  of  20  min. 
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Figure  5.  (a)  Fir*t-order  minor  species  res^nse 
(solifl  line) ,  second-order  minor  species  response 
(large  dashes),  and  third-order  minor  species  response 
(small  dashes),  at  89  )ui,  for  a  gravity  wave  with  an 
intrinsic  period  of  2  hrs  and  a  horizontal  phase  trace 
speed  of  150  m/s.  The  minor  species  profile  is  that 
of  Figure  la.  (b)  Same  as  (a) ,  but  at  a  fixed  time  as 
'V  function  of  horizontal  phase  trace  speed  (Vph,^) ,  for 
a  gravity  wave  with  an  intrinsic  period  of  2  nrs.  (c) 
Same  as  (b) ,  but  for  a  gravity  wave  with  a  period  of 
20  min. 
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FINAL  REPORT 

Submitted  by:  MJ.  Taylor 

Department  of  Physics, 

University  of  Southampton, 

Southampton,  S09  5NH,  U.K 

Period  covered:  3  May  1990  to  4  April  1991. 

1.  This  report  covers  the  extension  (Modification  No  1)  of  the  subcontract  between  the 
University  of  Southampton  and  the  University  of  Cincinnati  .issued  as  part  of  the  US 
Airforce  Contract  F19628-87-K-0023  to  UC  for  research  under  the  project  entitled 
"Investigations  in  Atmospheric  Dynamics". 

2.  The  primary  aim  of  this  research  contract  was  to  image  wave  structure  in  the  OH 
nightglow  emission  io  determine  the  frequency  of  occurrence,  geographical  location, 
horizontal  wavelength,  horizontal  velocity  (and  hence  apparent  period)  of  short  period 
(<1  hour)  mesospheric  gravity  waves  present  over  an  ocean  site.  The  measurements 
were  successfully  made  in  late  March  and  early  April,  1990  from  Maui,  Hawaii  as  part 
of  the  ALOHA-90  campaign. 

3.  Due  to  the  success  of  the  ALOHA  campaign  it  was  decided,  at  short  notice,  to  leave 
the  cameras  set  up  on  Haleakala  and  to  make  a  second  set  of  measurements  during  the 
following  new  moon  period  18-28  April.  All  three  cameras  continued  to  function  well 
and  excellent  data  were  obtained  for  comparison  with  the  Grst  set  of  measurements.  The 
additional  funding  awarded  in  Modification  No  1  of  this  subcontract  have  been  used 
towards  the  costs  incurred  during  this  second  set  of  measurements  and  to  cover  the  initial 
data  analysis. 

4.  In  progress  repon  No  1  (submitted  22  June  1990)  details  of  the  instrumentation  were 
given  together  with  a  summary  of  the  measurements.  An  introduction  to  the  image  data 
was  also  presented  at  the  first  ALOHA  workshop  which  was  held  during  the  CEDAR 
meeting,  Boulder,  in  June  1990. 

5.  The  analysis  of  both  data  sets  has  progressed  well.  The  first  detailed  results  were 
presented  at  the  ALOHA  .session  of  the  AGU  meeting,  held  in  San  Francisco,  December 
1990.  During  this  meeting  it  was  decided  to  concentrate  on  two  nights  25  and  31  March 
so  that  the  image  data  could  be  compared  with  the  measurements  obtained  by  the  other 
instruments  (both  airborne  and  ground  based). 

6.  The  results  of  this  study  have  recently  been  submitted  for  publication  in  Geophysical , 
Research  Letters  as  part  of  a  special  section  on  the  ALOHA-90  campaign.  Copies  of 
these  papers  are  appended.  Although  there  is  still  much  data  to  be  analyzed  the 
campaign  was  a  great  success  and  we  have  already  achieved  several  of  our  stated  goals. 


(Dr  MJ.  Tayior) 
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Appendix  1 

NEAR  INFRARED  IMAGING  OF  HYDROXYL  WAV1E  STRUCTURE 
OVTIR  AN  OCEAN  SITE  AT  LOW  LATITUDES 

M  J.  Taylor  and  M  J.  Hill 

Physics  Department,  University  of  Southampton, 
Southampton,  S09  5NH,  U.K. 


Abstract.  Coordinated  observations  of  ’*ave  structure  in  the  near  infrared  hydroxyl 
(OH)  nightglow  emission  have  been  made  from  Maui,  Hawaii  using  a  suite  of  narrow 
angle  and  all-sicy  TV  cameras.  Two  sets  of  data  were  obtained,  the  first  in  conjunction 
\vith  the  ALOHA-90  campaign  and  the  second  during  the  su^'sequent  new  moon  period. 
Well  formed,  short  period  (<20  min)  wave  patterns  of  comparable  morphology,  djmaraics 
and  abundance  to  those  regularly  imaged  from  mid-latitude  mountain  sites  were  detected 
on  several  occasions.  Although  the  Hawaiian  islands  comprise  several  high  volcanic 
peaks,  the  patterns  were  not  consistent  with  gravity  waves  generated  by  the  interaction 
of  strong  winds  with  the  local  island  topography.  This  suggests  that  ether  mid-latitude 
wave  patterns  may  also  not  be  of  mountain  origin.  The  wave  patterns  imaged  during 
ALOHA-90  were  of  significantly  lower  contrast  than  those  detected  later.  This  effect 
may  be  related  to  changes  in  the  characteristics  of  the  middle  atmosphere  that  occur 
shonly  after  the  spring  equinox. 

Introduction 

Photographic  and  low  light  TV  observations  of  wave  structure  in  the  near  infrared  . 
(NIR)  OH  nightglow  emission  (peak  altitude  -87  km,  haifwidth  5-8  km)  have  been  made 
from  sever?.!  mid-latitude  sites  for  almost  two  decades  [Peterson  and  Kieffaber,  1973; 
Moreek  and  Herse,  1977;  Armstrong,  1986;  Taylor  et  al.,  1987].  Measurements  have  also 
been  obtained  at  aurora!  latitudes  [Qairemidi  et  al.,  1985]  and  from  within  the  arctic 
polar  cap  [Taylor  and  Henrilcsen,  1989].  The  data  frequently  show  extensive  wave-like 
patterns  exhibiting  horizontal  wavelengths  shorter  than  100  km. 

To  date  there  is  very  little  information  available  on  mesospheric  wave  motions  that 
occur  at  low  latitudes.  Peterson  [1979]  reported  visual  and  near  infrared  photographic 
observations  from  Hawaii  of  two  small  scale  wave  events  (termed  ripples),  but  no  details 
were  given  of  any  other  wave  patterns  that  may  have  been  detected.  Subsequent  attempts 
to  image  wave  structure  at  low  latitudes  from  Puerto  Rico  (liS^N),  during  the  AIDA  ’89 
campaign  were  restricted  by  poor  observing  conditions. 

These  observations  of  the  NIR  and  visible  nightglow  emissions  from  Hawaii  allow  the 
properties  of  short  period  mesospheric  wave  motions  that  occur  at  low  latitudes  over  an 
oceanic  site  to  be  investigated  and  compared  with  those  regularly  observed  at  mid¬ 
latitudes. 

Observations  and  Results 

Observations  were  made  from  the  Department  of  Energy  (DOE)  building,  Haleakala 
Crater,  Maui  (20.8"N,  156.2*W,  297Qm),  over  two  consecutive  new  moon  periods;  18 
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March-4  April  during  the  ALOHA-90  campaign,  and  18-27  April  to  provide  additional 
data  on  mesospheric  wave  activity  shortly  after  the  spring  equinox.  The  high  altitude  and 
favourable  weather  during  these  periods  provided  excellent  conditions  for  imaging 
structure  in  the  faint  nightglow  emissions. 

Four  cameras  were  operated,  three  of  which  were  broad  band  NIR  OH  imagers;  a 
narrow  field  Isocon  (30®),  a  medium  field  CCD  (55“),  an  all-slcy  Isccon  (I80“),  and  an  all¬ 
sky  "monochromatic  imager"  filtered  to  observe  structure  in  the  OI  (557.7  nm)  nightglow 
emission  (peak  altitude  -95  km).  Further  details  of  the  instrumentation  are  given  in  the 
ALOHA-^  overview  paper  [Gardner,  1991}.  On  several  occasions  during  this  campaign, 
remarkably  coherent  but  faint  NIR  wave  patterns  were  imaged.  Comparison  with 
spectrometer  observations  confirmed  that  these  displays  were  due  primarily  to  variations 
in  the  OH  emission  [Taylor  et  al.,  1991]. 

Figure  1  shows  four  examples  of  wave  structure  imaged  by  the  CCD  camera  at  low 
elevations  (-15®)  during  ALOHA-90.  The  wave  patterns  were  clearest  in  the  CCD  data 
which  had  the  highest  signal-to-noise  ratio.  Measurements  of  the  horizontal  scale  size 
and  motion  of  each  of  these  wave  patterns  are  listed  in  Table  1.  Figure  la  shows  several 
faint  coherent  wave  crests  curving  upwards  from  the  eastern  horizon.  The  image  was 
obtained  on  21  March  at  11.34  UT  and  formed  part  of  an  extensive  display  that  lasted 
for  about  2  hours.  The  pattern  had  an  average  horizontal  wavelength,  of  -25  km  and 
moved  with  a  speed,  Vj^,  of  -47  ms’',  indicating  an  apparent  period,  r,  of  -9  min.  A 
similar  short  period  display  (-12  min)  is  shown  in  Figure  lb.  The  pattern  was  detected 
on  24  March,  also  around  noon  UT,  but  to  the  south  and  west  of  the  site.  Unfortunately 
meteorological  cloud  (seen  in  the  low  elevation  foreground  of  the  picture)  limited  the 
observations  this  night. 

Figure  Ic  shows  another  well  developed  wave  display.  The  data  were  obtained  at  12:16 
UT  on  30  March  looking  northwards.  On  this  occasion  images  of  wave  structure  were 
also  made  in  the  zenith  indicating  an  extensive  spatially  coherent  display  over  much  of 
the  sky  (confirmed  by  the  NIR  all-sky  observations).  Initially  faint  structure  was  detected 
around  08:00  UT  but  it  was  not  until  later  that  the  display  became  as  conspicuous  as 
shown  in  the  figure.  A  similar  enhancement  in  wave  activity  during  the  latter  half  of  the 
night  was  notii  cd  on  several  other  occasions. 

Sometimes  much  smaller  scale  wave  patterns  termed  "ripple  events"  [Peterson^  1979] 
were  detected.  Figure  Id  is  an  example  and  it  can  be  seen  that  the  morphology  of  the 
wave  forms  is  quite  different.  In  particular,  ripples  are  significantly  smaller  in 
geographical  extent  and  have  shorter  lifetimes  (typically  <1  hour)  than  the  displays 
discussed  above.  In  this  example  the  average  separation  of  the  wave  crests  was  14  j  km 
and  the  wave  packet  moved  with  a  uniform  velocity  of  46  ms'^  indicating  an  apparent 
period  of  5.3  min,  close  to  the  local  Brunt-Vaisala  period.  Unlike  the  larger  scale 
motions,  ripples  are  thought  to  be  caused  by  localised  regions  of  strong  wind  shear 
generating  Kclvin-Helmholtz  instabilities  in  the  mesosphere. 

Discussion 

Wave  displays  of  comparable  morphology  and  dynamics  to  the  examples  shown  in 
Figure  1  were  detected  on  many  occasions  during  ALOHA-90,  and  the  subsequent  new 
moon  period,  indicating  a  wealth  of  short  period  mesospheric  wave  structure  over  the 
Hawaiian  islands.  In  particular  the  patterns  shown  in  Figure  la,b,c  exhibit  good 
correspondence  with  the  mid-latitude  wave  patterns  termed  "stripes"  [Taylor  and 
Henriksen,  1989).  The  principal  features  of  this  type  of  display  are: 
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Horizontal  wavelength:  20-70  km. 

Horizontal  velocity:  10-50  ms*'- 
Lateral  extent  >500  km, 

Duration:  0.5-2  hours. 

The  majority  of  the  reports  detailing  measurements  of  mid  and  high  latitude  wave 
structure  are  of  this  type.  Experimental  and  theoretical  investigations  indicate  that  many 
of  the  patterns  are  induced  by  the  passage  of  short  period  (10-20  min)  gravity  waves 
through  the  nightglow  layers  e.g.  Taylor  el  al.  [1987].  As  similar  wave  motions  were  also 
detected  over  the  Hawaiian  islands  this  suggests  that  the  sources  of  short  period  gravity 
waves  are  plentiful  and  can  occur  at  any  latitude. 

Nearly  all  of  the  mid-latitude  observations  were  made  either  from  high  altitude 
mountain  sites  (to  facilitate  the  seeing  conditions)  or  in  the  vicinity  (<500  km)  of  large 
mountainous  regions.  Wave  activity  in  the  upper  troposphere  and  lower  stratosphere  is 
considerably  enhanced  over  mountain  regions  compared  with  the  oceans.  Much  of  this 
enhancement  has  been  attributed  to  the  interaction  between  strong  tropospheric  winds 
and  the  surface  topography  generating  gravity  waves  (mainly  in  the  lee  of  the  mountains). 
Schoeberl  [1985]  has  modelled  the  propagation  of  mountain  waves  through  the 
stratosphere  and  mesosphere  up  to  an  altitude  of  70  km.  He  determined  that 
orographically  induced  gravity  waves  with  horizontal  scale  sizes  greater  than  30  km  may 
penetrate  freely  into  the  upper  atmosphere. 

If  mountains  are  a  primary  source  of  the  short  period  wave  motions,  seen  in  the 
nightglow  emissions  then  it  is  reasonable  to  expect  considerably  less  wave  activity  over 
a  remote  oceanic  site  such  as  the  Hawaiian  islands  (which  are  over  3,000  km  from  any 
sizeable  land  mass)  than  over  a  large  mountainous  region  such  as  the  Rockies.  During 
both  observing  periods  NIR  wave  structure  was  imaged  on  a  remarkably  high  percentage 
of  nights  (nearly  every  night,  but  not  all  night).  This  compares  with  typically  30-80%  of 
the  clear  s!iy  observing  time  at  mid-latitudes,  depending  on  site  and  season.  TTiis  suggests 
that  there  are  similar  type  sources  for  the  waves  observed  over  the  Hawaiian  islands  and 
over  large  mountainous  region. 

However,  as  the  Hawaiian  islands  themselves  consist  of  several  high  volcanic  mountains 
then  prevailing  winds  blowing  over  local  island  peaks  (such  as  Haleakala  Crater)  may  be 
a  rich  source  of  gravity  waves.  Balsley  and  Carter,  [1989]  and  Hines,  [1989]  have 
investigated  the  properties  of  mountain  waves  observed  in  the  upper  troposphere  and 
lower  stratosphere.  Whether  these  waves  propagate  significant  amounts  of  energy  into 
the  upper  mesosphere  is  not  proven.  Due  to  the  isolated  lo(^tion  of  the  Hawaiian 
islands  any  gravity  waves  generated  by  topographic  forcing  should  be  distinguishable  from 
waves  of  meteorological  (or  other)  origin  as  they  should  be  observed  primarily  downwind 
of  the  source  mountain(s)  and  should  feature  near  zero  phase  speeds  (since  mountain 
waves  are  stationary),  and  curved  wave  crests  (as  the  wave  patterns  are  viewed  at  a  range 
of  only  a  few  hundred  km).  Although  curved  wave  patterns  were  observed  on  some 
occasions  (Figure  la)  the  majority  of  the  displays  were  highly  linear  and  often  extended 
over  much  of  the  sky  including  the  zenith  (Figure  Ic).  Furthermore,  all  of  the  displays 
examined  so  far  have  exhibited  substantial  horizontal  motion  (see  Table  1)  and  no 
preferred  direction  of  propagation  indicating  other  sources  for  the  waves. 

Compared  with  data  recently  obtained  using  the  same  instruments  at  a  mountain  site 
near  Nederland,  Colorado  (40.(rN,  105.6‘’W)  the  wave  patterns  imaged  during  ALOHA- 
90  were  of  significantly  lower  contrast  (<10%)  and  proved  quite  difficult  to  detect. 
During  the  subsequent  new  moon  period  the  contrast  of  the  displays  was  considerably 
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higher.  In  the  OI  (557.7  nm)  emission  similar  morphology  wave  patterns  were  imaged 
on  several  occasions  during  this  period,  but  only  once  during  ALOHA-90.  Unfortunately 
the  CCD  camera  was  not  available  for  the  second  observing  interval  and  a  quantitative 
comparison  of  the  intensity  and  contrast  of  the  patterns  can  not  be  made.  However,  the 
marked  improvement  in  contrast  manifested  itself  as  an  apparent  increase  in  wave 
activity  as  detected  by  the  OH  all-sky  camera.  Figure  2'plots  the  number  of  hours  wave 
structure  was  detected  somewhere  within  the  all-sky  field  (-500  km  circular)  for  both  new 
moon  periods.  Due  to  the  good  seeing  conditions  observations  were  possible  on  most 
nights.  Wave  structure  was  detected  on  nearly  every  occasion  but  many  more  high 
contrast  displays  were  observed  after  the  spring  equinox. 

As  well  as  this  temporal  increase,  the  all-sky  data  indicate  wave  forms  over  a  larger 
spatial  area.  This  effect  may  be  directly  related  to  the  higher  contrast  of  the  structures 
or  it  may  indicate  an  increase  in  the  number  of  potential  sources  in  the  vicinity  of  the 
Hawaiian  islands.  For  tropospheric  w^ave  sources  an  alternative  possibility  is  that  the 
background  winds  in  the  stratosphere  and  lower  mesosphere  changed  significantly  in  the 
interval  between  the  two  observing  periods.  The  COSPAR  atmospheric  model  for  20®N 
indicates  a  mean  zonal  wind  component  m  March  around  40  ms'*  at  -60  km  height. 
Depending  on  the  direction  of  propagation,  waves  with  phase  speeds  lower  than  this  prak 
value  will  be  strongly  attenuated  in  the  lower  mesosphere.  However,  in  April  this  wind 
drops  to  less  than  20  ms'*  and  waves  of  lower  phase  speed  (simOar  to  those  regularly 
seen  in  the  nightglow  emissions)  may  propagate  freely  into  the  upircr  atmosphere.  The 
fact  that  wave  structure  was  detected  frequently  in  both  the  OH  and  the  OI  (557.7  nm) 
emissions  during  the  post  equinox  period  supports  the  argument  for  freely  propagating 
wave  motions  possibly  of  meteorological  origin.  However,  during  ALOHA-90  there  is 
also  good  evidence  for  the  in  situ  generation  of  short  period  waves  [Taylor  and  Edwards, 
1991]. 

Summary 

Short  period  (<20  min)  mesospheric  wave  structure  frequently  exists  at  low-latitudes 
over  the  Hawaiian  islands  (at  least  around  the  spring  equinox).  Many  of  the  wave 
patterns  imaged  were  highly  coherent  and  exhibited  similar  morphology  and  dynamics  to 
the  nightglow  displays  seen  at  mid  and  high  latitudes. 

The  search  for  mountain  waves  induced  by  the  interaction  of  surface  winds  with  the 
local  island  topography  is  continuing.  However,  the  evidence  obtained  so  far  suggests 
(but  is  by  no  means  conclusive)  that  gravity  waves  generated  by  this  mechanism  may  not 
be  an  important  source  of  the  short  period  wave  motions  seen  in,  the  upper  atmosphere. 
For  waves  of  meteorological  origin  wind  filtering  in  the  intervening  atmosphere  may  play 
a  critical  role  in  governing  the  visibility  of  wave  pattcita  of  low  phase  speed. 
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List  of  Figures 


Fig.  1.  Four  examples  of  wave  structure  in  the  near  infrared  OH  nightglow  emission 
imaged  by  the  CCD  camera  during  the  ALOHA-90  campaign.  Each  image  has  been  Cat 
fielded  to  enhance  the  visibility  of  the  faint  displays.  Note  the  similar  morphology  of  the 
wave  patterns  except  in  picture  (d)  which  shows  a  small  scale  ripple  event.  Each  image 
was  obtained  using  an  integration  time  of  20s. 

Fig.  2.  Histogram  plots  the  number  of  hours  NIR  wave  structure  was  detected 
somewhere  within  the  all-sky  field  for  (a)  the  ALOHA-90  campaign  and  (b)  the 
subsequent  new  moon  period.  The  solid  areas  indicate  when  structuic  was  detected  over 
the  whole  sky  while  the  shaded  areas  signify  the  presence  of  some  structure  within  the 
camera  field  but  npt  the  amount.  (Note  CCD  data  were  obtained  for  longer  periods  of 
time  than  indicated  in  (a)  due  to  the  higher  sensitivity  of  the  camera.) 

TABLE  1.  Measurements  of  the  average  horizontal  wave  parameters  for  the  four  wave 
patterns  shown  in  Figure  1  (assuming  an  emission  altitude  of  87  km). 
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Figure  2  -  Histogram  plots  the  number  of  hours  NIB  wave  •ttructure 
was  detected  somewhere  within  the  all-sky  field  for  (a)  the 
AtOHA-90  camf>aign  and  (b)  the  subsequent  new  fiKwn  jxiriod.  Thi; 
solid  areas  indicate  when  structure  was  detected  over  the  whole 
sky  while  the  shaded  areas  signify  the  presence  of  some  structure 
within  the  camera  field  but  not  the  amount.  (Note  CCD  data  were 
obtained  for  longer  periods  of  time  than  indicated  in  (a)  due  to 
the  higher  sensitivity  of  the  camera.) 


176 


Appendix  2 

OBSERVATIONS  OF  SHORT  PERIOD  MESOSPHERIC  WAVE  PATTERNS: 

IN  SmJ  OR  TROPOSPHERIC  WAVE  GENERATION? 

M J.  Taylor,  and  R.  Edwards 

Physics  Department,  University  of  Sr..:  ihampton,. 

Southampton,  S09  5NH,  U.K 

Abstract  Near  infrared  images  showing  wave  structure  in  the  hydroxyl  (OH)  nrghtglow 
emission  have  been  obtained  from  Maui,  Hawaii  during  the  j\LOHA-90  campaign. 
Analysis  of  two  nights  during  this  campaign  (25  and  31  March)  indicate  extensive,  highly 
coherent  linear  wave  patterns  of  very  short  apparent  period  (-5  and  ]0  min 
respectively).  Both  displays  exhibited  several  features  characteristic  of  the  in  situ 
breakdown  of  a  large  scale,  long  period,  upper  atmospheric  wave  disturbance.  Data  in 
support  of  this  mechanism  was  found  by  other  ALOHA  instruments  which  detected 
concurrent  long  period  (1-2  hour)  mesospheric  wave  disturbances  on  beth  occasions. 
However,  a  tropospheric  source  for  these  waves  cannot  be.  ruled  out.  At  least  on  25 
March  a  weather  front  occurred  at  -1400  km  range  with  a  favourable  orientation  and 
location.  Although  its  range  was  relatively  large,  back-ground  winds  may  have 
substantially  increased  the  path  length  of  the  waves  through  the  intervening  atmosphere. 

Introduction 

The  presence  of  coherent  wave  structure  in  the  mid-latitude  mesospheric  nightglow 
cmis-sions  is  well  established  (Moreels  and  Hersc,  1977;  Armstrong,  1986;  Taylor  et  aL 
1987].  Detailed  observations  have  shown  that  these  patterns  often  exhibit  properties 
associated  with  shon  period  (10-20  min)  gravity  waves.  During  the  ALOHA-90  campaign 
similar  morphology  wave  forms  were  detected  at  low  latitudes  over  an  ocean  site  [Taylor 
and  Hill,  1991]. 

Many  sources  for  these  waves  have  been  suggested.  They  include  tropospheric 
disturbances  such  as  jet  streams,  storms,  weather  fronts  and  the  interaction  of  strong 
winds  with  large  mountains  (generating  lee-type  gravity  waves).  Several  upper 
atmospheric  sources  have  also  been  postulated;  they  include  the  auroral  electrojet,  strong 
particle  precipitation  and  the  in  situ  breakdown  of  large-scale,  long  period  atmospheric 
waves  such  as  tides.  It  is  often  difncull  to  identify  individual  wave  sources  from  the 
wealth  of  potentiar  disturbances.  For  tropospheric  sources  this  problem  is  further 
complicated  by  the  effects  of  background  winds  in  the  intervening  atmosphere,  which  can 
signiticantly  modify  the  wave  propagation.  Images  of  structure  in  the  nightglow  emissions 
give  unique  2-dimensional  ini^orma’  wji  on  the  horizontal  wavelength  and  apparent  phase 
velocity  of  the  gravity  wave  pc'  urbation.  These  data  together  with  the  location,  shape, 
orientation  and  time  of  O’  .  jrrcnce  of  the  display  have  proved  valuable  for  finding 
possible  wave  sources  [T  iylor  and  Hapgood.  1988].  By  imaging  stniaure  in  the  nightglow 
emissions  from  the  Hawaiian  islands  the  sources  of  short  period  gravity  waves  over  a  lov/ 
latitude  oceanic  site  can  be  investigated. 
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Observations  and  Results 


Observations  were  made  from  the  Department  of  Energy  (DOE)  building,  Haleakala 
Crater,  Maui  (20.8°N,  156.2®W,  2970m).  The  high  altitude  and  favourable  weather  during 
the  campaign  (18  March-4  April)  pro-vided  excellent  conditions  for  imaging  structure  in 
the  nightglow  emissions.  The  measurements  reported  here  were  made  with  a  medium 
field  (55®)  CCD  camera  filtered  to  observe  the  near  infra-red  (NIR)  OH  nightglow 
emission  (wavelength  range  780-1000  nm).  Further  details  of  the  instrumentation  are 
given  in  the  ALOHA-90  overview  paper  [Gardner,  1991].  Wave  patterns  were  imaged 
on  several  occasions  during  the  campaign  [Taylor  and  Hill,  1991].  In  this  report  the 
properties  of  two  faint  wave  displays  imaged  on  25  March  and  37  March  are  examined 
and  compared  with  simultaneous  airborne  and  ground  based  measurements. 

25  March  Display 

Observations  were  made  from  06:00  LIT  until  dawn  twilight  at  15:40  UT.  A  well 
formed  set  of  steeply  inclined  wav:-like  structures  was  detected  at  low  elevation  to  the 
northeast.  Figure  1  shows  a  CCD  image  of  the  NIR  display  at  07:10  UT.  Several  faint 
forms  emanating  from  the  northeastern  horizon  are  evident.  The  stniaures  arc  clearest 
on  the  right  of  the  image  but  the  data  shows  that  similarly  orientated  wave  crests  exist 
across  the  whole  field.  The  uniformity  of  the  wave  panem  is  quite  remarkable.  The 
display  was  clearest  around  06:30-08:00  UT,  but  was  still  evident  at  11:00  UT  nearly  5 
hours  lat-r. 

A  similarly  shaped  but  fainter  wave  pattern  was  also  detected  to  the  southwest  of  the 
site.  Figure  2  shows  a  composite  ground  map  of  the  nightglow  display  imaged  to  the 
northeast  at  07:10  UT  (Figure  1)  and  an  hour  later  to  the  southwest.  The  flight  path  of 
the  aircraft  during  this  night  is  also  marked  on  the  figure.  Each  image  was  calibrated 
using  k.nown  stars  as  reference  points  [Hapgood  and  Taylor,  1982]  and  mapped  onto  the 
ground  assuming  an  emission  altitude  of  87  km.  The  geographical  location  and 
orientation  of  the  two  wave  sets  indicates  an  extensive,  spatially  and  temporally  coherent 
display .  containing  many  wave  crests  >  800  km  in  length.  To  the  south,  in  the  direction 
of  the  airborne  measurements,  there  was  no  obvious  wave  structure,  but  this  direction 
was  affected  by  light  from  the  Milky  Way.  The  average  separation  of  the  wave  crests, 
determined  from  several  images,  gave  a  horizontal  wavelength  (A,,)  of  15.2  at  0.6  km. 
The  displacement  of  the  wave  crests  (measured  over  an  interval  of  90  min)  yielded  a 
mean  velocity  (v^)  of  51  a  5  ms'*  towards  the  southeast,  indicating  an  unusually  short 
apparent  wave  period  (r)  of  5.0  ±  0.5  minutes. 

31  Vfcrch  Display 

During  this  night  wave  structure  was  detected  a;  most  azimuths  but  was  clearest  to  the 
west  and  north.  To  ihe  west  the  wave  p.attcrn  consisted  of  several  steeply  inclined  crests 
of  similar  morphology  to  those  observed  on  25  March  (Figure  1)  To  the  north  the  wave 
crests  appeared  near  horizontal.  Figure  3  shows  a  composite  ground  map  of  the 
structures  at  09:36  UT  (to  the  west)  and  18  minutes  later  to  the  north.  The  coherence 
of  the  structures  is  not  as  marked  as  those  imaged  on  25  March  (Figure  2). 
.Nevertheless,  the  map  clearly  shows  two  sections  of  an  extensive  wave  display  viewed  at 
orthogonal  azimuths.  A  series  of  images  obtained  at  elevation  angles  ranging  from  17* 
to  90'  confirmed  the  existence  of  coherent  .wave  crests  extending  from  the  northern  and 
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western  horizon  tiirough  the  Icxal  zenith.  The<  pattern  was  present  all  night  indicating 
a  duration  of  over  6  hours.  Measurements  of  the  horizontal  parameters  of  the  wave 
pattern  during  the  interval  06:30-08:00  UT  gave  X,  =  19.8  z:  0.8  km,  v^j  =  32  i  4  ms*' 
towards  the  north  signifying  an  apparent  wave  period  r  =  10.3  i  1.4  min. 

Discussion 

Both  wave  patterns  were  extensive  (>250,000  km^),  remarkably  coherent,  and  lasted 
for  several  hours.  Their  horizontal  wavelengths  and  linear  morphology  were  also 
comparable,  indicating  similar  sources  for  the  di.sturbances.  However,  the  apparent  (i,e. 
observed)  wave  periods  were  significantly  different  at  5.0  and  10.3  min  on  25  and  31 
March  respectively.  In  the  presence  of  a  background  wind  the  intrinsic  period  of  the 
gravity  wave  is  Doppler  shifted  with  respect  to  its  apparent  period: 

^  “))  (') 

where  u  is  the  horizontal  background  wind  component  in  the  direction  of  motion  of  the 
wave,  T  is  the  apparent  period  as  measured  by  an  observer  on  the  ground  and  tj  is  the 
intrinsic  (or  true)  period  of  the  wave  measured  in  a  frame  of  reference  moving  with  the 
background  wind.  The  intrinsic  period  of  the  wave  motion  may  therefore  be  significantly 
different  from  the  apparent  period  of  the  wave  pattern. 

Tfic  motion  of  the  gravity  wave  through  the  atmosphere  depends  on  its  intrinsic  period 
[Hines,  1960]  which  is  usually  indeterminate.  However,  as  the  apparent  horizontal 
velocity  of  the  wave  will  be  independent  of  the  background  wind  regardless  of  height 
(Hines,  1968],  measurements  of  X,  and  of  the  nightglow  wave  patterns  (and  hence  r) 
are  very  useful  for  studies  of  the  wave  sources.  On  the  other  hand,  if  the  waves  are 
generated  by  a  disturbance  that  is  moving  with  respect  to  the  observer,  the  apparent 
period  of  the  wave  pattern  will  vary  with  Q  at  the  height  of  the  source. 

Tropospheric  Source 

The  extensive  nature  and  linear  form  of  both  OH  displays  would  suggest  cither  a  point¬ 
like  source  such  as  a  thunderstorm  [Taylor  and  Hapgood,  1988]  but  at  a  range  of  at  least 
1000  km  from  the  observing  site,  or  an  extended  "line  type"  source  such  as  a  weather 
front  or  a  jet  stream  close  to  the  observing  site  [Hines,  1968].  Without  a  measurement 
of  the  background  wind  in  the  intervening  atmosphere  it  is  difficult  to  estimate  the  range 
over  which  the  wave  energy  can  propagate.  Ho'^'cvcr,  gravity  waves  generated  nearby 
by  strong  winds  blowing  over  local  island  peaks  (such  as  Haleakala  Crater)  can  be 
effectively  ruled  out  on  both  occasions  because  of  the  large  lateral  extent  of  the  displays 
and  their  close  proximity  and  relative  motion  with  respect  to  the  Hawaiian  islands 
(Figures  2,3). 

Surface  weather  charts  for  24  and  25  March  indicate  an  extensive  high  pressure  area 
to  the  northeast  of  Hawaii  with  a  well  developed  low  to  the  northwest  both  at  a  large 
range  (>2000  km).  An  extended  cold  front,  associated  with  these  pressure  centres, 
moved  steadily  towards  the  southeast,  and  came  to  a  halt  between  00:00  and  12:00  UT 
on  25  March  at  a  distance  of  -1400  km  from  the  observing  site.  Both  the  orientation 
and  location  of  the  front  agree  well  with  the  observed  wave  motion  (Figure  2)  but  its 
range  was  relatively  large. 


In  the  absence  of  background  winds  wave  energy  from  a  tropospheric  source  of 
apparent  period  -  5  min  would  reach  the  OH  layer  at  a  quite  short  ground  range  of 
<100  km  (assuming  the  vertical  extent  of  the  source  region  was  no  greater  than  the 
atmospheric  scale  height)  [Hines,  1967],  If  the  background  winds  in  the  intervening 
atmosphere  substantially  increased  the  intrinsic  period  of  the  waves  on  this  occasion  they 
could  travel  very'  large  horizontal  distaners.  This  situation  would  arise  if  u  -  v^^  (-  50 
ms’’)  over  an  extended  height  inter^'a!  between  the  source  and  the  region  of  observation 
(equation  1). 

On  30  and  31  March  surface  weather  charts  show  a  cold  front  associated  with  a 
developing  low  and  a  high  pressure  region  to  the  north  of  Hawaii  at  a  range  of  >1500 
km.  Although  aligned  in  a  direction  similar  to  that  of  the  OH  wave  pattern  its  location 
is  not  consistent  with  the  nonhward  wave  motion  observed  in  the  nightglow  data.  To  the 
south  of  Hawaii  there  were  no  obvious  weather  sources  within  30“  of  latitude. 

In  Situ  Source 

If  the  short  period  waves  were  generated  in  the  upper  atmosphere  in  the  vicinity  of  the 
OH  layer  then  they  would  most  probably  have  been  vertically  evanescent-  The  long 
lifetimes  and  the  large  geographic  extent  of  the  wave  patterns  also  suggests  that  the 
source  regions  were  extensive  and  long  lasting.  One  possible  mechanism  that  could 
account  for  the  observations  has  been  postulated  by  Tuan  et  al.  [1979J.  Under  certain 
conditions  long  period  waves  propagating  through  the  upper  atmosphere  may  partially 
break,  generating  coherent,  small-scale  waves  of  intrinsic  period  close  to  the  local  Brunt- 
Vaisala  period.  On  both  25  and  31  March  was  found  to  be  about  4.6  min  (using 
spectrometer  data  of  Turnbull  and  Lowe,  [1991]).  If  both  wave  patterns  were  generated 
with  an  intrinsic  period  of  this  value  then  the  difference  in  the  apparent  wave  periods 
may  have  been  the  result  of  differing  local  wind  conditions  at  the  source  height  on  each 
night.  The  apparent  period  on  25  March  was  very  close  to  Tj,  indicating  a  background 
wind  component  u  -  3  ms  *  in  the  direction  of  motion  of  the  wave  (equation  1). 
However,  on  31  March  a  neutral  wind  u  -  38  ms**  would  have  been  needed  to  account 
for  an  apparent  period  of  - 10  min.  Neutral  winds  of  this  magnitude  and  variability  are 
frequently  observed  at  nightglow  altitudes  but  unfortunately  there  were  no  coincident 
measurements  over  Maui  with  which  to  compare  these  observations. 

If  the  short  period  waves  were  generated  in  situ  by  dm  mechanism  then  their 
occurrence,  speed  and  direction  of  motion  should  be  associated  with  the  presence  of  a 
large  scale  wave  disturbance  on  each  oeqasion.  Evidence  for  the  presence  of  long  period 
wave  disturbances  on  both  occasions  was  found  in  the  n^suremenis  of  the  other 
ALOHA-90  instruments.  On  25  March  the  airborne  observations  were  made  to  the 
south  of  Maui  (Figure  2).  Stxiium  lidar  measurements  were  obtained  from  08:50  to  15:00 
UT  (mainly  outside  the  limits  of  the  camera  field  of  view).  However,  close  to  the 
equator  a  strong  disturbance  of  horizontal  wavelength  400-550  km  was  measured  (CS. 
Gardner,  personal  communication,  1991).  Zenith  spectrometer  measurements  from 
Haleakala  show  a  near  linear  decrease  in  OH  (3,1)  band  intensity  during  most  of  the 
night  with  a  small  wave-like  disturbance  of  quasi-period  -1.2  hour  [Turnbull  and  Lowe. 
1991).  Coincident  observations  by  the  Aero.spacc  group  using  a  zenith  pointing  CCD 
imager  show  similar  period,  low  amplitude  variations  in  the  OH  (6,2)  band  emission  and 
the  Oi  (0.1)  Atmospheric  band.  The  disturbance  moved  towards  the  southeast  at  a 
speed  of  50-100  ms  ’,  in  a  direction  comparable  with  that  qf  the  short  period  wave 
pattern  (Hecht  and  Walterscheid,  1991]. 
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On  31  March  the  OH  emission  was  brighter  and  the  spectrometer  measured  large 
amplitude  wave-like  disturbance  with  dominant  periods  2.2  and  1.2  hours  in  both  intensity 
and  temperature.  The  Aerospace  images  showed  that  the  disturbance  was  moving  at  a 
velocity  of  30-50  ms"’  similar  in  magnitude  to  that  of  the  short  period  wave  pattern  but 
in  a  general  eastward  direction.  The  aircraft  flew  almost  due  west  (Figure  3)  and 
measurements  over  the  period  08:00-15:30  UT  show  peaks  in  the  horizontal  spectra  with 
wavelengths  of  170  and  340  km  again  indicating  a  large  scale  wave  motion.  However,  as 
the  aircraft  flight  path  was  almost  parallel  to  the  OH  wave  pattern  the  true  horizontal 
wavelength  may  have  been  much  shorter. 

Summary 

Observations  of  extensive  (>250,000  km^)  linear  wave  patterns  with  short  apparent 
periods  (<10  min)  and  long  lifetimes  (several  hours)  are  most  unusual.  Whether  these 
displays  are  peculiar  to  low  latitude  or  oceanic  sites  remains  unknown.  On  25  March 
there  was  a  favourably  orientated  weather  front  where  linearly  extensive  short  period 
(<  10  min)  gravity  waves  may  have  been  generated.  However,  if  the  waves  originated  in 
this  region  then  background  winds  in  the  intervening  atmosphere  must  have  played  an 
important  role  in  considerably  increasing  their  horizontal  range.  No  obvious  source  was 
found  for  the  31  March  display. 

It  appears  more  likely  that  the  waves  reported  here  were  generated  in  the  upper 
atmosphere  in  the  vicinity  of  the  OH  layer.  Both  patterns  exhibited  many  of  the 
properties  characteristic  of  the  sustained  breakdown  of  a  large  scale  upper  atmospheric 
wave  disturbance  (possibly  ol  tidal  origin).  Supporting  evidence  showing  concurrent  long 
period  (1-2  hour)  mesospheric  wave  disturbances  was  found  on  both  occasions.  Several 
other  short  period  (<10  min)  displays  were  also  imaged  during  the  ALOHA-90  campaign 
but  it  is  not  yet  known  if  they  were  associated  with  long  period  upper  atmosphere  wave 
disturbances.  Measurements  of  gravity  waves  at  low  latitudes  are  few  and  these  data  may 
offer  new  evidence  for  the  in  situ  generation  of  short  period  waves  in  the  upper 
atmosphere.  However,  the  possibility  that  the  waves  were  generated  by  other  weather 
disturbances  such  as  jet  streams  [sec  Hines,  1968]  or  distant  thunderstorms  [Taylor  and 
Hapgood,  1988]  has  yet  to  be  investigated. 

Acknowledgements.  We  arc  most  grateful  to  the  director  J,  Drydcn  and  staff  of  Holmes 
and  Narver  Inc.  fpr  arranging  the  use  of  the  DOE  facility  on  Mt.  Haleakala  for  the 
measurements.  We  thank  M  J.  Hill  for  his  considerable  help  both  prior  to  and  during 
the  campaign  and  B.S.  Lanchester  for  many  asetu!  di.scussions.  We  acknowledge  our 
ALOHA-90  colleagues  for  the  use  of  their  data.  Funding  for  this  project  was  provided 
through  the  University  of  Cincinnati  as  part  of  the  U.S.  Air  Force  Office  of  Scientific 
Research  MAPSTAR  programme. 


181 


References 


Armstrong,  E.B.,  Irregularities  in  the  80-100  km  region:  A  photographic  approach.  Radio 
ScK,  21,  313-318,  1986. 

Gardner,  C.S.,  Introduction  to  ALOHA-90:  the  airborne  Hdar  and  observations  of  the 
Hawaiian  airglow  campaign,  Geophvs.  Res.  Lett.,  this  issue,  1991. 

Hapgood,  M.A.,  and  M.J.  Taylor,  Analysis  of  airglow  image  data,  Ann.  Geophvs..  38, 805- 
813,  1982. 

Hecht,  J.H.,  and  R.L  Walterscheid,  Observations  of  the  OH  Meinel  (6,2)  and  O2 
Atmospheric  (0,1)  nightglow  emissions  from  Maui  during  the  ALOHA-90  campaign, 
Geophvs.  Res.  Lett.,  this  issue,  1991. 

Hines,  CO.,  Internal  atmospheric  gravity  waves  at  ionospheric  heights.  Can.  J.  Phvs..  3§, 
1441-1481,  1960. 

Hines,  C.O.,  On  the  nature  of  travelling  ionospheric  disturbances  launched  by  low- 
altitude  nuclear  explosions,  J.  Geophvs.  Res..  72,  1877-1882,  1967. 

Hines,  C.O.,  A  possible  source  of  waves  in  noctilucent  clouds,  J.  Atmos.  5ci..  25,  937-942, 
1968.  . 

Moreels,  G.,  and  M.  Herse,  Photographic  evidence  of  waves  around  the  85  km  level, 
Planet.  Snace  Sci..  21<  265-272,  1977. 

Taylor,  M.J.,  M.A.  Hapgood,  and  P.  Rothwell,  Observations  of  gravity  wave  propagation 
in  the  01  (557.7  nm),  Na  (589.2  nm)  and  the  near  infrared  OH  nightglow  emission, 
Planet.  Snace  Sci..  35.  413-427,  1987. 

Taylor,  M.J.,  and  M.A.  Hapgood,  Identification  of  a  thunderstorm  as  a  source  of  short 
period  gravity  waves  in  the  upper  atmospheric  nightglow  emissions,  Planet.  Space  Sci.. 
36.  975-985,  1988. 

Taylor,  M.J.,  and  M.J.  Hill,  Near  infrared  imaging  of  hydroxyl  wave  structure  over  an 
ocean  site  at  low  latitudes,  Geophvs.  Res.  Lett.,  this  issue,  1991. 

Tuan.  T.F..  R.  Hedinger,  S.M.  Silverman,  and  M.  Okuda,  On  gravity  wave  induced  Brunt- 
Vaisalii  oscillations,  J.  Geophvs.  Res..  84,  393-398,  1979. 

Turnbull,  D.N.,  and  R.P.  Lowe,  Temporal  variations  in  the  hydroxyl  nightglow  observed 
during  ALOHA-90,  Geonhvs.  Res.  Lett.,  this  issue,  1991. 

(Received  April  5,  1991; 

revised  June  5,  1991; 
accepted  June  10,  1991.) 

Copyright  1991  by  the  American  Geophysical  Union. 


102 


list  of  Figures 


Fig.l  CQ?  image  showing  NIR  wave  structure  imaged  to  the  N£  on  25  March  at  07:10 
UT.  The  camera  was  aimed  at  an  elevation  of  20”  and  had  a  field  of  view  cf  55”  by  30”. 

Fig  J  Compxjsite  map  showing  the  scale  size  and  geographic  location  of  the  OH  wave 
forms  imaged  to  the  NE  and  SW  of  Maui  on  25  March,  assuming  an  emission  altitude 
of  87  km  (only  the  brightest  wave  forms  have  been  plotted).  The  x  marks  the  aircraft 
location  at  09:10  UT  on  the  outward  leg  and  15:30  UT  on  the  return  leg. 

Fig3  Composite  map  showing  the  OH  wave  display  of  31  March.  For  comparison  the 
data  are  plotted  on  the  same  scale  as  lugure  2.  The  aircraft  flew  due  west  almost 
parallel  to  the  wave  structures. 
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F iquro  1  -  CCD  image  showing  NIR  wave  structure  imaged  to  the  NE  on  25  March 
at  07:10  UT.  The  camera  was  aimed  at  an  elevation  of  20°  and  had  a, field  of 
vi.jw  of  55°  by  30°.  '  ■ 


Figure  3  -  Composite  map  showing  the  OH  wave  display  of  31  March 
comparison  the  data  are  plotted  on  the  same  scale  as  Figure  2.  ' 
aircraft  flew  due  west  alnwst  parallel  to  the  wave  structures. 
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COINCIDENT  IMAGING  AND  SPECTROMEmiC  OBSERVATIONS 
OF  ZENITH  OH  NIGHTGLOW  STRUCTURE 
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Abstract.  During  the  ALOHA-90  campaign  a  novel  comparative  study  was  made 
between  near  infrared  wave  structure  imaged  in  the  zenith  using  a  CCD  camera  and  that 
detected  at  infrared  wavelengths  by  a  Fourier  Transform  Spectmneter.  Coincident 
measurements  were  made  briefly  on  several  occasions  and  for  an  extended  period  on  31 
March.  The  temporal  variations  imaged  in  the  near  infrared  structure  dunng  this  night 
almost  completely  matched  those  detected  in  the  OH  (3,1)  band  spectrometer  data  when 
similar  viewing  fields  were  compared.  However,  the  image  data  also  displayed  small 
.scale  wave  forms  that  were  not  resolved  by  the  larger  field  instrument  These  structures 
exhibited  significant  changes  in  brightness  and  position  on.  a  time  scale  much  shorter  than 
the  local  Brunt- Vaisala  period  indicating  that  very  high  resolution  measurements  arc 
necessary  to  investigate  short  period  (<20  min)  upper  atmospheric  wave  motions. 


Introduction 

For  many  years  spectrometric  measurements  of  the  night  sky  have  been  made  to 
ascertain  the  nature  and  origin  of  the  upper  atmospheric  night^ow  emissions.  More 
recent  investigations  have  been  into  quasi-periodic  variations  which  often  exhibit 
characteristics  of  freely  propagating  internal  gravity  waves  [Noxon,  1978;  Viereck  and 
Deehr,  1989;  Swenson  et  al.,  1990].  Observations  of  the  bright  infraiia]  OH  Meinel  band 
emissions,  which  originate  from  a  well  defined  layer  centred  at  -87  km,  have  proved  most 
useful  in  these  studies?  As  the  lower  rotational  levels  of  the  OH  emission  are  usually  in 
thermodynamic  equilibrium  with  the  local  atmosphere,  measurements  of  both  the 
emission  intensity  and  its  rotational  temperature  can  be  used  as  a  tracer  of  the  wave 
motion.  Furthermore,  long  exposure  photographic  and  low  light  videoimages  of  the  near 
infrared  (NIR)  OH  emission  can  be  used  to  investigate  the  2-dimensiona]  horizontal 
properties  of  the  wave  disturbance  [Moreels  and  Herse,  1977;  Armstrong,  1986;  Taylor 
et  al.,  1987]. 

In  general,  spectral  observations  suffer  fiom  difficulties  in  distinguishing  between  the 
induced  temporal  and  spatial  variations.  Recently  Taylor  et  al.  [1991]  have  reported  the 
results  of  a  simultaneous  imaging  and  interferometric  study  of  short  period  (-14  min) 
wave  structure  in  the  OH  nightglow  emission.  The  measurements  were  made  at  low. 
elevations,  -15*,.  (to  aid  the  video  observations),  and  show  h  detafl  the  relationship 
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between  a  coherent  gravity  wave  pattern  and  the  induced  intensity  and  temperature 
perturbations.  During  ALOHA-90  a  similar  comparison  wa;'  made.  Here  we  report  novel 
observations  of  a  small  scale  wave  pattern  imaged  in  the  zen  th  and  compare  the  intensity 
records  with  those  detected  by  the  spectrometer. 

Observations 

The  instrumc  .;ts  used  were  a  CCD  camera  (University  of  Southampton)  and  a  Fourier 
transform  spe  'trometer  (The  University  of  Western  Ontario),  both  of  which  were  located 
on  Haleakala  Crater,  Maui  (20.8°N,  156.2“W,  2970m).  A,  short  description  of  the 
instruments  and  their  operational  ch.aracteristics  is  given  by  Gardner  [1991].  The  CCD 
camera  imaged  structure  in  the  submicron  nightgiow  emissions  (half  bandwidth  780~1(KX) 
nm)  while  the  spectrometer  measured  the  zenitl;  intensity  of  several  OH  Meinel  bands 
in  the  wavelength  range  1000-1650  nm. 

Most  of  the  time  the  camera  was  used  to  investigate  the  low  elevation  sky  in  the 
direction  of  the  aircraft  flight  path  (Gardner,  1991].  However,  on  several  occasions  a 
search  for  structure  in  the  zenith  sl^  was  made.  Table  1  lists  the  dates  and  times  of 
these  measurements.  Structure  was  imaged  on  every  occasion  that  zenith  measurements 
were  attempted  (total  of  six  nights)  but  the  observations  were  usually  brief  enough  to 
determine  only  the  scale  size  and  orientation  of  the  wave  forms.  However,  on  31  March 
a  series  of  images  were  recorded  for  detailed  comparison  with  the  spectrometer 
measurements. 

31  March  Disniav 

A  summary  of  the  nightgiow  display  imaged  at  low  elevations  during  this  night  is  given 
by  Taylor  and  Edwards,  [1991].  To  the  west  and  north  of  the  site  a  uniform  wave  pattern 
of  horizontal  wavelength  -20  km  was  detected  over  a  period  of  several  hours  (see  their 
Figure  3).  Images  of  wave  structure  in  the  zenith  were  obtained  around  10:00  UT  and 
again  over  three  hours  later.  Figure  1  is  a  collection  of  six  CCD  images  showing 
examples  of  the  zenith  nightgiow  structure  recorded  at  4  minute  intervals  from  13:14  UT. 
The  pictures  are  quite  exceptional  and  show  considerable  evolution  and  movement  of  the 
structures.  Each  image  has  been  flat  fielded  (to  remove  the  effects  of  lens  vignetting  etc 
that  are  present  in  the  original  records),  and  color  enhanced  to  highlight  the  spatial 
features  of  the  structures  which  exhibited  low  contrast  (s5%).  After  13:18  UT  several 
wave-like  forms  are  evident;  the  east-west  orientation  and  separation  of  these  structures 
is  similar  to  those  observed  in  the  low  elevation  pattern,  conftrming  its  large  geographic 
extent  (>250,000  km‘).  A  total  of  32  images  of  the  zenith  sky  were  recorded  over  the 
period  13:00-14:51  UT.  The  field  of  view  of  the  camera  (55®  by  30")  was  larger  than  the 
spectro.meter  sample  area  (3U°circular),  but  smaller  than  the  Aerospace  Corporation 
camera  (60®  by  40°).  The  variation  in  responsivity  of  the  spectrometer  over  its  viewing 
field  has  not  been  measured. 

Results  and  Discussion 

Images  showing  nightgiow  wave  structure  at  high  elevations  and  in  the  zenith  are  rare 
[Peterson  and  Adams,  1983].  To  date  most  of  tlie  video  and  photographic  measurements 
reported  in  the  literature  have  been  made  at  low  elevation  angles  to  benefit  from  the  2-3 
fold  increase  in  emission  intensity  which  occurs  due  to  line  cf  sight  integration  through 


the  nightglow  layer.  The  oh'ervations  reported  here  vert  tr  ade  using  a  new  bread  band 
nightglow  imager  and  clearly  demoa'trate  the  existent  vt  of  NIR  structure  in  the  zenith 
as  well  as  at  low  elevations.  Although  similarities  in  the  separation  and  orientation  of 
the  wave  forms  are  evident  in  each  image,  the  intensity  distribution  is  viifferent .  indicating 
considerable  evolution  of  the  wave  fielii  in  a  relanv-ily  small  ti.me  iine  val  min). 

By  summing  the  video  signal  in  a  30“  square  centred  on  tl,»;  optic  axis,  the  ;.jteus;ty  of 
the  NIR  TV  data  can  be  compaied  with  the  OH  (3,1)  band  intensity  n^ea-cured  by  the 
spectrometer.  Figure  2a  plots  the  relative  intensity  derived  fi  am  the  snestrometer  data 
(averaged  over  1  min  intervals)  during  the  period  13:00-15:00  UT  [TurnbuH  and  Lov/e, 
1991].  During  the  night  both  the  OH  relative  intensity  and  the  rotation;;  ’  temperature 
exhibited  considerable  wave-like  variations.  In  particular  a  large  scale  perturbation  of 
quasi-period  -2.2  hours  was  detected.  This  wave  motion  is  present  in  the  data  from  the 
larger  field  of  view  Aerospace  camera  [Hecht  and  Walterscheid,  1991].  However,  the 
spectrometer  also  registered  much  smaller  scale  structures  which  are  not  seen  in  the 
Aerospace  data,  possibly  due  to  averaging  effects.  These  variations  are  significant  and 
should  be  present  in  the  image  data. 

Figure  2b  plots  the  relative  intensity  derived  from  the  32  CCD  images.  The  data  were 
recorded  at  intervals  ranging  from  1  to  4  minutes  using  a  fixed  exposure  time  of  30s. 
Each  point  represents  the  average  intensity  of  the  video  signal  over  a  30®  square  field. 
(Note,  the  dashed  line  linking  the  points  serves  only  as  a  guide  to  the  eye.)  It  is 
immediately  apparent  that  the  small  scale  intensity  variations  in  the  spectrometer  data 
are  also  in  the  CCD  data  and  with  no  obvious  time  shift.  The  video  data  were  the  sum 
of  several  signal  components  existing  within  the  broad  pass  band  of  the  camera  (half¬ 
width  780-1000  nm)  and  include  the  OH  (5,1),  (6,2),  (7,3),  (8,4),  (3,0)  and  (9,5)  Meinel 
bands,  the  02(0,1)  Atmospheric  band  near  865  nm,  the  nightglow  continuum  and 
integrated  starlight.  However,  as  the  temporal  variations  in  the  OH  (3,1)  band  intensity 
are  almost  exactly  duplicated  in  the  image  data  this  indicates  that  the  variations  in 
structure  in  the  images  were  primarily  due  to  changes  in  the  OH  emission.  This  is  to  be 
expected  as  the  OH  emission  dominates  this  region  of  the  nightglow  spectrum  having  an 
integrated  emission  intensity  of  typically  12kR. 

Alternatively,  the  intensity  variations  detected  by  both  instruments  could  have  been  due 
to  thin  meteorological  cloud  (although  none  was  visible).  However,  as  the  induced 
intensity  and  rotational  temperature  perturbations  were  highly  correlated  (Turnbull  and 
Lowe,  1991],  with  no  significant  scatter  that  would  result  from  thin  cloud,  there  can  be 
little  doubt  that  the  structures  were  caused  by  changes  in  the  OH  emission.  Any  similar 
perturbations  in  the  NIR  continuum  and/or  the  O2  (0,1)  band  emissions  were  therefore 
comparatively  faint  and/or  were  in  near  phase  with  the  OH  signal.  The  primary  effect 
of  these  emissions  was  to  reduce  the  contrast  of  the  OH  structures.  The  wave  pattern 
imaged  by  the  TV  camera  shows  several  features  of  smaller  dimensions  than  the 
spectrometer  sample  area.  (An  indication  of  the  limits  of  the  circular  spectrometer  field 
is  shown  in  Figure  1  by  the  large  white  box  which  shows  the  30®  square  video  sample 
area:)  Thus,  on  this  occasion  the  spectral  signal  was  the  spatial  average  of  more  than 
one  nightglow  feature.  (Note  the  spectrometer  data  will  be  of  higher  spatial  resolution 
than  indicated  by  this  box  due  to  the  .somewhat  smaller  area  of  the  30®  circular  field.) 
To  investigate  the  temporal  properties  of  the  display  in' more  detail  the  CCD  data  have 
been  re-analyzed  using  a  smaller  sample  area  of  5®  square  (marked  by  the  small  white 
box)  which  corresponds  to  a  zenith  "footprint”  of  -8  km  square.  As  the  wavelength  of 
the  OH  pattern  was  -20  km  individual  bright  and  dark  structures  of  this  dimension 
should  be  rfesolved. 


Figure  3  shows  the  relative  intensity  data  for  this  sample  field.  For  comparison  the 
data  have  been  plotted  using  the  same  intensity  and  time  scales  as  Figtne  2b.  Although 
similar  features  are  present  in  both  plots  the  short  wavelength  structures  are  now  evident 
and  the  overall  contrast  of  the  variations  is  increased.  More  critically,  although  a  4 
minute  sample  interval  (13:00-14:00  UT)  is  ample  for  tracking  the  large  scale 
perturbations  present  in  the  figure,  it  is  too  long  for  registering  accurately  the  small  scale 
variations.  After  14:15  UT  the  sample  interval  was  reduced  to  1  min  and  the  variations 
are  clearly  detected.  At  the  time  of  recording  the  data  little  or  no  change  in  the  images 
were  observed  by  eye.  However,  Figure  3  shows  that  significant  changes  in  the  integrated 
nightglow  brightness  took  place  on  a  minute-by-minute  time  scale. 

Studies  of  wave  structure  in  the  upper  atmosphere  usually  assume  the  local  Brunt- 
Vaisala  (buoyancy)  period,  r^,,  as  a  natural  lower  limit  for  the  temporal  resolution  of  the 
measurements.  At  nightglow  altitudes  is  typically  5-6  min.  Wave  motions  with  phase 
velocities  of  a  few  tens  ms'^  can  propagate  significant  distances  (several  km)  during  this 
time.  These  measurements  demonstrate  the  need  for  high  temporal  and  spatial 
resolution  when  investigating  the  properties  of  short  period  (<20  min)  waves  that 
frequently  occur  in  the  nightglow  emissions. 

Summary 

The  good  seeing  conditions  at  this  site  permitted  OH  nightglow  wave  structure  to  be 
imaged  in  the  zenith.  The  Close  similarity  in  the  relative  intensity  data  measured  by  the 
spectrometer  and  the  CCD  camera  indicates  no  observable  difference  in  the  wave 
structure  arising  primarily  from  high  OH  vibrational  levels  (v'>5)  compared  with  that 
observed  in  the  v'-3  level  of  the  OH  (3,1)  band.  To  date  most  image  data  on  the 
nightglow  emissions  have  been  obtained  at  low  elevations,  and  comparative  studies  with 
other  optical  devices  are  rare  [Taylor  et  al.,  1991].  The  extension  of  high  resolution 
imaging  measurements  into  the  zenith  has  several  advantages  but  is  made  at  the  expense 
of  signal  level. 

During  ALOHA-90  NIR  wave  patterr.s  of  horizontal  wavelength  <  <  100  km  were 
observed  on  several  occasions.  Large  field  instruments  will  naturally  integrate  out  much 
of  this  structure  which  appears  to  be  relatively  common  over  both  mountain  and  oceanic 
sites.  Coincident  imaging  and  spectrometric  measurements  of  the  OH  emission  are 
highly  complementary:  the  spectrometer  establishes  the  spectral  integrity  of  the  emission 
and  gives  accurate  intensity  and  temperature  data  while  the  imager  records  the  2- 
dimensional  spatial  structure  of  the  wave  perturbation,  its  motion  and  temporal  evolution. 
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List  of  Figures 


Fig.l.  Series  of  CCD  images  showing  small  scale  wave  structure  in  the  near  infrared 
zenith  nightglow  emissions.  Each  image  has  been  flat  fir^ded  and  color  enhanced  to 
highlight  the  wave  features.  The  structures  app>ear  mainly  yellow-green  but  brighter 
features  are  red.  The  Camera  was  aimed  at  an  azimuth  of  2^  and  the  near  eust'west 
aligned  waves  are  almost  vertical  in  the  images. 

FigJ.  Relative  intensity  of  (a)  the  OH  (3,1)  band  recorded  by  the  spectrometer  using 
a  30°  circular  field  of  view  and,  (b)  the  broad  band  NIR  structure  derived  from  the  CCD 
images  using  a  30”  square  sample  area. 


FigJ.  Relative  intensity  of  the  CCD  data  determined  using  a  5“  square  sample  area.  For 
comparison  the  data  are  plotted  on  the  same  scales  as  Figure  2.. 

TABLE  1.  Dates  and  times  when  NIR  structure  was  imaged  in  the  zenith  by  the  CCD 
camera. 


DAY  DATE  TIME  (UT) 


081  22/3  15:30 

082  23/3  14:53 

087  28/3  15:16 

089  30/3  12:22  -  12:59 

090  31/3  09:50  -  10:00 

13:00  -  14:51 
12:19  -  12:56 


092 


02/4 


13-1 4  UT 


13-18 


13-30  13-34 


Figure  1  -  Series  of  CCD  images  showing  small  scale  wave  structure,  in  the  near 
infrared  zenith  nightglow  emissions.  Each  image  has  been  flat  fielded  and  color 
enhanced  to  highlight  the  wave  features.  The  structures  appear  mainly  yellow-green 
but  brighter  features  are  red.  The  camera  was  aimed  at  an  azimuth  of  280®,  and 
the  near  east-west  aligned  waves  are  almost  vertical  in  the  images. 
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Figure  3  -  Relative  intensity  of  the  CCD  data 
determined  using  a  5*  square  sample  area.  For 
comparison  the  data  are  plotted  on  the  same  Scales 
at  Figure  2. 
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INTRODUCTION 

This  sub-contract  was  set-up  to  pay  partial  field  costs  (shipping,  travel  and 
living)  of  University  of  Western  Ontario  participation  in  the  ALOHA  '90  campaign  of 
co-ordinate  observations  of  mesopause  region  waves  under  the  direction  df  Professor 
Chet  Gardner,  University  of  Illinois.  Additional  funding  was  provided  by 
MAPSTAR/Utah  State  (data-processing  and  interpretation)  and  the  Institute  for  Space 
and  Terrestrial  Sciences  (salaries  and  additional  travel  and  living). 


SUMMARY  OF  FIELD  ACTIVITIES 

UWOMI-I,  a  Fourier  transform  interferometer  was  shipped  to  Maui  for 
installation  at  the  summit  of  Haleakelaon  March  22,  1991.  It  was  seriously  damaged 
in  shipment  and  required  extensive  repairs  before  being  put  in  operation  on  March  25. 
It  remained  on  the  summit  until  April  11,1 991 .  The  system  computer  was  damaged 
in  a  lightning  storm  on  April  5,  and  attempts  to  repair  it  were  unsuccessful.  Eight 
nights  of  data  were  obtained  consisting  of  over  7000  individual  airglow  spectra. 

D.  N.  Turnbull  (UWO/ISTS)  was  in  charge  of  field  operations  and  remained  on 
site  throughout  the  period  of  operations.  P.  Meier  (UWO/ISTS)  was  present  for  the 
first  week  of  operations  to  aid  in  initial  ser-up  and  service  the  damaged  instruments. 
R.  P.  Lowe  (UWO)  acted  as  observing  assistant  for  the  last  ten  days  of  the  operation. 


RESULTS 

All  data  has  been  processed  and  interpreted.  Three  papers  based  on  the  results 
obtained  combined  with  those  of  other  campaign  participants  will  appear  In  a  special 
issue  of  Geophysical  Research  Letters  vto  be  published  in  late  1991. 


